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DISPLAY3D 


A  Graphics  Preprocessor  for  CHIEF 


INTRODUCTION 

DISPLAY3D  processes  CHIEF  defined  geometries  into  a  format  that  can  be 
easily  displayed  The  computer  program  CHIEF  [1]  is  designed  to  obtain 
F.pproximate  solutions  to  exterior  steady-state  acoustic  radiation  problems  for 
surfaces  of  arbiti aiy  shapes  vibrating  with  a  prescribed  normal  velocity. 
Being  unable  to  view  the  shapes  defined  by  CHIEF  can  lead  to  many  difficult 
debugging  problems  as  well  as  incorrect  solutions.  DISPLAY3D  provides  points 
as  though  the  data  were  locations  on  a  three-dimensional  surface  or  solid  ar  d 
a  description  of  how  these  points  are  connected.  A  graphics  program  commonly 
called  a  graphics  driver,  is  required  to  display  these  data  point'.  Since 
here  are  so  many  graphics  standards  (protocols)  and  graphics  devices 
av_ilable  it  would  be  virtually  impossible  to  develop  a  generalized  graphics 
program;  however  two  examples  of  graphics  drivers  are  provided.  If  these 
drivers  cannot  be  used  on  available  graphics  devices,  the  user  may  write  a 
graphics  program  thct  can  read  DISPLAY3D  output  files,  or  use  one  of  the 
commercial  plotting  packages  available.  The  format  of  the  data  and  an  example 
of  a  graphics  driver  can  b  found  later  in  this  report.  Because  most  research 
laboratories  have  graphics:  capabilities,  the  task  of  taking  the  results  of 
DISPLAY3D  and  plotting  it  on  an  available  graphics  device  usually  requires 
little  effort. 

DISPLAY3D  has  several  opt'nns  that  makes  it  a  versatile  package.  The 
code  permits  the  user  to  change  the  orientation  from  which  the  data  is  viewed. 
DISPLAY3D  also  provides  the  option  of  viewing  the  direction  of  the  normal 
velocity  vector  as  defined  by  CHIEF.  These  options,  as  well  as  others,  will 
be  explained  in  detail. 


CODE  DESCRIPTION 

The  DISPLAY3D  program  has  two  options  for  reading  the  CHIEF  geometry 
data.  The  first  option  uses  the  main  program  PLOTCID  to  read  the  data  from 
the  output  file  CID  (Chief  Interactive  Driver)  [2]  which  is  a  computer  program 
that  interactively  generates  the  control  routine  for  CHIEF.  While  creating  a 
CHIEF  control  program,  CID  captures  the  information  from  the  user’.s 
interactive  session  and  .saves  it  in  a  separate  data  file  which  Is  read  the 
next  time  CID  is  executed-  This  data  file  is  .also  used  for  ''ntion  1  in  th'’ 
DISPLAY3D  package.  The  second  opt  ion  uses  t.he  subroutine  PI  TiTlin  J  KV  wh.ich  's 
placed  in  the  CHIEF  driver  prugr:.iin  after  the  l;>.st  LDSURR  call  ,  Id.OTCI!  I  HI' 
shares  the  same  variables  as  the  CHIEF  di-iver  pr'igra.m  throig’n  the  use  of  a 
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FORTRAN  COMMON  statement.  BoLh  of  these  options  generate  two  output  files 
which  include  all  translations,  rotations,  and  reflections  as  prescribed  in 
the  CHIEF  program. 

The  hierarchy  of  the  two  options  is  given  in  Fig.  1.  The  program  names 
in  quotes  indicate  the  programs  from  which  the  CHIEF  surface  data  are  obtained 
for  each  option.  The  advantage  of  option  1  is  that  the  CHIEF  driver  does  not 
have  to  be  compiled,  linked,  and  executed  in  order  to  create  the  graphing  data 
files.  However,  if  for  some  reason  CID  was  not  executed,  the  second  option  is 

provided. 

DISPLAY3D 


(Option  1) 
"CID" 


(Option  2) 
"CHIEF  DRIV'ER" 


PLOTCI 


LOAD  CID  DATA 


PLOTCHIEF 


RECTANGULAR_PLANE 

CIRCULAR_PLANE  - 

ELLIPTICAL_PLANE  - 
CIRCULAR  CYLINDER 


PUT_P0INT 
Pirr_POINT 
PUT_P0INT 
PUT  POINT 


ELLIPTICAL  CYLINDER - PUT  POINT 


SPHEROID 


PUT  POINT 


PROLATE  OBLATE  SPHEROID  -  PUT  POINT 


QUADRILATERAL  - 
AXISYMMETRIC  — 

TRm'GLE  - 

CONNECT  ARROWS 


PUTPOINT 
PUT_P0INT 
PUT  POINT 


PUT  POINT 


DRAW  ARROWS 


Fig.  1  -  The  program  hierarchy  of  DI3PLAY3D. 
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PLOTCID  (First  Option) 

PLOTCID  is  the  main  program  of  the  computer  code  and  is  designed  as  the 
control  program  for  ten  different  geometry  subroutines.  This  program 
determines  the  necessary  rotations  or  reflections  for  each  surface  in  a  block. 
This  prograa  generates  two  files  that  are  used  for  displaying  the  three- 
dimensional  wire  mesh  of  the  CHIEF  surfaces.  The  names  and  descriptions  of 
the  subroutines  are  given  below. 

LOAD_CID_DATA  Reads  necessary  data  from  a  CID  save 

file  into  variables  used  by  PLOTCID. 

PUT_P0INT  Translates  ?  3D  point  from  a  local 

coordinate  system  to  the  global 
coordinate  system. 

Depending  upon  the  symmetry  condition 
chosen  in  CHIEF,  this  subroutine 
reflects  or  rotates  the  point. 
PUTPOINT  stores  these  points  in  ai 
array  which  is  rotated  by  subroutine 
R0T3D. 


RECTANGULAR_PLANE  Generates  poii  ts  and  lines  for 

drawing  a  rectangular  planar  region. 

CIRCULARPLANE  Generates  points  and  lines  for 

drawing  a  circular  planar  region  or 
part  of  one. 

ELLIPTICALPLANE  Generates  points  and  lines  for 

drawing  an  elliptical  planar  region 
or  part  of  one. 

CIRCULAR_CYLINDER  Generates  points  and  lines  for 

drawing  a  cylinder  or  section  of  a 
cylinder. 

ELLIPTICAL_CYLINDER  Generates  points  and  lines  for 

drawing  an  elliptical  cylinder  or 

section  of  cylinder. 

SPHEROID  Generates  points  and  lines  for 

drawing  a  spheroid  or  .sectl'm  of  a 
spheroid . 
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PROLATEOBLATESPHERuID  Generates  nnints  and  lines  for 

drawing  a  prolate  spheroid  or  section 
of  a  prolate  spheroid.  This  same 
routine  is  also  used  for  generating 
points  and  lines  for  an  oblate 
spheroid . 

TOROID  Generates  points  and  lines  for 

drawing  a  toroid  or  part  of  a  toroid. 

QUADRILATERIAL  Generates  points  and  lines  for 

drawing  these  finite-element-type 
inputs  using  linear  or  quadratic 
interpolation  over  a  quadrilaterial . 

AXISYMMETRIC  Generates  points  and  lines  for 

drawing  these  f inite-element-type 
inputs  using  linear  or  quadratic 
axisymiuetric  interpolation. 

TRIANGLE  Generates  points  and  lines  for 

drawing  these  finite-element-type 
inputs  using  linear  or  quadratic 
interpolation  over  a  triangle. 

R0T3D  Rotates  coordinates  by  a  user 

specified  angle  for  different  viewing 
perspectives  and  stores  these  points 
in  a  data  -file. 


CONNECTARROWS  Generates  points  used  to  describe  the 

connectivity  of  the  reference  axis 
and  associated  labels. 


DRAWARROWS  Generates  the  coordinates  of  the 

points  needed  for  the  reference  axes 
and  associated  labels. 

PLOTCHIEF 

This  is  a  subroutine  that  is  placed  in  the  CHIEF  control  program  after 
the  last  call  to  LDSURR.  PLOTCHIEF  and  PLOTCID  are  identical  after  the 
i.aitial  variables  are  determined. 
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CODE  OPERATION 

Depending  upon  which  option  of  DISPLAY3D  is  used  (see  section  entitled 
CODE  DESCRIPTION)  the  user  must  access  the  appropriate  graphics  code.  As  a 
first  example,  option  2  will  be  exercised  on  a  free-flooded  cylinder  given  in 
sample  run  2  of  the  CHIEF  [1]  users  manual.  Appendix  A  contains  the  CHIEF 
driver  program  with  the  subroutine  call  to  PLOTCMEP  as  follows: 


CALL  PL0TCHIEF(RUNID,NBLKS-,S711IYP.,SUBDIV,AaE',AY,AZ).. 


ThiS'  subroutine  c’all  is  highlighted  with  an  arrow  in  the  left-hand,  margin. 
The  variables  RUNID,  NBLKS-,.  and  SYMTYP  are  all  local  to  the  CHIEF  driver; 
however,  the  user  has  the  ability  to  vary  RUNID  and  NBLKS  independently  of  the 
CHIEF  driver.  For  example,  the  user  might  want  to  set  NBLKS  to  1  in  the 
PLOTCHIEF  call  in  order  to  display  only  the  user-defined  section  regardless  of 
symmetry  type.  If  NBLKS  is  set  to  a  positive  value,  the  direction  of  the 
normal  velocities  (SIGN(IZAX))  will  be  identified  by  dashed  and  solid-lines. 
But  if  NBLKS  is  negative,  only  solid-line  drawings  will  be  produced.  The 
variable  SUBDIV  is  a  user-specified  flag.  When  the  flig  is  set  (SUBDIV=1) , 
PLOTCHIEF  will  honor  the  subdivisions  (NSU  and  NSV)  used  in  the  CHIEF  driver. 
If  the  flag  is  cleared  (SUBDIV=0) ,  PLOTCHIEF  uses  default  values  for  the 
number  of  subdivisions  which  should  be  optimal  for  viewing.  The  user  might 
set  the  SUBDIV  flag  equal  to  one  in  order  to  determine  visually  whether  the 
proper  number  of  subdivisions  was  chosen.  The  variables  AX,  AY,  and  AZ  are 
the  angles  in  degrees  used  to  rotate  the  figure  about  its  X,  Y,  and  Z  axis, 
respectively.  See  Figure  2. 


X2  (Y) 


X3  (Z) 


Figure  2  illustrates  the  direction  of  the  rotation  angles. 
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Since  rotations  are  not  associative,  the  user  must  know  the  order  in  which  the 
rotations  are  performed.  The  first  rotation  is  around  the  object’s  X  axis, 
next  around  its  Y  axis  and  then  about  its  Z  axis.  It  is  important  to  note 
that  after  a  rotation  around  the  X  axis  has  occurred  the  object’s  Y  and  Z  axes 
are  no  longer  the  global  Y  and  Z  axes. 

This  CHIEF  driver  must  be  compiled,  linked,  and  executed  and  once 
executed,  two  data  files  are  created  for  graphing  purposes.  The  extensions  of 
these  files  are  CRD  and  CON,  coordinate  and  connectivity  data,  respectively. 
For  this  example,  the  files  created  are  TEST2.C0N  and  TEST2.CRD.  The  format 
of  the  coordinate  data  is  3(E15.0),  and  the  format  of  the  connectivity  of  the 
points  is  3(112).  The  connectivity  file  contains  an  array  of  three  numbers. 
The  first  two  entries  determine  the  start  and  stop  of  a  line  segment.  The 
third  entry  is  the  color.  The  sign  o'  the  color  entry  determines  if  the  line 
segment  will  be  solid  (positive)  or  dashed  (negative)  .  The  value  of  this 
entry  determines  the  color  of  the  line  segment  depending  on  the 
characteristics  of  the  graphics  output  device.  These  colors  can  be  changed  by 
changing  the  parameters  COLOR_USER_DEFINED  and  C0L0R_TRANSF0RMED  at  the  top  of 
the  PLGTCID  and  PLOTCHIEF  programs.  These  parameters  represent  the  user 
defined  section  of  the  drawing  and  the  transformed  portion.  An  example  of 
both  of  these  data  files  are  listed  in  Appendix  B. 

These  data  files  can  be  plotted  using  almost  any  available  plotting 
package  as  long  as  the  protocol  has  "move  to  a  point*  and  "draw  a  line  from 
one  point  to  another"  commands.  Two  graphics  drivers  are  included  in  Appendix 
C  for  the  user’s  convenience.  These  drivers  provide  the  basic  graphing 
capability.  The  first  driver,  DRAW240,  can  be  used  on  any  VT240  or  compatible 
terminal.  This  program  reads  in  the  coordinate  and  connectivity  data,  scales 
it,  and  displays  the  drawing  on  the  screen.  The  user  can  modify  this  basic 
driver  by  adding  zoom,  rotation,  scaling,  or  more  advanced  color  capability. 
Also,  an  option  to  save  the  data  to  file  rather  than  displaying  the  drawing  on 
the  screen  could  easily  be  added.  The  second  driver  provided  is  written  in 
PY-WAVE  [4]  and  uses  PV-WAVE’s  powerful  graphing  and  data  manipulation 
functions.  This  driver  creates  a  ReGIS  (remote  graphics  instruction  set) 
display  file  for  viewing,  or  a  QMS  (Quality  MicroSystens)  file  for  dumping  to 
a  Talaris  or  QMS  printer. 

Figure  3  is  the  display  of  the  example  described  above  with  a  60,10,0 
viewing  perspective  with  the  negative  normals  illustrated  by  dashed  lines  and 
positive  normals  with  solid  lines.  Figure  3  was  plotted  using  the  PV-WAVE 
driver  on  a  Talaris  800  printer.  The  reference  axes  in  the  lower  left  corner 
of  Fig.  3  are  rotated  in  the  same  manner  as  the  CHIEF  display. 
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Fig.  -  A  free  flooded  cylinder  using  the 
DISPLAY3D  plotting  routine. 


If  the  CID  software  is  usc'^  to  create  the  CHIEF  control  program,  the  user 
can  access  DISPLAY3D  by  implementing  option  1  by  the  command  RUN  PLOTCID. 

The  program  prompts  for  the  name  of  the  file  that  was  created  by  CID  and  the 
rotation  angle  (see  Fig.  4).  The  output  of  PLOTCID  is  identical  to  files 
created  by  PLOTCHIEF  provided  that  the  CHIEF  surfaces  are  identical  for  both 
program  runs. 
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D  ^er  the  filename  used  in  the  CID  run:  TEST2 


DISPLAY: 

1. )  Only  the  user  def.ned  surfaces 

2. )  User  defined  surfaces  with  selected  transformations. 

Enter  1  or  2  (default  2):  2 

THE  NUMBER  OF  SUBDIVISIONS  PLOTTED: 

1. )  The  optimal  number  for  visualization. 

2. )  The  number  defined  in  the  CHIEF  driver. 

Enter  1  or  2  (default  1):  1 


DISTINGUISH  BETWEEN  POSITIVE  AND  NEGATIVE  NORMAL  VELOCITY: 

1. )  Yes.  (solid  lines  for  positive,  dashed  for  negative) 

2 .  )  No.  (only  solid  I ines) 

Enter  1  or  2  (default  2):  1 


Enter  the  angles  (in  degrees)  for  rotating  the  model 
about  the  X,  Y,  and  Z  axes,  (default  0.0,  0.0,  0.0) 

(XROT,  YROT,  ZROT)  :  60,10,0 


Coordinate  data  has  been  stored  in  TEST2.CRD 
Connection  data  hrs  been  stored  in  TEST2.C0N 


Fig.  4  -  An  interactive  run  fro-  * iiOTCID  program. 


A  second  example  is  an  array  of  hemispheres  generated  using  the  CHIEF 
program.  The  display  is  shown  in  Fig.  5.  DISPLAY3D  can  handle  all  of  the 
CHIEF  surfaces  including  the  finite-element-type  inputs;  however,  the  number 
of  subdivisions  in  DISPLAY3D  is  chosen  to  best  view  the  three-dimensional 
surface . 
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Fig.  5  -  An  array  of  hemispheres  drawn  using  DISPLAY3D  program. 


CODE  INSTALLATION 

The  authors  have  been  using  the  DISPLAY3D  program  on  the  VAX  and  MICROVAX 
computers  running  under  the  VMS  5.1  operating  system.  The  source  code  is 
device  independent,  but  any  graphics  driver  is  device  dependent.  A  listing  of 
the  source  code  is  given  in  Appendix  D. 

To  receive  more  information  on  DISPLAY3D  or  the  Fortran  computer  code, 
contact  Tina  Siders  at  P.O.  Box  568337,  Orlando,  FL  32856-8337. 
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APPENDIX  A  -  A  CHIEF  Driver  of  a  Free-Flooded  Cylinder 


This  driver  appears  in  the  CHIEF  user’s  manual  as  example  2. 

C  PROGRAM  TEST2 

C  Book  Examples 

C  VARIABLE  DECLARATION 

C  *****  CONTROL  88  ***** 

C  PROGRAM  CHIEF88  DRIVER 

C  MXSREG  -  MAXIMUM  NUMBER  OF  SURFACE  REGIONS 

C  MXIPS  -  MAXIMUM  NUMBER  OF  INTERIOR  POINTS 

C  MXARS  -  MAXIMUM  NUMBER  OF  SURFACE  SUBDIVISIONS/SYM  BLK 

C  MXGAUS  -  MAXIMUM  ORDER  OF  GAUSSIAN  QUADRATURE 

C  MXQPTS  -  MAXIMUM  NUMBER  OF  QUADRATURE  POINTS/SUBDIVISION 

C  MXBLKS  -  MAXIMUM  NUMBER  OF  SYMMETRY  BLOCKS 

C  MXFFP  -  MAXIMUM  NUMBER  OF  FAR-FIELD  POINTS 

C  MXNFP  -  MAXIMUM  NUMBER  OF  NEAR-FIELD  POINTS 

C  MXPTSC  -  MAXIMUM  NUMBER  OF  POINT  SOURCES 

C  MAXCOR  -  MAXIMUM  NUMBER  OF  FINITE  ELEMENT  NODES 

C  MXFPS  -  MAXO(MXARS-MXIPS, MXFFP, MXNFP) 

PARAMETER  (MXSREG=500) 

PARAMETER  (MAXCOR=1000) 

PARAMETER  (MXIPS=20) 

PARAMETER  (MXARS=500) 

PARAMETER  (MXGAUS=64) 

PARAMETER  (MXQPTS=512) 

PARAMETER  (MXBLKS=100) 

PARAMETER  (MXFFP=361) 

PARAMETER  (MXNFP=361) 

PARAMETER  (MXPTSC=20) 

PARAMETER  (MXFPS=520) 

C  PARAMETER  (MXFPS=MAX0(MXARS-.MXIPS, MXFFP, MXNFP)) 

C  PARAMETER  (N»DVEC=2«MXARS) 

C 

C  INPUT  COMMONS 

C 

COMMON/CONST/RHO,C 
COMMON/PRTCOM/NUNPRT, NUNERR 
COMMON/PRTRD/RUNID,OATE 
CHARACTER»4  RUNID 
CHARACTER'S  DATE 

COMMON/NDASG/NDQPTS , NDPMXS , NDVMXS , NODECM, NOVB-S , NDSPS , 

1  NDPMXF , NDVMXF , NOPMXN , NDVMXN , NOPSSP, NDEXPR , NDCOMV , 

1  NDTEMP.NDZRDB.NUrATB 

COMMON/SVALS/NSREG, NSEQNS (MXSREG) , SUL (MXSREG) ,SUU(MXSREG) , 

1  SVL (MXSREG) ,SVU (MXSREG) ,NSU (MXSREG) ,NSV (MXSREG) , 

1  CCS ( 10 , MXSREG) . TRNSS (3 , MXSREG) , IZAX (MXSREG) , 

1  IORDU(MXSREG),IOROV(UXSREG),NCCEQS 

COMMON/CORD/COQRDS (MAXCOR, 3) 

COMMON/IPTS/NUMIPS,IPXS(3, MXIPS) 

REAL  IPXS 

COMMON/PTSINP/NUMPTS , PTSRCS (4 , MXPTSC) , PTSWT (MXPTSC) , 

1  lOPTSC (MXPTSC) 

COMPLEX  PTSWT 

COMMON 'PLfINP/AINC,THTINC,PHIINC,ISCATR 

COMMON/BAFFLE/INFBAF 

C0MM0N/FFINP/NUMTHP,THTPHI(2,MXFFP) 

COMMON/NFINP/NUMFPN , NFPXS (3 , MXNFP) 

REAL  NFPXS 
C 

C  OUTPUT  COMMONS 

C 

COMWON/TAPREC/RECRD(10) , IRFCRD(30) 

CCWMON/TAPRC1/ARECRD(10) 

CHARACTERM  ARECRO 

CGUMON/PRGVLS/NDIMPV,NUyARS,Ni,VSFP,NUMf  FP,NUMN'-'‘>,  w:vL(.; 

C0MMQN/SURAR3/AREAS (MXARS) 

COMMON/aDSVEC/IVECT(MXARS),B(MXARS).IPIVTR(MXFPS) 
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COMPLEX  TVECT 

COMMON/VELSPS/VEL (UXARS) , SP (MXARS) 

COMPLEX  VEL,SP 

COMMON/PDISL/POWER,OIRIND,SRCLVL 

COMMON/FFVALS/FFP (MXFFP) , PNRMFF (MXFFP) . IFFNRM , RMFNRM 
COMPLEX  FFP 

COMMON/ rSCQM/TGTSTH (MXFFP) 

COMMON/NFVALS/NFP(MXNFP),PNRMNF(MXNFP),INFNRM,RMNNRM 

COMPLEX  NFP 

COMMON/PTSCOM/PTSSP(MXARS) 

COMPLEX  PTSSP 

COMMON/EXTCOM/EXTPRS (MXFPS) , lEXTFG 
COMPLEX  EXTPRS 

COMMON/NBPRTC/ IRHSPT, NARSPT, NPTBLK , FRQPT 

COMMON/NBPRTS/SYMTPT 

CHARACTER*3  SYMTPT 

DIMENSION  CC (10),  TRNS(3),  I£LTS(8,300) 
real  xl(1000) ,yl(1000) 

CHARACTER«3  SYMTYP 
CHARACTERS  FLDTYP.TAPEID.PRTTYP 

INTEGER  XIZ AX . XNSEONS , XIRG . XNSU, XNSV . XIORDU, XI ORDV 

INTEGER  ICOOR.  lELEM 

COMPLEX  PMATX (342  ),VMATX(342  ) 

UDSIZE  =  342 
RUNIO  =  ’TES’ 

DATE=’ 14- APR-89’ 

CALL  INITCM 
CALL  OPNSFL 
RHQ  =  1000 
C  =  1500 

OPEN (UNIT=NUNPRT, FILEsRUNID// ’ . OUT ’ , STATUS= ’ NE» ’ , 

1  FORM='FORMATTED’) 


C  FREQUENCY  AND  SYMMETRY  INPUTS 

PI=ACOS(-1.0) 

FREQ=238.7 

SYMTYP=’REF' 

NBLKS=2  »♦  3 

IRHSYM  =  1 

CONVERT=1,0 

C  SURFACE  REGION  INPUTS 

ROTLIM=PI/NBLKS 
NSREG=  3 

DO  1  1=1,10 

1  CC(I)=0.0 

DO  2  I  =  1,3 

2  TRNS(I)=0.0 

C  TOP 

XIRG=  1 
XNSEQNS=  2 

CC{  1)=.6  *  CONVERT 

TRNS(1)=0.00E*00»  CONVERT 
TRNS(2)=0.00E«00«  CONVERT 
TRNS(3)=0.00E*00*  CONVERT 
XIZAX=»3 

XSUL=1,68  »  CONVERT 

XSUU=2.04  ♦  CONVERT 

XSVL=0. 000000 
XSVU=1.57 
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XNSU=2 

XNSV=3 

XI0RDU=4 

XI0RDV=16 

CALL  LDSURR (XIRG , XNSEQNS , CC , TRNS . X IZAX , 

1  XSUL, XSUU, XSVL . XSVU, XNSU, XNSV, XIORDU, XIORDV) 


C  OUTSIDE 
XIRG=  2 
XNSEQNS=  4 

CC(  1)=2.04  »  CONVERT 

TRNS(l)=0.00E-00«  CONVERT 
TRNS(2)=0.00E.00»  CONVERT 
TRNS(3)=0.00E*00»  CONVERT 
XIZAX=.3 

XSUL=0.0  »  CONVERT 

XSUU=.6  »  CONVERT 

XSVL=0.0 

XSVU=1.57 

XNSU=2 

XNSV=3 

XI0RCJ=4 

XI0RDV=13 

CA J.  LDSURR (XIRG, XNSEQNS , CC , TRNS , XIZAX, 

1  XSUL , XSUU , XSVL , XSVU , XNSU , XNSV , XIORDU , XIORDV) 


C  INSIDE 

XIRG=  3 
XNSEQNSi  4 

CC(  1)=1.68  »  CONVERT 

TRNS(l)r0.00E.00»  CONVERT 
TRNS(2)=0.00E-00«  CONVERT 
TRNS(3)=0.00E*00«  CONVERT 
XIZAX=-3 

XSUL=0  .  CONVERT 

XSUU=.8  •  CONVERT 

XSVL=0 

XSVU=1.57 

XNSU=2 

XNSV=:3 

XI0RDU=4 

XI0RDV=18 

CALL  LDSURR (XIRG , XNSE^NS , CC . TRNS , XIZAX , 

1  XSUL , XSUU, XSVL , XSVU, XNSU, XNSV, XIORDU, XIORDV) 


PLOTCHIEF  paraneters  are; 


RUNID  -  naee  for  graphics  files  created. 

NBLKS  1  -  displays  only  user-defined  surfaces. 

NBLKS  -  (defined  by  CHIff  driver)  displays  the  surfaces 
including  synaetries. 

SYUTYP  -  defined  in  CHIEF. 

ISUBDIV  1  -  displays  surfaces  using  CHIEF  defined 
paraaeters  NSU  and  NSV. 

0  -  displays  surfaces  using  optiaal  subdivisions 
for  viewing. 

AX  -  (REAL)  angle  in  degrees  to  rotate  view  about  X  aos 

AX  -  (REAL)  angle  in  degrees  to  rotate  vew  about  X  ax’.s 

A/  -  (REAL)  angle  in  degree.s  to  rotate  view  about  Z  a«-s 


CALL  R_0TCH;ET\’TESr2’ ,  NaXS,  Er«lXP,  0,  60  0,  10  ?,  0  0) 
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CALL  PRT0UT(’GE0«’,1) 

C  GENERATE  SURFACE  P  AND  V  MATRICES 

CALL  SURMAT (FREQ , SYMTYP, NBLKS . PMATX, VMATX, MDSIZE) 

CALL  PRTOUTCAREA',1) 

C  DECOMPOSE  MATRICIES 

CALL  DECOMM(SYMTYP, NBLKS, IRHSYM, PMATX, VMATX, MDSIZE) 

CALL  IOSUB(NDVELS,10,VEL,0) 

C  BLOCK  f  1 

va(  1)  =  (0.0000,0.0000) 

VEL(  2)  =  (0  0000,0  0000) 

VEL(  3)  =  (0  0000,0.0000) 

VEL(  ♦)  =  (0.0000,0.0000) 

VEL(  5)  =  (0  0000,0.0000) 

VEL(  6)  -  (0.0000,0.0000) 

Va(  7)  =  (0.0000,0.0000) 
va(  8)  =  (0  0000,0.0000) 

var  9)  =  (0  0000,0  0000) 

va(  10)  =  (0.0000,0.0000) 
va(  11)  =  (0.0000,0.0000) 
va(  12)  =  (0.0000,0.0000) 
va(  13)  =  (0.0000,0,0000) 
va(  14)  =  (0.0000,0,0000) 
va(  15)  =  (0.0000,0,0000) 
va(  16)  =  (0.0000,0.0000) 
va(  17)  =  (0.0000,0.0000) 
va(  18)  =  (0.0000,0.0000) 

CALL  I0SU8(NDVaS,l,Va,N»DVEC) 

C  GENERATE  SURFACE  PRESSURES 

C AU  SURPRS  (FREq ,  SYMTYP ,  NBLKS ,  I RHSYM ,  PMATX ,  VMATX ,  MDS IZE) 

CALL  PRTOUTCSP  ’,1) 

aOTYP  =  ’FAR’ 

NUMTHP  =  19 
DO  25  I  =  1, NUMTHP 
THTPHI(1,I)  =  (I-l)«5 
THTPHI(2,I)  =  0.0 
25  CONTINUE 

C  CALCULATE  FAR  FiaO  MATRICIES 

C AU  UDMAT  (FREq ,  SYMTYP ,  UDTYP ,  N^KS ,  PMATX ,  VMATX ,  MD5I ZE) 

IFFNRM  =  19 

C  CALCULAIE  FAR-FiaO  PRESSURES 

C  AU  FLDPRS  (FREq ,  aDTYP ,  NBLKS ,  I  RHSYM ,  PMATX ,  VMATX ,  MDS  IZE) 

CAU  PRTOUT(FLDTYP,0) 

STOP 

END 
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APPENDIX  B  -  TEST2.CRD  and  TBST2.C0N 


TEST2.CRD  is  the  coordinate  file. 


X 

Y 

Z 

X 

Y 

Z 

-l.SEe288 

-0.7643659 

-0,4413067 

-1.550288 

-0.7843659 

-0.4413087 

-1.536148 

-0.8717934 

-0.2502503 

-1.536148 

-0.8717934 

-0.2502503 

-1.493970 

-0.9730662 

-5.9966244E-02 

-1.493970 

-0.9730662 

-5.9966244E-02 

-1.424474 

-1.066453 

0.1262928 

-1.424474 

-1.066453 

0.1262928 

-1.328849 

-1.150358 

0.3053430 

-1.328849 

-1.150358 

0.3053430 

-1.208730 

-1.223347 

0.4741242 

-1.208730 

-1.223347 

0.4741242 

-1.066168 

-1.284172 

0.6297510 

-l;066160 

-1.284172 

0.6297510 

-0.9036023 

-1.331794 

0.7895637 

-0.9036023 

-1.331794 

0.7695637 

-0.7238098 

-1.365398 

0.8911722 

-0.7238098 

-1.365398 

0.8911722 

-0.5298645 

-1,384410 

0.9924980 

-0.5298645 

-1.384410 

0.9924980 

-0.3250613 

-1.388506 

1.071809 

-0.3250813 

-1.388506 

1.071809 

-0.1129605 

-1.377614 

1.127750 

-0.1129605 

-1.377614 

1.127750 

0.1028719 

-1.351922 

1.159364 

0.1026719 

-1.351922 

1.1S9364 

-1.904818 

-0.8185040 

-0,4725634 

1.758665 

-0.2590765 

-0.1495778 

-1.887648 

-0.9489513 

-0.2405663 

1,744526 

-0.3708224 

3.8985357E-02 

-1.836432 

-1.071926 

-9,5071S97E-03 

1.702347 

-0.4849768 

0.2218322 

-1.752046 

-1.136325 

0.2166646 

1.832852 

-0.5995883 

0.3958373 

-1.635929 

-1.287209 

0.4340828 

1.537227 

-0.7126979 

0.5580264 

-1.490070 

-1.375838 

0.6390312 

1.417108 

-0.8223723 

0.7056271 

-1.318959 

-1.449697 

0.8280068 

1.274546 

-0.9267365 

0.8361166 

-1.119557 

-1.507524 

0.9977794 

1.111980 

-1.024007 

0.9472845 

-0.9012382 

-1.548329 

0.145447 

0.9321876 

-1.112521 

1.037171 

-0.6657331 

-1.571415 

1.268485 

0.7382423 

-1.190766 

1.104299 

-0.4170677 

-1.676388 

1.364792 

0.5334591 

-1.257403 

1.147501 

-0.1594928 

-1.583163 

1.432720 

0.3213384 

-1.311295 

1.166039 

0.1025897 

-1.531965 

1.471108 

0.1055059 

-1.351520 

1,159596 

-2.009007 

-0.3067828 

-0.1771211 

2.113198 

-0.2049384 

-0.1183212 

-1.991837 

-0.4372306 

5.4a76041E-02 

2.096028 

-0.3406299 

0.1106484 

-1.940621 

-0.5602046 

0.2859351 

2.044810 

-0.4792459 

0.3326767 

-1.856234 

-0.6736033 

0.5121069 

1.960423 

-0.6184171 

0.5439687 

-1.740118 

-0,7754880 

0.7295250 

1.844307 

-0.7557644 

0.7409126 

-1.594258 

-0.8641173 

0.9344736 

1.698447 

-0.8689402 

0.9201420 

-1.421148 

-0.9379763 

1.123449 

1.525337 

-1.015668 

1.078594 

-1.223746 

-0.9958025 

1.293222 

1.327935 

-1.133783 

1.213559 

-1.005427 

-1.036608 

1.4408.:) 

1.109618 

-1.241264 

1.322731 

-0.7699220 

-1.059694 

1,563928 

0.8741109 

-1.336275 

1.404243 

-0.5212566 

-1.064667 

1.660234 

0.6254455 

-1.417192 

1.456703 

-0.2636815 

-1.051442 

1.728162 

0.3678704 

-1.482633 

1,479214 

-1.S992371E-03 

-1.020244 

1.766550 

0.1057881 

-1.531477 

1.471390 

-1.904818 

-0.8185040 

-0.4725634 

2.009007 

0,3067828 

0.1771211 

-1.887648 

-0.9489518 

-0.2405663 

1.991837 

0.1710913 

0.4060907 

-1.836432 

-1.071928 

-9.5071597E-03 

1.940621 

3.247B308E-02 

0.6281190 

-1.752045 

-1.185325 

0.2166648 

1.856234 

-0.1066958 

0.8394110 

-1.635929 

-1.287209 

0.4340628 

1.740118 

-0.2440431 

1.036355 

-1.490070 

-1.375838 

0.6390312 

1.594258 

-0.3772191 

1.215534 

-1.316959 

-1.449697 

0.8280068 

1.421148 

-0.5039472 

1.374036 

-1.119557 

-1.507524 

0,9977794 

1.223748 

-0.6220614 

1.509001 

-0.9012382 

-1.548329 

1.145447 

1.005427 

-0.7295428 

1.618173 

-0.6857331 

-1.571415 

1.268485 

0.7699220 

-0.8245541 

1.699686 

-0.4170677 

-1.576388 

1.364792 

0.5212566 

-0.9054713 

1.752146 

-0.1594926 

-1.563163 

1.432720 

0.2636815 

-0.9709113 

1.774656 

0.1025897 

-1.631965 

1.471108 

1.5992371E-03 

-1.019755 

1.766832 

-1.664477 

-0.2526447 

-0.1458644 

2.113198 

-0.2049384 

-0.1183212 

-1.640337 

-0.3600722 

4.6192011E-02 

2.096026 

-0.3406299 

0.1106484 

-1.598158 

-0.4613450 

0.2354780 

2.044810 

-0.4792459 

0.3326767 

-1.528663 

-0.5547321 

0.4217350 

1.960423 

-0.6184171 

0.5439687 

-1.433038 

-0.6386372 

0.6007853 

1.844307 

-0.7557644 

0,7409126 

-1.312919 

-0.7116260 

0.7695664 

1.698447 

-0.8889402 

0.9201420 

-1.170357 

-0.7724610 

0.9251933 

1.625337 

-1.015668 

1.078594 

-1.007791 

-0.8200727 

1.065003 

1.327935 

-1.133783 

1.213559 

-0.8279988 

-0.8536767 

1.186615 

1.109616 

-1.241264 

1.322731 

-0.6340534 

-0.8726891 

1.287940 

0.8741109 

-1.336275 

1.404243 

-0.4292701 

-0.8767845 

1.367252 

0.6254455 

-1.417192 

1.456703 

-0.2171495 

-0.B85B929 

1.423192 

0.3678704 

-1.4B2633 

1.479214 

-l,3170240E-03 

-0.8402009 

1.454806 

0.1057881 

-1.531477 

1.471390 
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1.6S4477 

0  2526447 

0.1458644 

-1.119557 

0  1103402 

-1.804443 

1.640337 

0.1408987 

0.3344277 

-0  9012381 

0  2178216 

-1.913615 

1.598158 

2  6744332E-02 

0.5172745 

-0.6657332 

0  3128329 

-1.995128 

1  528663 

-8.7867148E-02 

0.6912796 

-0  4170677 

0.3937501 

-2.047588 

1  433038 

-0.2009787 

0.8534633 

-0.1594928 

0  4591901 

-2.070098 

1  312919 

-0.3106510 

1.001069 

0  1025897 

0.5080343 

-2.062275 

1.170357 

-0  4160153 

1.131559 

-1  654477 

-0  2526447 

-0,1458644 

1.007791 

-0.5122858 

1.242707 

-1  640337 

-0.1408987 

-0.3344277 

0.8279988 

-0.6007999 

1  332613 

-1.598158 

-2,6744332E-02 

-0  5172745 

0.6340534 

-0.6790445 

1.399741 

-1.528663 

8.7867148E-02 

-0  6912796 

0.4292702 

-0  7456822 

1.442943 

-1.433038 

0  2009767 

-0.8534686 

0.2171495 

-0.7995740 

1.461481 

-1.312919 

0.3106510 

-1.001069 

1.3170240E-03 

-0.8397987 

1.455039 

-1.170557 

0.4150153 

-1.131559 

1.758665 

-0  2590765 

-0  1495778 

-1  007791 

0.5122858 

-1.242707 

1.744526 

-0.3708224 

3.8985357E-02 

-0  8279988 

0.6i4iJ/999 

- 1 . 33iui3 

1.702347 

-0  4849768 

0.2218322 

-0.6340534 

0.6790445 

-1.399741 

1.632852 

-0.5995883 

0.3958373 

-0  4292702 

0.7456822 

-1.442943 

1  537227 

-0.7125979 

0.SS80264 

-0  2171495 

0.7995740 

-1.461481 

1.417100 

-0.8223723 

0.7056271 

-1.3i73240E-03 

0  8397987 

-1.455039 

1.274546 

-0.9267365 

0  8361166 

-1.550288 

-0.7643659 

-0.4413067 

1.111980 

-1.024007 

0.9472645 

-1.536148 

-0.6526198 

-0  6298699 

0.9321876 

-1.112521 

1.037171 

-1.493970 

-0.5384656 

-0.8127167 

0.7382423 

-1.190766 

1.104299 

-1.424474 

-0.4238540 

-0  9867219 

0.5334591 

-1.257403 

1.147501 

-1.328849 

-0.3107444 

-1.148911 

0.3213384 

-1.311295 

1.166039 

-1.208730 

-0.2010702 

-1.296512 

0.1055059 

-1.351520 

1.159596 

-1.066168 

-9.6705869E-02 

-1.427001 

-1.550288 

-0.7643659 

-0.4413067 

-0.9036024 

5.6463428E-04 

-1.538149 

-1.536148 

-0.6526198 

-0  6298699 

-0.7238098 

8  9078687E-02 

-1.628055 

-1.493970 

-0.5384658 

-0.8127167 

-0.5298645 

0  1673233 

-1.695183 

-1.424474 

-0.4238540 

-0.9867219 

-0.3250813 

0.2339610 

-1.738386 

-1.328849 

-0.3107444 

-1.148911 

-0.1129606 

0  2678528 

-1.756924 

-1.208730 

-0.2010702 

-1  296512 

0.1028719 

0.3280776 

-1.750481 

-1.066168 

-9.6705869E-02 

-1.427001 

1.758665 

-0.2590765 

-0.1495778 

'0.9036024 

5.6463428E-04 

-1.538149 

1.744526 

-0.1516490 

-0.3406343 

'0.7238098 

8.9078687E-02 

-1.628055 

1.702347 

-5  0376166E-02 

-0.5309183 

'0.5298645 

0.1673233 

-1.695183 

1.632852 

4.3010946E-02 

-0.7171774 

-0.3250813 

0.2339610 

-1.738386 

1.537227 

0.1269160 

-0.8962277 

'0.1129608 

0.2878528 

-1.756924 

1.417108 

0.1999048 

-1.065009 

0.1028719 

0.3280776 

-1.750481 

1.274546 

0.2607298 

-1.220636 

'1.904818 

-0.8185040 

-0.4725634 

1.111980 

0.3083515 

-1.360448 

'1.887648 

-0.6828125 

-[B7015330 

0.9321876 

0.3419556 

-1.482057 

-1.836432 

-0.5441964 

-0.9235613 

0.7382423 

0.3609679 

-1.583383 

-1.752045 

-0.4050253 

-1.134853 

0.5334591 

0.3650633 

-1.662694 

-1.635929 

-0.2676781 

-1.331797 

0.3213384 

0.3541718 

-1.718634 

-1.490070 

-0.1345021 

-1.511027 

0.1055059 

0.3264796 

-1.750249 

-1.316959 

-7.7739796E-03 

-1.669478 

2.113196 

-0.2049384 

-0,1183212 

-1.119557 

0.1103402 

-1.804443 

2.096026 

-7.4490622E-02 

-0.3503184 

-0.9012381 

0.2178216 

-1.913615 

2.044810 

4.8483498E-02 

-0.5813774 

-0.6657332 

0.3128329 

-1.995128 

1.960423 

0.1618821 

-0.8075492 

-0.4170677 

0.3937501 

-2.047588 

1.844307 

0.2637668 

-1.024967 

-0.1594928 

0.4591901 

-2.070098 

1.698447 

0.3523961 

-1.229918 

0.1025897 

0.5080343 

-2.062275 

1.525337 

0.4262551 

-1.418891 

-2.009007 

-0.3067828 

-0.1771211 

1.327935 

0.4840814 

-1.588664 

-1.991837 

-0.1710913 

-0.4060907 

1.109616 

0.5248864 

-1.736331 

-1.940621 

-3.2475308E-02 

-0.6281190 

0.8741109 

0.5479728 

-1.859370 

-1  856234 

0.1066958 

-0.6394110 

0.6254455 

0.5529457 

-1.955676 

-1.740118 

0.2440431 

-1.036355 

0.3678704 

0.5397204 

-2.023604 

-1.594258 

0.3772191 

-1.215584 

0.1057881 

0.5085227 

-2.061993 

-1.421148 

0.5039472 

-1.374036 

2.009007 

0.3067828 

0.1771211 

-1.223746 

0.6220614 

-1.509001 

1.991837 

0.4372306 

-5.4876041E-02 

-1.005427 

0.7295428 

-1.618173 

1.940621 

0.5602046 

-0.2859351 

-0.7699220 

0.8245541 

-1.699686 

1.856234 

0.6736033 

-0.5121069 

-0.5212566 

0.9054713 

-1.752146 

1.740118  0 

7754880  -0 

.7295250 

-0  2636815 

0.9709113 

-1.774656 

1.594258 

0.8641173 

-0.9344736 

-1.5992371E-03 

1.019755 

-1.766832 

1  421148 

0  9379763 

-1.123449 

-1.904818 

-0.8185040 

-0.  •tr  25634 

i . 223746 

0.9958025 

-1.293222 

-1.887648 

-0.6828125 

-0.7015330 

1.005427 

1.036608 

-1.440889 

-1.836432 

-0.5441964 

-0.9235613 

0  7699220 

1.059694 

-1.563928 

-1.752045 

-0.4050253 

-1.134853 

0.5212566 

1.064667 

-1.660234 

-1.635929 

-0.2676781 

-1.331797 

0.2636815 

1  051442 

-1.728162 

-1  490070 

-0  1345021 

-1.511027 

1  5992371E-03 

1.020244 

-1.766550 

-1.316959 

-7.7739796E-03 

-1.669478 

2  113198 

-0.2049384 

-0.1183212 
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2  096226 

-7.4490622E-02 

-0.3503184 

-1.005427 

-1.036608 

1.440889 

2  044810 

4  8483490E-02 

-0,5813774 

-0.7699220 

-1,059694 

1.563928 

1.960423 

0  1618821 

-0,8075492 

-0  5212566 

-1.064667 

1  660234 

1,844307 

0.2637668 

-1.024967 

-0.2636815 

-1.051442 

1.728162 

1.698447 

0  3523961 

-1.229916 

-1.5992371E-03 

-1.020244 

1.766550 

1  525337 

0  4262551 

-1.418891 

-2  113196 

0.2049384 

0.1183212 

1.327935 

0,4840814 

-1.588664 

-2.096026 

7.4490622E-02 

0.3503184 

1.109616 

0  5248864 

-1.736331 

-2.044810 

-4.8483498E-02 

0.5813774 

0  8741109 

0.5479728 

-1.859370 

-1.960423 

-0.1618821 

0.8075492 

0  6254455 

0.5529457 

-1.955676 

-1.844307 

-0  2637668 

1  024967 

0.3678704 

0,5397204 

-2.023604 

-1.698447 

-0  3523961 

1.229916 

0.1057881 

0  5085227 

-2.061993 

-1.525337 

-0  4262551 

1.418891 

1  654477 

0.2526447 

0.1458644 

-1  327935 

-0.4840814 

1 . 588664 

1.540337 

0,3600722 

-4,5192011E-02 

-1.109616 

-0.5248864 

1.736331 

1.598158 

0.4613450 

-0.2354760 

-0  8741109 

-0  5479728 

1.859370 

1.528663 

0  5547321 

-0.4217350 

-0.6254455 

-0.5529457 

1.955676 

1.433038 

0  6386372 

-0  6007853 

-0  3678704 

-0  5397204 

2.023604 

1  3' 2919 

0.7116260 

-0  7695664 

-0  1057881 

-0  5085227 

2.061993 

1  173.8.57 

-'*.9251933 

-1.654477 

-0.2526447 

-0.1458644 

1,007791 

0  8200727 

-1.066006 

-1.640337 

-0.3600722 

4.5192011E-02 

0.8275980 

0  853676' 

-1.186615 

-1.598158 

-0  4613450 

0.2354760 

0.6340534 

0.8726?91 

-1.287940 

-1.528663 

-0  5547321 

0.4217350 

0.4292701 

0.876'/045 

-1.367252 

-1  433038 

-0.6386372 

0.6007853 

0,2171495 

0.8553929 

-1.423192 

-1.312919 

-0.7116260 

0.7695664 

1.3170240E-03 

0  8402009 

-1.454806 

-1.170357 

-0.7724510 

0.9251933 

1.758665 

-0.2590765 

-0.1495778 

-1.007791 

-0  8200727 

1.065006 

1  744526 

-0.1516490 

-0.3406343 

-0.8279988 

-0.8536767 

1.186615 

1.702347 

-5.0376166E-02 

-0.5309183 

-0.6340534 

-0.8726891 

1.287940 

1  632852 

4.3010946E-02 

-0.7171774 

-0.4292701 

-0.8767845 

1.367252 

1.537227 

0.1269160 

-0.8962277 

-0,2171495 

-0.8658929 

1.423192 

1.417108 

0.1999048 

-1.065009 

-1.3170240E-03 

-0.8402009 

1.454808 

1.274546 

0  2607298 

-1,220636 

-1.758665 

0.2590765 

0.1495778 

1.111980 

0.3083515 

-1.360448 

-1.744526 

0.1516490 

0.3406343 

0  321878 

0.3419556 

-1.482057 

-1.702347 

S.0376166E-02 

0.5309183 

0.7382423 

0.3609679 

-1.583383 

-1.632852 

-4. 3010941 E-02 

0.7171774 

0  5334591 

0.3650633 

-1.662694 

-1.537227 

-0.1269160 

0.8962277 

0  3213304 

0.3541718 

-1.718634 

-1.417108 

-0.1999048 

1.065009 

0.1055059 

0.3284798 

-1.750249 

-1.274546 

-0.2C 07298 

1.220636 

-1.750665 

0  259076S 

1495778 

-1.111980 

-0.308351S 

1.360448 

-1.744526 

0.1516490 

0.3406343 

-0.9321876 

-0.3419556 

1.482057 

-1.702347 

5.0376166E-02 

0.5309’' J 

-0.7382423 

-0,3609679 

1.583383 

-1  632852 

-4.3310946E-02 

0.71717,4 

-0.5334591 

-0.3650633 

1.66269*. 

-1.537227 

-0,1269160 

0.896227/ 

-0.3213384 

-0.3541718 

1.718634 

-1.417108 

-0,1999048 

1,065009 

-0.1055059 

-0.3284798 

1.750249 

-1.274546 

-0.2607298 

1.220636 

1.550288 

0.7643659 

0.4413067 

-1.111980 

-0.3083515 

1.360448 

1.536148 

0.6526198 

0.6298699 

-0.93218/6 

-0.3419556 

1.482057 

1.493970 

0.5384656 

0.8127167 

-0.7302423 

-0.3609679 

1 . 583383 

1.424474 

0.4238540 

0.9867219 

-0.5334591 

-0.3650633 

1.662694 

1.328849 

0,3107444 

1.148911 

-0.3213304 

-0,3541718 

1.718634 

1.208730 

0.2010702 

1.296512 

-0.1055059 

-0.3284798 

1.750249 

1.066168 

9.6705869E-02 

1.427001 

-2.113196 

0.2049384 

0.1183212 

0.9036024 

-5.6463428E-04 

1.538149 

-2.096026 

7,4490622E-02 

0.3503184 

0.7238098 

-8.9078687E-02 

1.628055 

-2.044810 

-4.8483498E-02 

0.5813774 

0.5298645 

-0.1673233 

1.695183 

-1.960423 

-0.1618821 

0.8075492 

0.3250813 

-0.2339610 

1.738386 

-1.844307 

-0.2637668 

1.024967 

0.1129606 

-0.2878528 

1.756924 

-1.698447 

-0.3523961 

1.229916 

-0.1028719 

-0.3280776 

1.750481 

-1.525337 

-0.4262551 

1.418891 

1.904818 

0.8165040 

0.4725634 

-1.327935 

-0.4840814 

1  588664 

1.887648 

0.6828125 

0.7015330 

-1.109616 

-0.5248864 

1.736331 

1.838432 

0.5441964 

0  9235613 

-0.8741109 

-0.5479728 

1  859370 

1.752045 

0.4050253 

1.134853 

-0.6254455 

-0.5529457 

1.955676 

1.635929 

0.2876781 

1.331797 

-0  3678704 

-0.5397204 

2.023604 

1  490070 

0.1345021 

1.511027 

-0  1057881 

-0  5085227 

2  061993 

1.316959 

7.7739796E-03 

1.669478 

-2.009007 

-0,3067828 

-0  1771211 

1  119557 

-0  1103402 

1  804443 

-1  .^(91837 

-0.4372306 

6  4876041E-02 

0,9012381 

-0.2178216 

1.913615 

-1  940621 

-0.560.;046 

0.2859351 

0.6657332 

-0  3128329 

1.995128 

-1  856234 

-0,6736033 

0,5121069 

0  4170677 

-0. 393/501 

2  0475S8 

-1.740110 

-0  7754880 

0  7295250 

0  1594926 

-0.4591901 

2.070098 

-1  594258 

-0  8641173 

0  9344736 

-0  1025897 

-0  5080343 

2.062275 

-1  421148 

-0  9379763 

1  123449 

2  009007 

0.3087628 

0.1771211 

-1  223746 

-0.9958025 

1  293222 

1  991837 

0  1710913 

0.4060907 
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1.940621 

3.2475308E-02 

0  6281190 

-0  8741109 

1.336275 

-1.404243 

1  856234 

-0.1066958 

0  8394110 

-0.6254455 

1.417192 

-1.456703 

1.740118 

-0.2440431 

1.036355 

-0  3678704 

1.482633 

-1.479214 

1  594258 

-0.3772191 

1.215584 

-0.1057881 

1.531477 

-1.471390 

1.421148 

-0.5039472 

1.374036 

-2.009007 

-0.3067828 

-0.1771211 

!  223746 

-0.6220614 

1.509001 

1.091837 

-0.1710913 

-0.4060907 

1.005427 

-0.7295428 

1.618173 

-1.940621 

-3.247S308E-02 

-0.6281190 

0  7699220 

-0.8245541 

1  699686 

-i. 856234 

0.1066958 

-0  8394110 

0.5212566 

-0.9054713 

1.752146 

-1.740118 

0.2440431 

-1.036355 

0  2636815 

-0  9709113 

1  774656 

-1.594258 

0.3772191 

-1.215584 

1.5992371E-03 

-1.019755 

1.766832 

-1.421148 

0.5039472 

-1.374036 

1.904818 

0.8185040 

0  4725634 

-1.223746 

0.6220614 

-1.509001 

1.887648 

0.6828125 

0.7015330 

-1.005427 

0.7295428 

-1.618173 

1.836432 

0.5441964 

0.9235613 

-0.7699220 

0  8245541 

-1.699686 

1.752045 

0.4050253 

1.134853 

-0.5212566 

0.9054713 

-1.752146 

1.635929 

0.2676781 

1.331797 

-0.2636815 

0.9709113 

-1.774656 

1.490070 

0.1345021 

1  511027 

-1.S992371E-03 

1.1)19755 

-1.766832 

1.316959 

7.7739796E-03 

1.669478 

-2.113196 

0.2049384 

0.1183212 

1.119557 

-0.1103402 

1.804443 

-2.096028 

0.3405299 

-0  1106484 

0.9012381 

-0.2178216 

1.913615 

-2.044810 

0.4792459 

-0.3326767 

0.6657332 

-0  3128329 

1.995128 

-1.960423 

0.6184171 

-0.5439687 

0.4170677 

-0.3937501 

2.047588 

-1.844307 

0.7557644 

-0.7409126 

0.1594928 

-0.4591901 

2.070098 

-1.698447 

0.8889402 

-0.9201420 

-0.1025897 

-0.5080343 

2.062275 

-1.525337 

1.015668 

-1.078594 

1.654477 

0  2526447 

0.1458644 

-1.327935 

1.133783 

-1.213559 

1.640337 

0.1408987 

0.3344277 

-1.109616 

1.241264 

-1.322731 

1.598158 

2.6744332E-02 

0  5172745 

-0.8741109 

1.336275 

-1.404243 

1.528663 

-8.7867148E-02 

0  6912796 

-0.6254455 

1.417192 

-1.456703 

1.433038 

-0.2009767 

0.8534686 

-0.3678704 

1.482633 

-1.479214 

1.312919 

-0.3106510 

1.001069 

-0.1057881 

1.531477 

-1.471390 

1.170357 

-0.4150153 

1.131559 

-1  654477 

-0  2526447 

-0.1458644 

1  007791 

-0.5122858 

1.242707 

-1.640337 

-0.1408987 

-0.3344277 

0.8279988 

-0.6007999 

1.332613 

-1.598158 

-2.6744332E-02 

-0.5172745 

0  6340534 

-0.6790445 

1.399741 

-1.528663 

8.7867148E-02 

-0.6912796 

0.4292702 

-0.7456822 

1.442943 

-1.433038 

0.2009767 

-0.8534686 

0.2171495 

-0.7995740 

1.461481 

-1.312919 

0.3106510 

-1.001069 

1.3170240E-03 

-0.8397987 

1.455039 

-1.170357 

0.4150153 

-1.131559 

1.550288 

0.7643659 

0.4413067 

-1.007791 

0.5122858 

-1.242707 

1.536148 

0.6526198 

0.6298699 

-0.8279988 

0.6007999 

-1.332613 

1.493970 

0.5384656 

0.8127167 

-0.6340534 

0.6790445 

-1.399741 

1.424474 

0.4238540 

0.9867219 

-0.4292702 

0.7456822 

-1.442943 

1.328849 

0.3107444 

1.148911 

-0.2171495 

0.7995740 

-1.461481 

1.208730 

0.2010702 

1.296512 

-1.3170240E-03 

0.8397987 

-1 . 45S039 

1.066168 

9.6705869E-02 

1.427001 

-1.758665 

0.2590765 

0.1495778 

0.9036024 

-5.6463428E-04 

1.538149 

-1.744526 

0.3708224 

-3.898S357E-02 

0.7238098 

-8.9078687E-02 

1.628055 

-1.702347 

0.4849768 

-0.2218322 

0.5298645 

-0.1673233 

1.695183 

-1.632852 

0.5995883 

-0.3958373 

0.3250813 

-0.2339610 

1.738386 

-1.537227 

0.7126979 

-0.5580264 

0.1129606 

-0.2878528 

1.756924 

-1.417108 

0.8223723 

-0.7056271 

-0.1028719 

-0.3280776 

1.750481 

-1.274546 

0.9267365 

-0.8361166 

-1.758665 

0.2590765 

0.1495778 

-1.111980 

1.024007 

-0.9472645 

-1.744528 

0.3708224 

-3.8905357E-02 

-0.9321876 

1.112521 

-1.037171 

-1.702347 

0.4849768 

-0.2218322 

-0.7382423 

1.190768 

-1.104299 

-1.632852 

0.5995883 

-0.3958373 

-0.5334591 

1.257403 

-1.147501 

-1.537227 

0.7126979 

-0.5580264 

-0.3213384 

1.311295 

-1.166039 

-1.417108 

0.8223723 

-0.7056271 

-0.1055059 

1.351520 

-1.159598 

-1.274546 

0.9267365 

-0.8361166 

1.550288 

0.7643659 

0.4413067 

-1.111980 

1.024007 

-0.9472645 

1.536148 

0.8717934 

0.2502503 

-0.9321876 

1.112521 

-1.037171 

1.493970 

0.9730662 

5.9966244E-02 

-0.7382423 

1.190768 

-1.104299 

1.424474 

1.066453 
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611 

-255 

521 

534 

-2 

612 

613 

-255 

525 

538 

-2 

613 

614 

-255 

529 

542 

-2 

614 

615 

-255 

533 

546 

-2 

615 

616 

-255 

547 

548 

255 

616 

617 

-255 

548 

549 

255 

617 

618 

-255 

549 

550 

255 

618 

619 

-255 

550 

551 

255 

619 

620 

-255 

551 

552 

255 

620 

621 

-255 

552 

553 

255 

621 

622 

-255 

553 

554 

255 

622 

623 

-255 

554 

555 

255 

623 

624 

-255 

555 

556 

255 

599 

612 

-255 

556 

557 

255 

603 

616 

-255 

557 

558 

265 

607 

620 

-255 

558 

559 

255 

611 

624 

-255 

560 

561 

255 

625 

626 

5 

561 

562 

255 

625 

627 

5 

562 

563 

255 

625 

628 

5 

563 

564 

255 

629 

628 

5 

564 

565 

255 

626 

630 

5 

565 

566 

255 

630 

629 

5 

566 

567 

255 

631 

627 

5 

567 

568 

255 

627 

632 

5 

568 

569 

255 

632 

631 

5 

569 

570 

255 

633 

628 

5 

570 

571 

255 

628 

634 

5 

571 

572 

255 

634 

633 

5 

547 

560 

255 

635 

636 

5 

551 

564 

255 

637 

638 

5 

555 

568 

255 

639 

640 

5 

559 

572 

255 

640 

641 

5 

573 

574 

255 

640 

642 

5 

574 

575 

255 

643 

644 

5 

575 

578 

255 

644 

645 

5 

576 

577 

255 

645 

646 

5 

577 

578 

255 

578 

579 

255 

579 

580 

255 

580 

581 

255 

581 

582 

255 

582 

583 

255 

583 

584 

255 

584 

585 

255 

586 

587 

255 

587 

588 

255 

588 

589 

255 

589 

590 

265 

590 

591 

255 

591 

592 

255 

592 

593 

255 

593 

594 

255 

594 

595 

255 

595 

598 

255 

596 

697 

255 
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APPENDIX  C  -  Two  Examples  of  Graphics  Drivers 


DRAW240  -  written  by  D.  Mitchell  Bolling  Jr. 

TRI/TESSCQ,  Inc. 

August  1989 

Purpose:  This  prograw  is  an  elementary  graphics  driver.  DRAW240 
loads  in  the  coordinate  and  connectivity  files  that  result  from 
PLOTCID  or  PLOTCHIEF  and  displays  the  figure  on  the  screen  in  REGIS 
format.  It  way  be  used  on  any  VT240  or  compatible  terminal  or  any 
term i na I /computer  using  a  VT240  emulator  program. 

Subroutines  used:  None. 

Modification  Log: 


PsRAMETFR  PLT  UNIT  =  6 
PARAMETER  CRD' UNIT  =  0 
PARAMETER  CQN'uNIT  =  9 
PARAMETER  MAXCOOROS  =  10000 


lun  number  of  output  (6  =  screen)  default 
lun  number  of  coordinate  file 
lun  number  of  connectivity  file 
maximum  number  of  coordinates  for  plotting 


C . 

C  The  following  variables  are  used: 

C . 

CHARACTER*50  TEMP  I  a  temporary  string  variable 
CHARACTER«20  RUNIO  !  driver  name  used  in  CID  or  CHIEF  run 


REAL 

COORDS (3, 

MAXCOOROS) 

coordinates  read  in  froa  coordinate  file 

REAL 

NCOOROS(3, 

,  MAXCOOROS) 

coordinates  modified  by  scaling 

REAL 

MAXX 

largest  x  coordinate 

REAL 

MAX'.' 

largest  y  coordinate 

REAL 

MINX 

smallest  x  coordinate 

REAL 

MINY 

seal  lest  y  coordinate 

REAL 

OISX 

largest  distance  between  2  x  values 

REAL 

OISY 

largest  distance  between  2  y  values 

REAL 

OIS 

maximum  distance 

REAL 

OFFSETX 

offset  for  centering  horizontally 

REAL 

OFFSETY 

offset  for  centering  vertically 

INTEGER  C0NNECTS(3,  MAXCOORDS)  !  connectivity  data  read  in  from  file 
INTEGER  MAXP  ' 

INTEGER  MAXC 
INTEGER  I 
INTEGER  H 
INTEGER  V 
INTEGER  COLOR 


numbs,,  of  points 
!  number  of  lines 
!  scratch  counter 

!  number  of  pixels  horizontally  (0-799) 

!  number  of  pixels  vertically  (0-479) 

I  number  associated  with  different  colors 


C - - - - - 

C  Prompt  the  user  for  filename. 

C . — . . . - - - 

WRITE(6,10)  'Enter  filename  used  in  CIO  or  CHIff  run: 
10  FORkiAT('l',A) 

READ(5,20)  LEN,  TEMP 
20  FORMAT (0, A) 

TYPE  », '  ' 

C . . - - - 

C  strip  off  the  file  extension. 

C - - - - - 

L  =  IN0EX(Ta(P,  ’  .  ’) 

IF  (L  GT.  0)  THEN 
LEN  =  L  -  1 
RUNIO  =  TEMP(1:  LEN) 

asE 

RUNIO  =  TEMP 
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END  IF 


Vw - 

C  Read 

in  the  coordinate  data. 

TEMP  =  RUNID(1:LEN)  //  ’.CRD* 

OPEN  (UNIT=CRD  UNIT,  FILE=TEMP,TYPE=’OLD’) 
MAXP  =  1 

100 

READ  (CRD  UNn,»,ERR-125)  CD0RDS(1,  MAXP), 

COORDS (2,  MAXP), 

♦ 

COORDS (2,  MAXP)  =  -COORDS (2,  MAXP) 

MAXP  =  MAXP  ♦  1 

COORDS (3,  MAXP) 

IF  (MAXP  .LE.  MAXCOORDS)  GOTO  100 

125 

MAXP  =  MAXP  -  1 

r _ 

CLOSE  (CRD  UNIT) 

C  Read  in  the  connectivity  and  color  data. 

C - - - - - - - 

TEMP  =  RUNI0(1:LEN)  //  '.CON* 

OPEN  (UNIT=CON_UNIT.  FILE=TEMP,TYPE=’OLD’) 

MAXC  =  1 

200  READ  (C0N_UNIT,v,ERR=225)  CQNNECTS(1.  MAXC), 

*  C0NNECTS(2,  MAXC),  CQNNECTS(3,  MAXC) 

MAXC  =  MAXC  *  1 

IF  (MAXC  .LE.  MAXCOORDS)  GOTO  200 

225  MAXC  =  MAXC  -  1 

CLOSE  (C0N_UNIT) 

C . . . . . . 

C  Copy  the  coordinates  into  a  aorl(in9  array  A  possible  enhancenent 

C  aight  be  to  pass  this  working  array  to  R0T3D(),  allowing  the  user 

C  to  interectively  rotate  the  view. 

C . - . 

DO  I  =  1.  MAXP 
DO  J  =  1,  3 

NCOORDS(J,I)  =  COORDS(J.I) 

END  DO 
END  DO 

300  MINX  =  99999999999.0 

MINY  =  MINX 
MAXX  =  -MINX 
MAXY  =  -MINX 

C . . - . . . . . . 

C  Find  the  eaxieua  and  einieue  values  of  the  figure  in  the  x  and  y 

C  directions. 

C . — - - - - - 

DO  I  =  1,  MAXP 

IF  (NC00RDS(1,I)  .GT.  MAXX)  MAXX  =  NC00RDS(1,I) 

IF  (NC00RDS(1,I)  .LT.  MINX)  MINX  =  NC00RDS(1,I) 

IF  (NCOORDS(2,I)  .GT.  MAXY)  MAXY  =  NC00RDS(2,I) 

IF  (NC00RDS(2,I)  .LT.  MINY)  MINY  =  NC00RDS(2,I) 

END  DO 


C . - - - - - . . . 

C  Find  the  eaxieua  width  and  height  of  the  figure. 

C - - - - 

DISX  =  MAXX  -  MINX 
DISY  =  MAXY  -  MINY 

C . . . - . - . . 

C  Find  out  which  is  greater,  the  eaxieua  width  or  eaxieua  height. 

C- . . . 

IF  (DISX  GT.  DISY)  THEN 
DIS  =  DISX 
ELSE 

DIS  =  DISY 
END  IF 

OFFSETX  =  400  -  (521  »  DISX)  /  (2  *  DIS) 

OFFSETY  =  240  -  (431  «  DISY)  /  (2  «  DIS) 
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C - - - - - - 

C  Scale  the  coordinates  so  that  the  drawing  will  fit  on  the  screen 
C - - - 

DO  I  =  1,  MAXP 

NC00RDS(1,I)  =  521  »  (NCQORDS(l.I)  -  MINX)  /  DIS  ♦  DFFSETX 
NC00RDS(2,I)  =  431  »  (NC00RDS(2, I)  -  MINY)  /  OIS  ♦  OFFSETY 
END  DO 


C - - - 

C  Set  Regis  graphics  eode  and  clear  the  screen. 

C - - - - - - - 

WRITE(PLT  UNIT,*)  CHAR(27)  //  'P0p’ 
WRITE(PLT  UNIT,*)  'S(E)' 

COLOR  =  0 


C - 

C  Create  the  picture. 

C - 

DO  I  =  1,  MAXC 

C - 

C  If  .  lere  is  a  color  change,  then  change  the  color. 
C - - — . — . 

IF  (CONNECTS (3, I)  NE.  COLOR)  THEN 
COl  OR  =  CQNNECTS(3,I) 

WRITE(PLT_ UNIT, 900)  ABS (COLOR) 

900  F0RUAT('JW(r,I3,’)’) 

IF  (COLOR  LT.  0)  THEN 

WRITE(PLT_UNIT,*)  ’W(P2)’ 

ELSE 

»RITE(PLT  UNIT,*)  ’»(P1)’ 

END  IF 
END  IF 

C . - . 

C  Move  the  cursor  to  first  coordinates. 

C . 

H  =  NCOORDS(1,CONNECTS(1,I)) 

V  *  NC00RDS(2,C0NNECTS(1,I)) 

IRITECPlT  unit, 1000)  H,  V 

1000  F0RMAT('J/[’,I4,',’,I4,’]’) 

C . — . 

C  Draw  a  line  to  the  seconc  coordinates. 

C . 

H  =  NC00RDS(1,C0NNECTS(2,I)) 

V  =  NC00RDS(2,C0NNECTS(2,I)) 
WRITE(PLT_UNIT,U00)  H,  V 

1100  FORMATC  V[’,I4,‘,’,I4,’r) 

END  DO 


C . - . 

C  Exit  ReGIS  aode  and  end. 
C . . . . 


WRITE(PLT_UNIT,*)  CHAR(27)  //  ’\’ 
END 
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DRAW3O.PR0  -  written  by  0.  Mitchell  Bolling  Jr. 

TRI/TESSCO,  Inc. 

August  1989 

Purpose;  This  program  is  an  elementary  graphics  driver.  DRAfSD.PRO 
is  a  PVWAVE  program  that  loads  in  the  coordinate  and  connectivity 
files  that  result  from  the  PLOTCID  or  PLOTCHIEF  programs  and  displays 
the  figure  on  the  screen  using  the  graphics  protocol  defined  in  the 
PVWAVE  software.  PVWAVE  is  a  data  visualization  and  color  graphics 
package  from  Precision  Visuals. 

Subroutines  used:  None. 


FILENAME  =  ’ ’ 

READ,  'Enter  filename  of  plot  data.  ’.FILENAME 


Read  in  the  coordinate  data. 


GET_LUN,  F 

OPENR,  F,  FILENAME^’.CRD’ 

PI  =  D0U8LE(3. 14159265358979) 

I  =  ^ 

RECORD  =  0BLARR(3) 

WHILE  NOT  EOF(F)  DO  BEGIN 
REAOF,  F,  RECORD 
IF  (I  ED  0)  THEN  BEGIN 
COORDS  =  RECORD 
I  =  1 

ENDIF  ELSE  BEGIN 

COORDS  =  [[COORDS], [RECORD]] 

endb.se 

ENDWHILE 
CLOSE,  F 


Read  in  the  connectivity  data. 


OPENR,  F,  FILENAME- ’.CON’ 

I  =  0 

RECORD  =  INTARR(3) 

WHILE  NOT  EOF(F)  DO  BEGIN 
READF,  F,  RECORD 
IF  (I  EH  0)  THEN  BEGIN 
CONNECTS  =  RECORD 
I  =  1 

ENDIF  ELSE  BEGIN 

CONNECTS  =  [[CONNECTS], [RECORD]] 
ENOaSE 
ENDWHILE 
CLOSE,  F 
FREE_LUN,  F 


;  Adjust  connectivity  data  due  to  PVWAVE’s  indexing  starting  at  0. 
CONNECTS(0,0)  =  CONNECTS (0:1,«)  -  1 


If  the  color  is  negative,  then  lines  should  be  dashed  to  represent 
inward  normal.  Outward  normal  is  shown  by  solid  lines. 


LINaVPES  =  CONNECTS (2,*)  LT  0 
CONNECTS (2,0)  =  A8S (CONNECTS (2,*)) 
NCOORDS  =  COORDS 


Determine  the  scaling  for  the  plot. 
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ITYPE  =  4  ♦  8  ♦  32  ‘  64  ;  inhibit  x  and  y  axis  rounding 

MINX  =  UIN(NCOORDS(0,»)) 

MAXX  =  MAX(NCOORDS(0,»)) 

MIDX  =  (MINX  -  MAXX)  /  2 
DISX  =  MAXX  -  MINX 

MINY  =  MIN(NC00RDS(1,«)) 

MAXY  =  MAX (NCOORDS (!,»)) 

MIDY  =  (MINY  ♦  MAXY)  /  2 
DISY  =  MAXY  -  MINY 

DIS  =  .55  »  MAX([DISX,DISY]) 


Use  the  whole  screen  for  plotting. 


!SC1  =  0 
!SC2  =  800 
ISC3  =  0 
ISC4  =  480 


,  Left  side  of  screen 
,  Right  side  of  screen 
;  Bottom  of  screen 
;  Top  of  screen 


Create  the  plot  array.  XCOORDS  will  contain  the  starting  and  ending 
X  coordinates  for  all  the  line  segments  to  be  drawn.  YCOORDS  will 
contain  the  starting  and  ending  Y  coordinates  for  ail  line  segments 
to  be  drawn. 


XCOORDS  =  NCOOROS(0,*) 

YCOORDS  =  NCOORDS  (1,») 

XCOORDS  =  XCOORDS (CONNECTS (0:1,*)) 
YCOORDS  =  YCOORDS (CONNECTS (0.1,*)) 


Draw  the  object 


SET  XY,  (MI0X-1.56*0IS),  (MIDX*1 . S6«0IS) ,  MIOY-DIS,  MIDY*DIS 
1C  =  0 

LINEPLOTXY,  XCOORDS,  YCOORDS,  C0NNECTS(2,») ,  LINETYPES 
END 
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(Blank  page) 
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APPENDIX  D  -  Listing  of  Programs  and  Subroutines  in  DISPLAY3D 


PLOTCID  -  written  by  D  Mitchell  Bolling  Jr. 

TRI/TESSCQ,  Inc. 

August  1989 


Purpose;  The  PLOTCID  program  is  a  graphics  routine  used  to  display 
three  dimensional  surfaces  defined  by  the  CHIEF  program. 

CHIEF  (Combined  Helmholtz  Integral  Equation  Formulation),  developed 
at  N05C,  computes  the  acoustic  radiation  from  arbitrary-shaped  bodies. 


PLOTCID  creates  a  plot  file  which  cai 
compiling,  linking  and  executing  the 
the  surface  information  from  the  CID 

Subroutines  used:  LOAD  CID  DATA() 
CIRCULAR  PLANE 0 

circular'cylinoer 

SPHEROIDO 

TOROIDO 

AXISYUMETRICO 

CONEO 

CONNECT^ ARROWS  0 


I  be  viewed  and  rotated  before 
CHIEF  driver.  PLOTCID  reads  in 
(Ciiief  Interactive  Driver)  program. 

RECTANGULAR  PLANEO 
ELLIPTICAL  FLANEO 
)  ELLIPTICAL_CYLINDER() 

PROL ATE^  OBLATE_  SPHERO I D ( ) 
OUAORILATERALC) 

TRIANGLEO 

R0T30() 

ORAW_ ARROWS 0 


Modification  Log. 

AUG  89  -  Added  OUADRILATERAL() ,  AXISYHMETRIC()  and  TRIANGLE() . 

Also  modified  program  by  adding  ROT30  so  that  the  user  can  rotate 
the  points  from  a i thin  PLOTCHIEF.  This  requires  that  the  points 
are  stored  in  an  array,  then  passed  to  R0T3D,  which  writes  them  out 
to  a  file.  Before,  the  points  were  sent  straight  to  a  file. 

JAN  90  -  Added  CONNECT_ARROWSp  and  ORAW_ ARROWS {)  in  order  to 
better  show  the  axis  orientation  during  rotations. 

JU.  90  -  Added  CONE()  geometry  to  PLOTCID. 


PARAMETER  PI  =  3.14159.  o36 

PARAMETER  COLOR  USERDEFINED  =  255  !  color  to  use  for  user  defined  part 
PARAMETER  COLOR_TRANSFORMED  =  2  !  color  for  transformed  part 

PARAMETER  CROUNIT  =  8  !  lun  number  of  coordinate  file 

PARAMETER  CONUNIT  =  9  !  lun  number  of  connectivity  file 

PARAMETER  MXSREG  =  500  !  maximum  number  of  surface  regions 

PARAMETER  MAXCOORDS  =  10000  !  ma'imum  number  of  coordinates  for  plotting 


ist  of  global  variables: 

INTEGER 

NSREG 

number  of  CHIEF  surfaces 

INTEGER 

NSERNS (MXSREG) 

types  of  CHIEF  surfaces 

REAL 

SUL (MXSREG) 

lower  1 imits  of  U 

REAL 

SUU (MXSREG) 

upper  limits  of  U 

REAL 

SVL (MXSREG) 

lower  1 imits  of  V 

REAL 

SVU (MXSREG) 

upper  limits  of  V 

REAL 

CCS (10, MXSREG) 

constants  needed  for  CHIEF  equat 

ions 

REAL 

TRNSS (3, MXSREG) 

translation  from  local  to  global 

origin 

INTEGER 

IZAX (MXSREG) 

global  axis  that  corresponds  to 

loca  1 

!  Z  axis 


COMMON/SVALS/NSREG, NSEONS (MXSREG) , SUL (MXSREG) ,SUU (MXSREG) , 

S VL (MXSREG) ,SVU (MXSREG) ,NSU (MXSREG) .NSV (MXSREG) 
CCS ( 10 , MXSREG) , TRNSS (3 , MXSREG) , IZAX (MXSREG) , 
lORDU(MXSREG) , lORDV (MXSREG) ,NCCEHS 

INTEGER  SYMTYP  !  type  of  symmetry:  reflective  (1), 

!  rotational  (2),  none  (0) 

INTEGER  N0LKS  !  number  of  symmetry  blocks 

!  for  reflective  symmetry: 
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-s'- 

SYMX 

!  if  -1,  then  reflects  about 

YZ 

plane 

SYMY 

I  if  -1,  then  reflects  about 

XZ 

p  lane 

SYMZ 

!  if  -1,  then  reflects  about 

XY 

p  1  ane 

!  for  rotational  symmetry: 

SYMANG 

1  radians  to  rotate  object 

plot  INFO  /  SYUTYP,  NBLKS,  SYMX,  SYMY,  SYMZ,  SYMANG 


■’T‘.  CG0RD3D  !  array  of  3d  points 

'•-.Au  ROTATED  !  3d  points  after  rotations 

common  /  ROTINFO  /  C00RD3D(3,  MAXCOORDS),  R0TATED(3,MAXC00RDS) 


of  iocal  variables: 


C‘-*,RACTER*S0  TEMP  !  a  teaporary  string  variable 
CmaRACTER»20  RUNID  !  driver  naae  used  in  CID  run 


'S^ECER  LEN,  L 
INTEGER  ISUBOIV 
INTEGER  COLOR 
INTEGER  POINTCNT 
INTEGER  SHOW_NORMALS 
INTEGER  REGION 
INTEGER  BLOCK 


!  holds  length  of  string  RUNID 
I  set  when  NSV  and  NSU  are  to  be  used 
!  nuaber  associated  with  different  colors 
I  keeps  track  of  next  point  to  be  created 
!  when  set,  uses  dotted  lines  for  inward  noraals 
!  counter  used  for  looping  through  each  region 
I  counter  used  for  looping  through  each  syaaetry 
I  block  (whether  reflective  or  rotational) 


’LAL  AX  !  number  of  degrees  to  rotate  view  about  x  axis 
REAL  AY  !  number  of  degrees  to  rotate  view  about  y  axis 
REAL  AZ  !  number  of  degrees  to  rotate  view  about  z  axis 


T.  me  or  all  of  these  constants  are  eultiplied  by  each  data  point 


REAL  REFLECTX(8)  !  array  of  8  constants  for  reflecting  about  yz  plane 
REAL  REFLECTY(8)  !  array  of  8  constants  for  reflecting  about  xz  plane 
REAL  REFLECTZ(8)  !  array  of  8  constants  for  reflecting  about  xy  plane 


C . - . 

C  Below  is  a  table  that  indicates  the  number  of  constants  needed  for 
C  each  reflective  symmetry  option.  Thus,  the  dashed  horizontal  line 
C  is  over  the  constants  needed  for  each  option. 

C . - . - . 

C  3  plane  symmetry  |< - >| 

C  2  plane  symmetry  |< - >|  j 

C  1  plane  symmetry  j< - )|  |  | 

C  III  I 

DATA  RERECTX  /  1,  -1,  1,  -1,  1,  -1,  1,  -1  / 

DATA  REFLECTY  /  1,  1,  -1,  -1,  1,  1,  -1,  -1  / 

DATA  RERECTZ  /  1,  1,  1,  1,  -1,  -1,  -1,  -1  / 


C - - - - 

C  Get  CID  run  name  from  user. 

r _ _ _ 

WRITE(6,100)  'Enter  the  filename  used  in  the  CIO  run: 
100  F0RMAT(’S’,A) 

REA0(5,125)  LEN,  TEMP 
125  FORMAT (Q, A) 


L  =  INDEX(TEMP, ' . ’) 

IF  (L  .GT.  0)  THEN 
LEN  =  L  -  1 
RUNID  =  TEMP(1  :  LEN) 

ELSE 

RUNID  =  TEMP 
END  IF 

C - - - 

C  Read  global  data  from  the  CID  output  file 
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C-  - - - 

TEWP  =  RUNID(1:LEN)  //  ’.DAT’ 

CALL  LOAD_CIO_DATA(TEMP) 

C - - - - - - - 

C  If  the  «odel  is  using  reflective  syneetry,  then  N8LKS  needs  to  be 
C  converted  from  the  number  of  planes  of  symmetry  to  the  total  number 
C  of  blocks. 

C - - - 

IF  (SYMTYP  .EQ.  1)  NBLKS  =  2  «•  NBLKS 

C - - - - - 

C  Ask  user  if  rotations  or  reflections  should  be  drawn. 

C- - - - - 

IF  (SYMTYP  .NE.  0  AND.  NBLKS  .GT.  1)  THEN 
IRITE(8,200)  ■  ' 

WRTTE{6,200)  ■  ’ 

WR1TE(6.200)  '  DISPLAY:' 

WRITE(6,200)  '  1.)  Only  the  user  defined  surfaces. ’ 
WRITE(6,200;  '  2.)  User  defined  surfaces  with  selected’  // 

♦  ’  transformations. ’ 

WRITERS  200)  '  ' 

*RriE(6,100)  'Enter  1  or  2  (default  2):  ’ 

PEAD(5,'50)  L 

150  format(bn:i0) 

IF  (L  Eq.  1)  NBLKS  =  1 
END  IF 


C  Find  out  if  NSU  and  NSV  subdivisions  should  be  used  in  plot. 


ISUBDIV  =  0 
WRTTE(6,200)  '  ’ 

WRITE(6,200)  ’  ’ 

»RIT£(8,200)  ’  THE  NUMBER  OF  SUBDIVISIONS  PLOHED:’ 
WRITE(6,200)  '  1.)  The  optimal  number  for  visualization.’ 
WRITE(8,200)  '  2.)  The  number  defined  in  the  CHIEF  driver.’ 
WRITE(6,200)  ’  ’ 

WRITE(6,100)  'Enter  1  or  2  (default  1);  ’ 

READ(5,150)  L 

IF  (L  .Eq.  2)  ISUBDIV  -  1 


C . - . 

C  Find  out  if  the  user  wishes  to  use  dotted  lines  to  represent 
C  negative  normals. 

C . . . . . . 

WRITE(6,200)  ’  ’ 

WRITE(8,200)  ’  ’ 

fRITE(6,200)  ’  DISTINGUISH  BETHEEN  POSITIVE  AND  NEGATIVE’  // 

♦  ’  NORMAL  VaOCITY:  ’ 

WRITE(6,200)  ’  1.)  Yes.  (solid  lines  for  positive,’  // 

*  ’  dashed  for  negative)’ 

IRITE{6,200)  ’  2.)  No.  (only  solid  lines)’ 
fRITE(8,200)  ’  ’ 

WRITE(6,100)  'Enter  1  or  2  (default  2):  ’ 

REAO(5,150)  L 

IF  (L  .Eq.  1)  THEN 
SHOf_NORMALS  =  1 
ELSE 

SHOW_NORMALS  =  0 
END  IF 


C . . . . . . 

C  Get  angles  of  rotation  (in  degrees)  about  the  3  axes  from  the  user. 

C . - - - - - - - 


200 


WRITE(6,200) 

WRITE(6,200) 

WRITE(6,200) 

* 

FORMAT (A) 

fRITE(8,100) 

fRITE(6,200) 


I  I 
I  I 

’  Enter  the  angles  (In  degrees)  for  rotating  '  // 

'the  model ’ 

'about  the  X,  Y,  and  Z  axes,  (default  00,  0  0,  00)’ 
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WRITE(6,100)  '(XROT,  YROT,  ZROT) 
R£AD(S,22S)  AX,  AY,  AZ 
225  FORMAT(3(E20  0)) 

C- - - - - 

C  Initialize  nan  ab  les 

C - 

POINTCNT  =  0 
SYMX  =  1 
SYMY  =  1 
SYUZ  =  1 


C - 

C  Open  the  connectivity  file  PLQTCID.CON  «i I  I  contain  the  line 
C  drawing  data  in  the  foreat:  STARTPOINT  ENDPOINT  LINECDLOR 

C  Note:  The  actual  coordinate  data  will  be  written  to  PLOTCID.CRD 
C  in  the  RQT3D()  subroutine  located  near  the  end  of  this  progra*, 

C  Since  the  coordinates  Bust  be  rotated  before  they  can  be  written  to 
C  a  f  i  le 

V,--'---------- - - - - - - - - - - - - - - - 

TEMP  =  RUNID(1:LEN)  //  '  CON’ 

OPEN  (UNI T=C0N  UN I T , F ILE=TEMP , TYPE= ’ NEW ‘ ) 


C . - - - 

C  Generate  points  in  each  syeietry  block. 

C . - . 

00  BLOCK  =  NBLKS,  1,  -1 


IF  (SYUTYP  EO.  1)  THEN 


C  For  current  block,  deterwines  which  planes  of  reflection 

C . 

SYMX  =  REFLECTX (BLOCK) 

SYMY  =  RERECTY  (BLOCK) 

SYMZ  =  REFLECTZ(BLOCK) 

ELSE  IF  (SYMTYP  .ER.  2)  THEN 


C  For  current  block,  deteriines  angle  to  rotate 

C . 

SYMANG  =  (2  a  PI  •  (BLOCK  -  1)  /  NBLKS) 

END  IF 


C . 

C  Loop  through  all  the  types  of  surfaces  defined  and  call  the  proper 
C  shape  drawing  subroutine. 

C . 

00  REGION  ^  1,  NSREG 


C .  - - - - - . 

C  User  defined  block  will  be  a  different  color  froa  the  coaputer 
C  generated  block(s).  If  IZAX  is  negative  (inward  noraal)  then  aake 
C  color  negative.  The  graphics  driver  should  draw  a  dotted  line  when 
C  colors  are  negative. 

C  -  . - . . . . . - . 

IF  (BLOCK  ER.  1)  THEN 

IF  (SH0W_N0RUALS  .ER.  1)  THEN 

COLOR  =  ISIGN(C0L0R_ USER, DEFINED, IZAX(REGION)) 

ELSE 

COLOR  =  C0L0R_USER_ DEFINED 
END  IF 
asE 

IF  (SH0f_N0RMALS  ER  1)  THEN 

COLOR  =  ISIGN(C0L0R_TRANSF0RMED,IZAX(REGI0N)) 
aSE 

COLOR  =  COLOR_TRANSFORMED 
END  IF 
END  IF 

IF  (NSERNS (REGION)  .ER.  1)  THEN 

CALL  RECTANGULAR  aANE(REGI0N,  POINTCNT, 

♦  COLOR,  CON  UNIT,  ISUBOIV) 
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ELSE  IF  (NSEQNS (REGION)  .EQ.  2)  THEN 
CALL  CIRCULAR  PLANE(REGIQN,  POINTCNT, 

COLOR,  CON  UNIT,  ISUBDIV) 

ELSE  IF  (NSEQNS (REGION)  .EQ.  3)  THEN 

CALL  ELLIPTICAL, PLANE(REGION,  POINTCNT, 

COLOR,  CON_UNIT,  ISUBDIV) 


ELSE  IF  (NSERNS (REGION)  .ER.  A)  THEN 

CALL  CIRCULAR  CYLINDER (REGION,  POINTCNT, 

COLOR,  CON_UNIT,  ISUBDIV) 

ELSE  IF  (NSERNS (REGION)  ER.  5)  THEN 

CALL  ELLIPTICAL  CYLINDER(REGION,  POINTCNT, 

COLOR,  CONUNIT,  ISUBDIV) 


ELSE  IF  (NSERNS (REGION)  ER.  8)  THEN 
CALL  SPHEROIO(R£GICN,  POINTCNT,  COLOR, 

CONUNIT,  ISUBDIV) 


ELSE  IF  (NSERNS (REGION)  .ER.  7)  THEN 
CALL  PROLATE  OBLATE_SPHEROIO(REGION, 

CONUNIT, 


POINTCNT,  COLOR, 
ISUBDIV) 


ELSE  IF  (NSERNS (REGION)  .ER.  8)  THEN 
CALL  PROLATE_OBLATE,SPHEROIO(REGION, 

CONUNIT, 


POINTCNT,  COLOR, 
ISUBDIV) 


ELSE  IF  (NSERNS (REGION)  .ER,  9)  THEN 
CALL  TOROIO(REGION,  POINTCNT,  COLOR, 

CON_UNIT,  ISUBDIV) 


ELSE  IF  (NSERNS (REGION)  ,ER.  10)  THEN 

CALL  RUAORILATERAL(REGION,  POINTCNT,  COLOR, 

CON  UNIT,  ISUBDIV) 


aSE  IF  (NSERNS (REGION)  .ER.  11)  THEN 

CALL  AXISYMMETRIC (REGION,  POINTCNT,  COLOR, 

CONUNIT,  ISUBDIV) 


aSE  IF  (NSERNS (RitlON)  ER.  12)  THEN 
CALL  TRIANGLE(REGION,  POINTCNT,  COLOR, 

CON_UNIT,  ISUBDIV) 

aSE  IF  (NSERNS (REGION)  .ER.  13)  THEN 
CALL  CONE(REGION,  POINTCNT,  COLOR, 

CON_UNIT,  ISUBDIV) 

END  IF 
END  DO 
END  DO 


CAU  CONNECT, ARROWS(CON_ UNIT,  POINTCNT) 
CLOSE(CON_UNIT) 


Coipleted  creating  the  shapes, 
and  Z  axes 

Rotated  coordinates  about  the  X,  Y, 

CALL  R0T3D(P0INTCNT,  AX,  AY, 

AZ) 

Write  the  rotated  coordinates  to  coordinate  file. 

TEVP  =  RUNI0(1:LEN)  //  ’.CRD’ 

OPEN  (UNIT=CRD_UNIT,FILE=TEMP,TYPE=’NEW’) 

DO  I  =  1,  POINTCNT 

WRITE(CRD,UNIT,«)  R0TATED(1 , I) ,  R0TATED(2, I) ,  R0TATFD(3,I) 
END  DO 

CALL  DRAW  ARROWS (CRD  UNIT,  POINTCNT,  AX,  AY,  AZ) 
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CLOSE(CRD_UNIT) 

WRITE(6,20e)  ■  ■ 

WRITE(6,200)  ’  ’ 

TYPE  »,  'Coordinate  data  has  been  stored  in  ’  //  TEMP 
TEHP  =  RUNID(1;LEN)  //  '.CON' 

TYPE  «,  'Connection  data  has  been  stored  in  '  //  TEMP 
TYPE  ♦ , ’  ' 

END 
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PLOTCHIEFO  -  written  by  D.  Mitchell  Bolling  Jr. 

TRI/TESSCO,  Inc. 

August  1989 

Purpose:  The  PLOTCHIEF  subroutine  is  a  graphics  routine  used  to  display 
three  dimensional  surfaces  defined  by  the  CHIEF  program. 

CHIEF  (Combined  Helmholtz  Integral  Equation  Formulation),  developed 
at  NPSC,  computes  the  acoustic  radiation  from  arbitrary-shaped  bodies 
PLOTCHIEF  creates  a  plot  file  which  can  be  viewed  and  rotated  before 
compiling,  linking  and  executing  the  CHIEF  driver.  PLOTCHIEF  should 
be  called  after  all  calls  to  LOSURR  are  made. 


Subroutines  used: 


RECTANGULAR_PLANE() 

aLIPTICAL_PLANE() 

ELLIPTICAL_CYLINDER() 

TOROIDO 

QUADRILATERAL  0 

TRIANGLEO 

AXISYMMETRICO 

CONEO 

CONNECT  ARROWS^) 


CIRCULARPLANEO 

CIRCULARCYLINDERO 

SPHEROIDO 

PROLATE_OBLATE_SPHEROID() 

AXISYMMETRICO 

QUAuRILATERALO 

TRIANGLEO 

R0T3DO 

DRAW  ARROWS  0 


Modification  Log: 


AUG  89  -  Added  QUADRILATERAL  0 ,  AXISYMMETRICO  and  TRIANGLEO. 

Also  modified  program  by  adding  R0T3D  so  that  the  user  can  rotate 
the  points  from  within  P'.OTCHIEF.  This  requires  that  the  points 
are  stored  in  an  array,  then  passed  to  RQT30,  which  writes  them  out 
to  a  file.  Before,  the  points  were  sent  straight  to  a  file. 


JAN  90  -  Added  CONNECT. ARROWS 0  and  DRAW.ARROWSO  in  order  to 
better  show  the  axis  orientation  during  rotations. 


JUL  90  -  Added  CONEO  geometry  to  PLOTCHIF. 


SUBROUTINE  PLOTC!IEF(RUNID,  NUMBLKS,  CSYMTYPE,  ISUBDIV, 
AX,  AY,  AZ) 


C 

c 

c 


rnRAMETER  PI  =  3. H 1595. 536 

PARAMETER  COLOR  USER  DEFINED  =  255  !  color  to  use  for  user  defined 
PARAMETER  cOLOR'TRANSFORMED  =  2  !  color  for  transformed  part 

PARAMETER  CRDUNIT  =  8  I  'un  number  of  coordinate  file 

PARAMETER  CONUNIT  =9  !  lun  number  of  connectivity  file 

PARAMETER  UXSREG  =  500  I  Maximum  number  of  surface  regions 

PARAMETER  MAXCOORDS  =  10000  !  m^.lmum  number  of  coordinates  to  plot 


Parameters  used: 


CHARACTER»(») 

INTEGER 

CHARACTER»3 

INTEGER 

REAL 

REAL 

REAL 


RUNID  !  title  of  CHIEF  run 

NUMBLKS  !  number  of  symmetry  blocks 

CSYMTYPE  !  symmetry  type  from  CHIEF 

ISUBDIV  !  if  1,  then  NSV  and  NSU  are  to  be  used 

AX  !  rotation  angle  in  degrees  about  x  axis 

AY  I  rotation  angle  in  degrees  about  y  axis 

AZ  !  rotation  angle  in  degrees  about  z  axis 


part 


C . - . - - - 

C  List  of  global  variables: 

C . - . - - - - - 

C  INTEGER  NSREG  I  number  of  CHIEF  surfaces 

C  INTEGER  NSEQNS(MXSREG)  I  types  of  CHIff  surfaces 

COMMON/SVALS/NSREG,NSEQNS(MXSREG) ,  SUL  (UXSREG)  ,SW(MXSREG) , 

.  SVL(MXSREG) ,SVU(MXSREG) , NSU (UXSREG) , NSV (MXSREG) , 

♦  CCS ( 10 , MXSREG) , TRNSS (3 , MXSREG) , IZAX (MXSREG) , 

♦  lORDU(MXSREG) ,IORDV(MXSREG),NCCEQS 

INTEGER  SYUTYP  !  typo  of  symmetry:  reflective  (1), 

I  rotational  (2),  none  (0) 
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INTEGER 

NBLKS 

!  number  of  symmetry  blocks 

!  for  reflective  symmetry: 

REAL 

SYMX 

!  if  -1,  then  reflects  about 

YZ 

p  1  ane 

REAL 

SYMY 

!  if  -1,  then  reflects  about 

XZ 

p  1  ane 

REAL 

SYMZ 

!  if  -1,  then  reflects  about 

XY 

p  lane 

!  for  rotational  symmetry: 

REAL 

SYMANG 

!  radians  to  rotate  object 

COMMON  / 

PL0T_INF0  /  SYMTYP,  NBLKS,  SYMX,  SYMY,  SYMZ, 

SYMANG 

REAL 

C00RD3D 

!  array  of  3d  points 

REAL 

ROTATED 

!  3d  points  after  rotations 

COMMON  /  R0T_INFQ  /  CaQRD30(3,  MAXCOORDS),  ROTATED (3, MAXCOORDS) 

C - - - - 

C  List  of  local  variables: 

C . . . . . . 

CHARACTER»S0  TEMP  !  a  tenporary  string  variable 

INTEGER  COLOR  I  nueber  associated  with  different  colors 

INTEGER  POINTCNT  !  keeps  track  of  next  point  to  be  created 

INTEGER  SHOWNORMALS  !  when  set,  uses  dotted  lines  for  inward  noreals 

INTEGER  REGION  I  counter  used  for  looping  through  each  region 

INTEGER  BLOCK  !  counter  used  for  looping  through  each  syeeetry 

!  block  (whether  reflective  or  rotational) 


C . - . 

C  Some  or  all  of  these  constants  are  euitiplied  by  each  data  point 

C . - . 

REAL  REFLECTX(8)  !  array  of  8  constants  for  reflecting  about  yi  plane 

REAL  REFLECTY(8)  !  array  of  8  constants  for  reflecting  about  xi  plane 

REAL  REFLECTZ(8)  !  array  of  8  constants  for  reflecting  about  xy  plane 


C . 

C  Below  is  a  table  that  indicates  the  nueber  of  constants  needed  for 
C  each  reflective  sysnetry  option.  Thus,  the  dashed  horizontal  line 
C  is  over  the  constants  needed  for  each  option. 

C . 

C  3  plane  symmetry  |< - >| 

C  2  plane  symmetry  |< - >|  j 

C  1  plane  symmetry  j< - >|  j  j 

C  III  I 

DATA  RERECTX  /  1,  -1,  1,  -1,  1,  -1,  1,  -1  / 

DATA  REFLECTY  /  1,  1,  -1,  -1,  1,  1,  -1,  -1  / 

DATA  REFLECTZ  /  1,  1,  1,  1,  -1,  -1,  -1,  -1  / 


C . 

C  Initial ize  variables 

C . - . . 

POINTCNT  =  0 
SYMX  =  1 
SYMY  =  1 
SYMZ  =  1 


If  NUMBLKS  is  positive,  then  display  of  normals  is  differentiated 
by  line  type.  Solid  for  positive  normals,  dashed  for  negative  normal. 


IF  (NUMBLKS  LT.  0)  THEN 
N8LKS  =  ABS (NUMBLKS) 
SHQW_NORMALS  =  0 
ELSE 

NBLKS  =  NUMBLKS 
SHOW_NCRMALS  =  1 
END  IF 

IF  (CSYMTYPE  EQ.  'REF')  THEN 
SYMTYP  =  1 


38 


NRL  Memorandum  Report  6709 


BLSE 

SYMTYP  =  2 
END  IF 


C-  - . - - - - - 

C  Open  the  connectivity  file.  <RUNID>.CON  will  contain  the  line 
C  drawing  data  in  the  format:  STARTPQINT  ENDPOINT  LINECOLOR 

C  Note:  The  actual  coordinate  data  will  be  written  to  (RUNID>.CRD 
C  in  the  RQT3D()  subroutine  located  near  the  end  of  this  program, 

C  since  the  coordinates  must  be  rotated  before  they  can  be  written  to 
C  a  f i le. 

C - - - - - - 

TEMP  =  RUNID  //  ’.CON’ 

OPEN  (UNIT=CON_UNIT,FILE=TEMP,TYPE=’NEI*) 


C . . . 

C  Generate  points  in  each  symmetry  block. 

C . . . . . . . 

DO  BLOCK  =  NBLKS,  1,-1 

IF  (S^YP  .EQ.  1)  THEN 

C . . - - - - - - 

C  For  current  block,  determines  which  planes  of  reflection 

C - - - 

SYMX  =  REFLECTX (BLOCK) 

SYMY  =  REFLECTY (BLOCK) 

SYMZ  =  REFLECTZ (BLOCK) 

ELSE  IF  (SYMTYP  EQ.  2)  THEN 


C  For  current  block,  determines  angle  to  rotate 

C . - . - . 

SYMANQ  =  (2  •  PI  »  (BLOCK  -  1)  /  NBLKS) 

END  IF 

C . 

C  Loop  through  all  the  types  of  surfaces  defined  and  call  the  proper 
C  shape  drawing  subroutine. 

C . 

DO  REGION  =  1,  NSREG 

C . - 

C  User  defined  block  will  be  a  different  color  from  the  computer 
C  generated  block(s)  If  IZAX  is  negative  (inward  normal)  then  make 
C  color  negative.  The  graphics  driver  should  draw  a  dotted  line  when 
C  colors  are  negative. 

C- . - . 

IF  (BLOCK  EQ.  1)  THEN 

IF  (SHOW_NORMALS  .El).  1)  THEN 

COLOR  =  ISIGN(COLOR_i;SER_DEFINED,  IZAX(REGION)) 

ELSE 

COLOR  =  COLOR_USER_DEFINED 
END  IF 
ELSE 

IF  (SHOI_NORUALS  .ER,  1)  THEN 

COLOR  =  ISIGN(COLOR_TRANSFORMED,  IZAX(REGION)) 
a.SE 

COLOR  =  COLOR_TRANSFORMED 
END  IF 
END  IF 

IF  (NSERNS (REGION)  ,ER.  1)  THEN 

CALL  RECTANGULAR_PLANE(REGION,  POINTCNT, 

•  COLOR,  CON_UNIT,  ISUBDIV) 

ELSE  IF  (NSEQNS (REGION)  ER.  2)  THEN 
CALL  CIRCULAR_PLANE(REGION,  POINTCNT, 

•  COLOR,  CONUNIT,  ISUBDIV) 

ELSE  IF  (NSERNS (REGION)  .ER.  3)  THEN 
CALL  aLIPTICAL^PLANE(REGION,  POINTCNT, 

•  COLOR,  CON  UNIT,  ISUBDIV) 

ELSE  IF  (NSERNS (REGIONl  ,ER.  A)  THEN 

CALL  CIRCUUR_CYLINDER(REGION,  POINTCNT, 
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COLOR,  CONUNIT,  ISUBDIV) 

ELSE  IE  (NSEQNS (REGION)  .EQ.  5)  THEN 

CALL  aLIPTICAL  CYLINDER(REGION,  POINTCNT, 

COLOR,  C0N_UNIT,  ISUBDIV) 

ELSE  IF  (NSEQNS (REGION)  .EQ.  6)  THEN 
CALL  SPHERQID(REGION,  POINTCNT,  COLOR, 

CON_UNIT,  ISUBDIV) 

ELSE  IF  (NSEONS (REGION)  .EO.  7)  THEN 

CALL  PRQLATE_OBLATE_SPHERQIO(REGION,  POINTCNT, 

COLOR,  CON_UNIT.  ISUBDIV) 

ELSE  IF  (NSEQNS (REGION)  EQ.  8)  THEN 

CALL  PR0LATE_08LATE  SPHEROIO(REGION,  POINTCNT, 

COLOR,  CON_UNIT,  ISUBDIV) 

ELSE  IF  (NSEONS (REGION)  ER.  9)  THEN 

CALL  TQROIO(REGIQN,  POINTCNT,  COLOR,  CON_UNIT,  ISUBDIV) 

ELSE  IF  (NSEQNS (REGION)  .ER.  10)  THEN 
CALL  0UAORILAT£RAL(REGION,  POINTCNT, 

COLOR,  C0N_UN1T,  ISUBDIV) 

ELSE  IF  (NSEONS (REGION)  ER.  11)  THEN 
CALL  AXISYMUETRIC (REGION,  POINTCNT, 

COLOR,  CON_UNIT,  ISUBDIV) 

ELSE  IF  (NSEONS (REGION)  ER.  12)  THEN 
CALL  TRIANGLE(R£GION,  POINTCNT, 

COLOR,  CON_UNIT,  ISUBDIV) 

aSE  IF  (NSEONS (REGION)  ER.  13)  THEN 
CALL  CONE(REGION,  POINTCNT, 

COLOR,  C0N_UNIT,  ISUBDIV) 

END  IF 
END  DO 
END  DO 

CALL  CONNECT_ARROWS(CON_UNIT,  POINTCNT) 

CLOSE(CON_UNIT) 


Conplefced  creating  the  shapes.  Rotate  coordinates  about  the  X,  Y, 
and  Z  axes. 


CALL  ROT30(POINTCNT,  AX,  AY,  AZ) 


C . - . . . . . - - - 

C  Write  the  rotated  coordinates  to  coordinate  file. 

C . - . . . . 

TEMP  =  RUNIO  //  ’.CRD’ 

OPEN  (UNIT=CRD_UN IT , FILE=TEJ(IP, TYPE=  ’ NEW ’ ) 

DO  I  =  1,  POINTCNT 

WRITE(CRD_UNIT,«)  ROTATED(l,I) ,  R0TATED(2,I) ,  ROTATED(3,I) 
END  DO 

CALL  DRAW_ARROWS(CRD_UNIT,  POINTCNT,  AX,  AY,  AZ) 
CLOSE(CRD_UNIT) 

WRITE  (6,.)  ’  ’ 

WRITE  (6,»)  ’Coordinate  data  has  been  stored  in  ’  //  TEMP 
TEMP  =  RUNID  //  ’  CON’ 

WRITE  (6,*)  'Connection  data  has  been  stored  in  ’  //  TEMP 
WRITE  (6,.)  ’  ’ 

RETURN 

END 
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LOAD_CID_DATA()  -  written  by  D.  Mitchell  Bolling  Jr. 

TRI/TESSCQ,  Inc. 

August  1989 

Purpose:  This  subroutine  was  designed  to  read  necessary  data  froe  a 
CIO  save  file  into  global  variables  used  by  PLOTCID. 

Subroutines  used;  None 

Modification  Log: 


SUBROUTINE  LOAD_CID_DATA(FILENAME) 

PARAMETER  MXSREG  =  500  I  Maxieui  nuober  of  surface  regions 

PARAMETER  MAXCOR  =  1000  '  Maxieu*  number  of  finite  element  nodes 

PARAMETER  PI  =  3.1415926S36 


List  of  parameters: 


CHAR.ACTERv(»)  FILENAME  !  name  of  CIO  data  file  to  load 
List  of  global  variables: 


INTEGER 

NSREG 

number  of  surfaces 

INTEGER 

NSEQNS (MXSREG) 

types  of  surfaces 

REAL 

SUL (MXSREG) 

lower  1 iaits  of  U 

REAL 

SUU (MXSREG) 

1 

upper  1 iaits  of  U 

REAL 

SVL (MXSREG) 

1 

lower  1  iaits  of  V 

REAL 

SVU (MXSREG) 

1 

upper  1  iaits  of  V 

REAL 

CCS (10, MXSREG) 

1 

constants  needed  for  CHIff  equations 

REAL 

TRNSS (3, MXSREG) 

1 

translation  froa  local  to  global  or  g 

INTEGER 

IZAX (MXSREG) 

global  axis  that  corresponds  to  local 
Z  axis 

COMMON/SVALS/NSREG, NSE9NS (MXSREG) , SUL (MXSREG) ,SUU(MXSREG) , 

SVL (MXSREG) ,SVU (MXSREG) ,NSU (MXSREG) ,NSV (MXSREG) , 
CCS (10 , MXSREG) , TRNSS (3 , MXSREG) , IZAX (MXSREG) , 
lORDU(MXSREG) , lOROV (MXSREG) .NCCEqS 


INTEGER 

SYMTYP 

!  type  of  syaaetry 

INTEGER 

NBLKS 

!  number  of  syaaetry  blocks 

!  for  reflective  syaaetry: 

REAL 

SYMX 

!  if  -1,  then  reflects  across 

YZ  plane 

REAL 

SYMY 

I  if  -1,  then  reflects  across 

XZ  plane 

REAL 

SYMZ 

!  if  -1,  then  reflects  across 

XY  plane 

!  for  rotational  syaaetry: 

REAL 

SYMANG 

!  radians  to  rotate  object 

COMMON  /  PLOT_INFO  /  SYMTYP,  NBLKS,  SYMX,  SYMY,  SYMZ,  SYMANG 


COMMON  /  CORD  /  COORDS (MAXCOR,  3) 

C - - 

C  List  of  local  variables. 

C . - . - - - - 

CHARACTER»8  COOR  FILE 
CHARACTER»3  COOR_EXT 
CHARACTER«55  COOR_FORMAT 
CHARACTERvB  B.EM^FILE 
CHARACTER*3  ELEM_EXT 
CHARACTERvSS  a.EM_FORMAT 

INTEGER  CCTEMP(10) 

C - - - - - 

C  The  following  variables  are  dusmy  variables  that  are  used  to 
C  read  past  unwanted  data  Some  of  thew  are  also  used  to  teaporarily 
C  hold  string  data  so  that  it  say  be  converted  into  nuaeric  data 
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INTEGER  lOUM 

tntcgcr  TDUM50(R0^ 

CHARACTER*3  C0UM3 
CHARACTER* 4  CDUM4 
CHARACTER*?  CDUM7 
CHARACTER* 8  CDUM8 
CHARACTER*9  CDUM9 
CHARACTER*58  CDUM58 
CHARACTER*!  COUM1_50(50) 
CHARACTER*2  ClXjM2_50(50) 
CHARACTER*8  CDUM8_100_3(100,3) 
CHARACTER*0  CDUM8_50_3(100,3) 
CHARACTER*10  CDUU10_ 50(50) 
CHARACTER*10  CDUM10_50  10(50,10) 
CHARACTER*20  CDUM20_50'(50) 


C . . . . - . . 

C  Open  the  CIO  interactive  driver's  save  file. 

C— . . . . 

OPEN (UNIT=20 , FILE=FILENA«E, STATUS= 'OLD' , 
-  F0RM=’ UNFORMATTED') 

READ (20)  CDUM8 
READ (20)  C0UM9 
READ (20)  C0UM58 
READ (20)  CDUU7 
READ (20)  CDUM7 
READ(20)  CDUM7 
READ (20)  CDUM7 
READ (20)  C0UM4 
READ (20)  CDUM4 


Depending  on  the  value  of  CDUU3,  SYMTYP  is  set  for  either  reflective 
or  rotational  syseetry. 


READ (20)  C0UM3 
IF  (C0UM3  .EQ.  'REF')  THEN 
SYMTYP  =  1 

B-SE  IF  (CDUM3  ER.  'ROT')  THEN 
SYMTYP  =  2 

asE 

SYMTYP  =  0 
END  IF 


Read  and  convert  the  string  C0UM3  into  a  nueber  and  store  in  NBLKS. 
NBLKS  »ill  be  either  1,  2  or  3  for  reflective  syeeetry,  representing 
the  nueber  of  planes,  or  nueber  of  blocks  for  rotational  syeeetry. 


READ (20)  C0UM3 

READ (UNIT=CDUM3 (1:3), FMT= ' (BNI3) ' )  NBLKS 


READ(20)  IDUM 
READ (20)  IDUM 

READ(20)  ((CDUM8_100  3(X,XX) ,X=1, 100) ,XX=1,3) 
READ(20)  101*1 


Reading  in  the  nueber  of  CHIEF  regions  defined 
READ (20)  NSREG 

READ (20)  (CDUM20_50(X),X=1,50) 


C . . . — 

C  Read  in  the  surface  type  for  each  CHIEF  region 

C . - . 

READ (20)  (NSERNS(X),X=1,B0) 


C 
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C  Reading  in  the  local-to-global  translations.  There  are  3  of  the* 
C  for  each  region. 


READ (20)  ( (CDUU8_60_3 (X , XX) , X=1 . 50) , Xy=l , 3) 

DO  X  =  1,  50 
DO  XX  =  1,  3 

READ(UNIT=CDUM8_50_3(X,XX) (1:8), 

FMT='(BNE8.0)’)  TRNSS(XX,X) 

END  DO 
END  00 


C - - - - - 

C  Read  in  IZAX  for  each  region.  IZAX  describes  the  global  orientation 
C  of  each  region. 

C - - - 

READ(20)  (CDUkl2_50(X),X=l,50) 

READ(20)  (CDUMl  50(X) ,X=1 , 50) 

DO  X  =  1,  50 

REAO('.'NIT=CDUM2_50(X)  (2:  2)  ,FMT=  ’  (BNIl)  ’)  IZAX(X) 

IF  (-OUM1_50(X)(1:1)  Eq.  ’-')  IZAX(X)  =  -iZAX(X) 

END  DO 

READ f 20)  (COUU2_50(X),X=1,50) 

READ(20)  (CDUU2_S0(X),X=1,50) 


C . . 

C  Read  in  the  lower  liait  of  U  for  all  the  regions 

C . . 

REAO(20)  (CDUU10  50(X),X=1.50) 

DO  X  =  1,  50 

READ(UNIT=CDUM10  50(X)(1:10), 

*  FMT=‘(BNE10.0)’)  SUL(X) 

END  DO 


Read  in  the  lower  liiit  of  V  for  all  the  regions 


REAO(20)  (COUM10_50(X),X=1,50) 

DO  X  =  1,  50 

READ(UNIT=C(XJU10_50  (X)  (1 : 10) , 

FMT='(BNE10.0)’)  SVL(X) 

END  00 


C . - . . 

C  Read  in  the  upper  liait  of  U  for  all  the  regions 

C . . 

READ(20)  (COUM10_50(X),X=1,50) 

DO  X  =  1,  50 

READ(UNIT=CDUM10_50(X) (1 : 10) , 

»  FMT=’(BNE10.0)’)  SUU(X) 

END  DO 


C . . . . 

C  Read  in  the  upper  liait  of  V  for  all  the  regions 

C . . . 

READ(20)  (CDUi(10_50(X),X=l,50) 

DO  X  =  1,  50 

READ(UNIT=CDOW10_50(X) (1:10) , 

*  FUT=’(BNE10.0)’)  SVU(X) 

END  DO 


C . - . . . . . -  - - - 

C  Read  in  the  constants  for  all  the  regions  used  in  CHIEF  subroutine  CCUMND 

C . . . . 

READ(20)  ((CDUM10_50_10(X,XX),X=1,50),XX=1,10) 

DO  X  =  1,  50 
DO  XX  =  1,  10 

REAO(UNIT=COUM10^60_10(X,XX) (1:10) , 

♦  FMT=’(BNE10  0)’,ERR=100)  CCS(XX,X) 

100  END  DO 

END  DO 
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c - - - - - 

C  Read  in  values  for  NSU,  the  number  of  subdivisions  in  the  U  direction 


REAO(20)  (CIXJ«2_50(X),X=1,50) 

DO  X  =  1,  50 

READ(UNIT=CDUW2_60(X) (1 :2) ,FMT=’ (8NI2) ’)  NSU(X) 
END  DO 


C - - - - - 

C  Read  in  values  for  NSV,  the  number  of  subdivisions  in  the  V  direction 

C - - - - - - - 

READ(20)  (CDUtt2_50(X),X=1.50) 

DO  X  =  1,  50 

READ(UNIT=:COUU2_50(X)  (1:2)  ,FHT=‘  (BN12)  ’)  NSV(X) 

END  DO 

REAO(20)  (CDUM1_50(X),X=1.50) 

RE>^D(20)  (CDUM10_50(X),X=1,S0) 

REXD(20)  (CDUM10_50(X),X=1,50) 


C - - - 

C  The  coordinate  and  element  files  are  used  rith  shapes  10,  11,  and  12. 

C . . - - - - - - - 

READ(20)  COQR  FILE 
READ (20)  C00P_EXT 
READ (20)  C0QR_FQRUAT 

READ (20)  ELEM  FILE 
READ (20)  ELEMIEXT 
READ (20)  ELEM_F0RMAT 


C 

C 

C 


Close  the  file. 
CLOSE (20) 


IF  ((NSREG  .EQ.  1)  .AND.  (NSE5NS(1)  .EQ.  10))  THEN 


If  there  is  only  one  surface,  and  the  surface  equation  is 
linear  or  quadratic  interpolation  over  a  quadrilateral,  then 
read  the  nodes  into  the  COORDS  array. 


CDUM58  =  COOR  FILE  //  ' . '  //  C00R_EXT 
OPEN(UNIT=20,FILE=CDgM58,TYPE=’OLD‘) 
lOUM  =  1 

CDUilSe  =’(■//  COOR  FORMAT  //  ’)’ 

200  REAO(20,FMT=COIAI58,END=22S)  (C00RDS(IDUM,X),X=1,3) 

IDUM  =  IDUM  ♦  1 
GOTO  200 
225  CLOSE(20) 

C-- . . . . . . 

C  Read  in  the  eleient  data  into  CCS(l)  through  CCS(8) . 

C . . . — . . 

CDUMS8  =  ELEM_FILE  //  //  a.EM_EXT 

OPEN (UNIT=20 , FILE=C0UM58 , TYPE= ’ OLD ’ ) 

NSREG  =  1 

CDOM58  ='('//  ELEM  format  //  ’)’ 

250  READ(20,FMT=COUM58,END^275)  (CCTEVlP(X)  ,X=1,8) 

DO  I  =  1,  8 

CCS(I, NSREG)  =  CCTEMP(I) 

END  DO 

NSEONS (NSREG)  =  10 
NSREG  =  NSREG  »  1 
GOTO  250 
275  CLOSE(20) 

NSREG  =  NSREG  -  1 

ELSE  IF  ((NSREG  EQ.  1)  AND.  (NSE5NS(1)  .ER.  11))  THEN 


C  If  there  is  only  one  surface,  and  the  surface  equation  is 
C  linear  or  quadratic  axisyisetric  interpolation,  then 
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C  read  the  nodes  into  the  COORDS  array, 

C . - - - 

rmiusn  -  rnno  pti  c  //  •  •  n  rnnp  pyt 

0  PEN (UNI T=20 . FILE=CDUW58 , TVpE= ’ OLD ' ) 

IDUM  =  1 

CDUM58  =  ■(’  //  COOR_FORMAT  //  ■)’ 

3ZS  READ(20.FMT=CDUM58,END=325)  (COORDS (IDUM, X) ,X=1, 2) 

lOUM  =  IDUW  •  1 
GOTO  oob 
325  CLOSE(20) 

C . . . - -  - - 

C  Read  in  the  element  data  into  CCS(l)  through  ^03(3). 

C- . - - - - - - - 

CDUM58  =  ELEM_FILE  //  //  ELai(_0(T 

0  PEN (UNIT=20 , F ILE=CDUM58 , TYPE= ’ OLD ’ ) 

NSREG  r  1 

CDUM58  =  ■('  //  EL0.(_FORMAT  //  ’)* 

350  READ(20,FMT=CDUM58,END=375)  (CCTEMP(X) ,X=1,3) 

DO  I  =  1,  3 

CCS (I, NSREG)  =  CCTEMP(I) 

PNH  00 

NSEqNS(NSREG)  =  11 
SVL (NSREG)  =  -PI  /  NBLKb 
SVU (NSREG)  =  PI  /  NBLKS 
NSREG  =  NSREG  *  1 
GOTO  350 
375  CLOSE(20) 

NSREG  =  NSREG  -  1 

a.SE  IF  ((NSREG  ER.  1)  .AND  (NSERNS(l)  .ER.  12))  THEN 

C . - . 

C  If  there  is  only  one  surface,  and  the  surface  equation  is 
C  linear  or  quadratic  interpolation  over  a  triangle,  then 
C  read  the  nodes  into  the  COORDS  array. 

C . 

CDUMS8  =  COOR  FILE  //  //  COOR  EXT 

OPEN (UN IT=20 , FILE=C0UMS8 , TYPE= ’ OLD ’ ) 

IDUM  =  1 

C0UM58  =  ’(’  //  COOR  FORMAT  //  ’)■ 

400  REAO(20,FMT=CDUM58,END=425)  (COOkOS(IDUM,X) ,X=1,3) 

IDUM  =  IDUM  ♦  1 
GOTO  400 
425  CLOSE(20) 

C . - . - - - 

C  Read  in  the  element  data  into  CCS(l)  through  CCS(6). 

C . . . - . - . - . - . . . 

COUU50  =  B.EM_FILE  //  ' .  ’  //  a.Bi(_EXT 
OPEN (UNIT=20 , FILE=CDUV58 , TYPE= ’ OLD ’ ) 

NSREG  =  1 

COWISa  =  ’(’  //  aEM_FORMAT  //  ')’ 

450  READ(20,FUT=COUM58,END=475)  (CCTEMP(X) ,X=1,6) 

DO  I  =  1,  6 

CCS (I, NSREG)  =  CCTEUP(I) 

END  DO 

NSERNS (NSREG)  =  12 
NSREG  =  NSREG  •  1 
GOTO  450 
475  CLOSE(20) 

NSREG  =  NSREG  -  1 
END  IF 

RETURN 

END 
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PUT  POINTO  -  written  by  D.  Mitchell  Bolling  Jr. 

TRI/TESSCO,  Inc. 

August  1989 

Purpose:  This  subprogras  takes  a  3D  point  and  translates  it  fro* 
the  local  coordinate  system  to  the  global  coordinate  system.  It 
then  reflects  or  rotates  the  point,  if  necessary  and  writes  it  to 
a  file. 

Subroutines  used;  none 
Modification  Log: 


SUBROUTINE  PUT_PQINT(IRG,  N.  X,  Y,  Z) 

PARAMETER  MXSREG  =  500  !  Maximum  number  of  surface  regions 

PARAMETER  MAXCOORDS  =  10000 


List  of  parameters  used: 


INTEGER  IRG 
INTEGER  N 

X,  Y,  2 


!  CHIEF  surface  region  identification  number 
!  current  coordinate  number 
I  Local  coordinates 


List  of  global  variables  used: 


REAL  TRNSS (3 , MXSREG)  !  translation  from  local  to  j,l:tal  origin 

INTEGER  IZAX(MXSREG)  !  global  axis  that  corresponds  to  local 

!  Z  axis 

COMMON/SVALS/NSREG.NSEGNS (MXSREG) , SUL  (MXSREG) ,SUU (MXSREG) , 

SVL (MXSREG) ,SVU(MXSREG) ,NSU(MXSREG) ,NSV (MXSREG) , 

CCS ( 10 , MXSREG) , TRNSS (3 , MXSREG) , I2AX (MXSREG) , 
IQROU(MXSREG) , lORDV (MXSREG) .NCCEHS 


INTEGER 

SYMTYP 

!  type  of  symmetry 

INTEGER 

N8LKS 

!  number  of  symmetry  blocks 

!  for  reflective  symmetry: 

REAL 

SYMX 

!  if  -1,  then  reflects  across 

YZ  plane 

REAL 

SYMY 

!  if  -1,  then  reflects  across 

XZ  plana 

REAL 

SYMZ 

!  if  -1,  then  reflects  across 

XY  plane 

!  for  rotational  symmetry: 

REAL 

SYMANG 

!  radians  to  rotate  object 

COMMON  / 

PL0T_INF0  /  SYMTYP,  NBLKS,  SYMX,  SYMY,  SYMZ,  SYMANG 

REAL 

C00RD3D 

!  array  of  3d  points 

REAL 

ROTATED 

!  3d  points  after  rotations 

COMMON  /  R0T_INF0  /  C00RD3D(3,  MAXCOORDS),  R0TAT£D(3, MAXCOORDS) 

C . . . . 

C  Local  variables  used. 

C . - . . . . . 


REAL  PX,  PY,  PZ  !  holds  global  coordinates  of  point 


C . - - - - - 

C  Translate  from  local  to  global  coordinates  according  to  IZAX. 

C - - - - 

IF  (IZAX(IRG)  .ER.  I)  THEN 
PX  =  Z  -  TRNSS (1,  IRG) 

PY  =  X  ♦  TRNSS (2,  IRG) 

PZ  =  Y  ♦  TRNSS(3,  IRG) 

ELSE  IF  (IZAX(IRG)  ,ER.  2)  THEN 
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PX  =  Y  ♦  TRNSS(1,  IRG) 

PY  =  Z  ♦  TRNSS(2,  IRG) 
rt  =  A  »  IRC) 

ELSE  IF  (IZAX(IRG)  .EQ.  3)  THEN 
PX  =  X  *  TRNSS(1,  IRG) 

PY  =  Y  ♦  TRNSS(2,  IRG) 

PZ  =  Z  *  TRNSS(3,  IRG) 

ELSE  IF  (IZAX(IRG)  .ER.  -1)  THEN 
PX  =  -Z  ♦  TRNSS(1,  IRG) 

PY  =  X  ♦  TRNSS(2,  IRG) 

PZ  =  -Y  .  TRNSS(3,  IRG) 

ELSE  IF  (IZAX(IRG)  ,Eq.  -2)  THEN 
PX  =  -Y  ♦  TRNSS(1.  IRG) 

PY  =  -Z  »  TRNSS(2.  IRG) 

PZ  =  X  ♦  TRNSS(3,  IRG) 

ELSE  IF  (IZAX(IRG)  .EQ.  -3)  THEN 
PX  =  X  .  TRNSS(1,  IRG) 

PY  =  -Y  .  TRNSS(2,  IRG) 

PZ  =  -Z  ■  TEiNSS(3,  IRG) 

END  IF 


C  Do  reflective  or  rotational  syeeetry  if  needed  on  the  GLOBAL  coordinates. 


IF  (SYMTYP  .ER.  1)  THEN 
PX  =  PX  »  SYMX 
PY  =  PY  «  SYUY 
PZ  =  PZ  »  SYMZ 
ELSE 

R  =  SqRT(PX*PX  ♦  PY»PY) 
IF  (R  .NE.  00)  THEN 
A  =  ATAN2(PY,  PX) 

A  =  A  ♦  SYUANG 
PX  =  R  «  COS (A) 

PY  =  R  ♦  SIN (A) 

END  IF 
END  IF 


C  Output  the  point  to  the  coordinate  file  for  plotting. 


N  =  N  ♦  1 

C00RD3D(1,  N)  =  PX 
COQR030(2,  N)  =  PY 
C00RD3D(3,  N)  =  PZ 
RETURN 
END 
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c - 

c 

RECTANGULAR  PLANEC)  -  written  by  0.  Mitchell  Bolling  Jr 

TRI/TESSCO,  Inc 
August  1989 

Purpose:  RECTANGULAR  PLANE  is  called  by  PLOTCID  or  PLOTCHIEF  It 
generates  the  points  and  lines  required  for  drawing  a  cylinder  or 
sect  ion  of  a  cy 1 inder . 

Subroutines  used:  PUT  PQINT() 

Modification  Log: 


SUBROUTINE  RECTANGULAR  PLANE(IRG,  N.  COLOR,  CON  UNIT,  SUBDIV) 

PARAMETER  MXSREG  =  500  !  Max  leu*  nueber  of  surface  regions 

C— - - - - - - - - 

C  List  of  subroutine  paraneters: 

C . . — - - - - - - - - 

INTEGER  IRG  !  surface  region  id  nuuber  (1  <  IRG  (  NSREG) 

INTEGER  N  !  current  coordinate  nuwber 

INTEGER  COLOR  !  number  representing  the  color  to  draw  plane 

INTEGER  CONUNIT  !  lun  number  for  connections  file 

INTEGER  SUBDIV  !  if  set,  then  use  NSU  and  NSV  for  subdivisions 


C . . 

C  List  of  global  variables: 


c 

REAL 

SUL (MXSREG) 

!  lower  1 imits  of  U 

C 

REAL 

SUU (MXSREG) 

!  upper  1 imits  of  U 

c 

REAL 

SVL (MXSREG) 

!  lower  limit?  of  V 

c 

REAL 

SVU (MXSREG) 

!  upper  1 imits  of  V 

c 

REAL 

CCS (10, MXSREG) 

!  constants  needed  for  CHIff^  equations 

COMMON/SVALS/NSREG.NSEQNS (MXSREG) , SUL (MXSREG) ,SUU (MXSREG) , 

SVL(MXSREJ) ,SVU(MXSREG) , NSU (MXSREG) , NSV (MXSREG) , 
CCS ( 10 , MXSREG) , TRNSS (3 , MXSREG) , IZ AX (MXSREG) , 
lORDU (MXSREG) , lORDV (MXSREG) .NCCERS 


C . 

C  List  of  local  variables: 

C . 

INTEGER  ICOUNT 
INTEGER  JCOUNT 
INTEGER  BEGINPOINT 
INTEGER  ENDPOINT 
INTEGER  NUMVPOINTS 
INTEGER  NUMUPOINTS 


!  counter  used  to  draw  segmenting 
!  counter  used  to  draw  segmenting 
!  used  for  connecting  segments 
!  used  for  connecting  segments 
!  number  of  points  used  when  drawing  U  lines 
!  number  of  points  used  when  drawing  V  lines 


REAL  I,  J 
REAL  X,  Y.  Z 
REAL  USTEP 
REAL  VSTEP 


!  scratch  variables 

!  current  calculated  x,  y,  z  values  (local) 
!  distance  between  each  U  segment  drawn 
!  distance  between  each  V  segment  drawn 


C . - - - 

C  If  SUBDIV  is  sot,  then  use  NSV  and  NSU  to  determine  segmenting, 
C  otherwise,  segmenting  defaults  to  4  for  U  and  V  direction. 

C . - - - - - . 

IF  (SUBDIV  .EQ.  1)  THEN 
NUMVPOINTS  =  NSV(IRG) 

NUMUPOINTS  =  NSU (IRG) 

a.SE 

NUMVPOINTS  =  4 
NUMUPOINTS  =  4 
END  IF 

USTEP  =  (SUU(IRG)  -  SUL(IRG))  /  NUMUPOINTS 
VSTEP  =  (SVU(IRG)  -  SVL(IRG))  /  NUMVPOINTS 
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C  Draw  paralell  lines  in  V  direction. 

C - 

Z  =  CCS(1,  IRG) 

X  =  SUL(IRG) 

DO  JCOUNT  =  0,  NUMUPQINTS 

Y  =  SVL(IRG) 

CALL  PUT_POINT(IRG,  N,  X,  Y,  Z) 

Y  =  SVU(IRG) 

CALL  PUT_PQINT(IRG,  N,  X,  Y,  Z) 
WRITE(CQN_UNIT,»)  N  -  1,  N,  COLOR 
X  =  X  .  USTEP 
END  DO 


C - 

C  Draw  lines  pareliel  to  U  direction. 

C - 

Y  =  SVL(IRG) 

DO  ICOUN^  0,  NUMVPOINTS 
X  =  jULCIRG) 

CALI.  PUT_POINT(!RG,  N,  X,  Y,  Z) 

X  =  SUU(IRG) 

CALL  PUT_PQINT(IRG,  N,  X,  Y,  Z) 
HRITECCON  UNIT,.)  N  -  1,  N,  COLOR 
Y  =  Y  .  VSTEP 
END  DO 

RETURN 

END 
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C-  - - 

c 

C  CIRCULAR^ PLANEO  -  written  by  D,  Mitchell  Bolling  Jr. 

C  TRI/TESSCO,  Inc. 

C  August  1989 

C 

C  Purpose:  CIRCULARPLANE  is  called  by  PLOTCID  or  PLOTCHIEF.  It 
C  generates  the  points  and  lines  required  for  drawing  a  circular 
C  planar  region  or  section  of  a  circular  planar  region 
C 

C  Subroutines  used:  PUT_PQINT() 

C 

C  Modification  Log: 

C 

C - - - 


SUBROUTINE  CIRCULAR_PLANE(IRG,  N,  COLOR,  CONUNIT,  SUBDIV) 

PARAMETER  PI  =  3.1415926536 

PARAMETER  MXSREG  =  500  !  Maxieun  number  of  surface  regions 

C - - - - - -  - 

C  List  of  subroutine  parameters: 

C - - - - - - - - - 

INTEGER  IRG  !  CHIEF  surface  region  id  number  (1  <  IRG  <  NSREG) 

INTEGER  N  !  current  coordinate  number 

INTEGER  COLOR  •  number  representing  color  to  draw  the  plane. 

INTEGER  CONUNIT  !  lun  number  for  connections  file 

INTEGER  SUBDIV  !  If  set,  then  use  NSU  and  NSV  for  suodi visions 


C . 

C  List  of  global  variables: 


c 

REnL 

SUL (MXSREG) 

!  lower  1 imits  of  U 

c 

REAL 

SUU (MXSREG) 

1  upper  1 imits  of  U 

c 

REAL 

SVL (MXSREG) 

!  lower  1 imits  of  V 

c 

REAL 

SVU (MXSREG) 

!  upper  1 imits  of  V 

c 

REAL 

CCS (10, MXSREG) 

!  constants  needed  for  CHia  equations 

COMMON/SVALS/NSREG, NSEONS (MXSREG) , SUL (MXSREG) .SUU (MXSREG) . 

♦  SVL(MXSREG) ,SVU (MXSREG) , NSU (MXSREG) , NSV (MXSREG) , 

♦  CCS ( 10 , MXSREG) , TRNSS (3, MXSREG) . IZAX (MXSREG) , 

♦  lORDU(MXSREG) ,  lORDV(MXSREG)  .NCCEIJS 

C . - . 

C  List  of  local  variables: 

C . - . . . . . 

INTEGER  ICOUNT  !  counter  used  to  draw  segmenting 

INTEGER  JCOUNT  !  counter  used  to  draw  segmenting 

INTEGER  BEGINPOINT  !  first  point  of  circular  plane 

INTEGER  ENDPOINT  !  last  point  of  circular  plane 

INTEGER  NUMVPOINTS  !  number  of  points  to  use  when  drawing  V  arcs 

INTEGER  NUMUFUINTS  !  number  of  points  to  use  when  drawing  U  lines 


REAL  I,  J 
REAL  X,  Y,  I 
REAL  VSTEP 
REAL  USTEP 


!  scratch  variable 

!  current  calculated  x,  y,  i  values  (local) 

!  angle  between  each  V  arc  point  drawn  (in  radians) 
'  distance  between  each  U  segment  drawn 


C . . - . . . . . - 

C  Determine  the  number  of  points  to  plot  for  arc  and  the  step  angle 
C  between  each  point.  For  a  full  2  *  PI  (380  degree)  arc,  48  points 
C  will  be  used.  The  minimum  number  of  points  used  is  4. 

C- - - - - - - - 

IF  (SUBDIV  .EQ.  1)  THEN 
NUMVPOINTS  =  4  »  NSV(IRG) 

NUMUPOINTS  =  NSU(IRG) 

asE 

T  =  24  0  «  (SVU(IRG)  -  SyL(IRG))  /  PI 
NUMVPOINTS  =  4  «  INT(I  /  4  0  ♦  99) 

NUMUPOINTS  =  1 
END  IF 
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VSTEP  =  (SVU(IRG)  -  SVL(IRG))  /  NUUVPQINTS 
USTEP  =  (SUU(IR0  -  SUL(IRG))  /  NUUUPOINTS 

J  =  SUL(IRG) 

BEGINPOINT  =  N  «  1 
DO  JCQUNT  =  0,  NUMUPQINTS 

C - - - - - 

C  Plot  the  first  point  o,"  the  outer  arc. 

C - - - . . 

X  =  J  .  G0S(SVL(IRG)) 

Y  =  J  »  SINfSVl^IRC)) 

Z  =  CCS(i,IRG) 

CALL  PUT  PQINT^IRG,  N,  X,  Y,  Z) 
ENDP0INT'=  N 


C  Plot  the  outer  arc. 


I  =  SVi.rRG)  ♦  "STEP 
DO  IlOUNT  =  1,  NUUVPOINTS 
),  =  J  .  CQS(I) 

Y  =  J  f  SIN(I) 

CALL  PUT_POINT(IRG,  N,  X,  Y,  Z) 
WRITb(CON  UNIT,*)  N  -  1,  N,  COLOR 
I  =  I  «  VSTEP 
END  DO 

J  =  J  *  USTEP 
END  DO 


C— . - . 

C  Subdivide  it  into  sections  to  show  solidity  and  exit. 

C . - . 

DO  ICOUNT  =  BEGINPOINT,  BEGINPOINT  ♦  NUMVPOINTS,  4 
WRITE(CON_UNIT,*)  ICOUNT,  ENDPOINT,  COLOR 
ENDPOINT  =  ENDPOINT  .  4 
END  DO 


REniRN 

END 
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C - - - - 

c 

C  ELLIPTICAL_PLANE()  -  written  by  D.  Mitchell  Bolling  Jr. 

C  TRI/TESSCO,  Inc. 

C  August  1969 

C 

C  Purpose:  ELLIPTICAL_PLANE  is  ca I  led  by  PUOTCID  or  PLOTCHIEF,  It 
C  generates  the  points  and  lines  required  for  drawing  a  elliptical 
C  planar  region  or  section  of  an  elliptical  planar  region. 

C 

C  Subroutines  used:  PUT_P0INT() 

C 

C  Modification  Log: 

C 

c - - - - - 


SUBROUTINE  aLIPTICAL_PLANE(IRG,  N,  COLOR,  C0N_UNIT,  SUBDIV) 
PARAMETER  PI  =  3.1416926536 

PARAMETER  MXSREG  =  500  !  Maxinui  nusber  of  surface  regions 

C  - . . . . . - - - 

C  List  of  subroutine  paraoeters: 

C . .  - - - - 

INTEGER  IRG  !  CHIEF  surface  region  id  nueber  (1  <  IRG  <  NSREG) 

INTEGER  N  !  current  coordinate  nuaber 

INTEGER  COLOR  !  number  representing  color  to  draw  the  plane 

INTEGER  CONUNIT  !  lun  number  for  connections  file 

INTEGER  SUBDIV  !  if  set,  then  use  NSU  and  NSV  for  subdivisions 


C . 

C  List  of  global  variables: 


c 

REAL 

SUL (MXSREG) 

!  lower  1 imits  of  U 

c 

REAL 

SUU (MXSREG) 

!  upper  limits  of  U 

c 

REAL 

SVL (MXSREG) 

!  lower  limits  of  V 

c 

REAL 

SVU (MXSREG) 

!  upper  1 imits  of  V 

c 

REAL 

CCS (10, MXSREG) 

!  constants  needed  for  CHIEF  equations 

COMMON/SVALS/NSREG, NSEONS (MXSREG) , SUL (MXSREG) ,SUU (MXSREG) , 

.  SVL (MXSREG) , SVU (MXSREG) , NSU (MXSREG) , NSV (MXSREG) , 

♦  CCS (10, MXSREG) , TRNSS (3, MXSREG), IZAX (MXSREG) , 

♦  lORDU(MXSREG) , lORDV (MXSREG) ,NCCEqS 

C . . 

C  List  of  local  variables: 

C . . . . 

INTEGER  ICOUNT  !  counter  used  to  draw  segmenting 

INTEGER  JCOUNT  !  counter  used  to  draw  segmenting 

INTEGER  BEGINPOINT  !  first  point  of  circular  plane 

INTEGER  ENDPOINT  !  last  point  of  circular  plane 

INTEGER  NUMVPOINTS  !  number  of  points  to  use  when  drawing  V  arc 

INTEGER  NUMUPOINTS  !  number  of  points  to  when  drawing  U  lines 

REAL  I,  J  !  scratch  variable 

REAL  X,  Y,  Z  !  current  calculated  x,  y,  z  values  (local) 

REAL  VSTEP  !  angle  between  each  V  arc  point  drawn  (in  radians) 
REAL  USTEP  !  distance  between  each  U  segment  drawn 


C . . . . . 

C  Determine  the  number  of  points  to  plot  for  arc  and  the  step  angle 
C  between  each  point.  For  a  full  2  «  PI  (380  degree)  arc,  48  points 
C  wi 1 1  be  used.  The  minimum  number  of  points  used  is  4. 

C . . . - 

IF  (SUBDIV  .EQ.  1)  THEN 
NUMVPOINTS  =  4  *  NSV(IRG) 

NUMUPOINTS  =  NSU(IRG) 

ELSE 

I  =  24.0  «  (SVU(IRG)  -  SVL(IRG))  /  PI 
NUMVPOINTS  =  4  »  INT(I  /  4.0  *  .99) 

NUMUPOINTS  =  1 
END  IF 
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VSTEP  =  (SVU(IRG)  -  SVL(IRG))  /  NUMVPOINTS 
USTEP  =  (SUU(IRG)  -  SUL(IRG))  /  NLMJPOINTS 

J  =  SUL(IRG) 

BEGINPQINT  =  N  ♦  1 
DO  JCOUNT  =  B,  NUMUPOINTS 

c  -  - .  - . - - - : - 

C  Calculats  the  first  point  of  outer  ellipse 
C - - - - - 

X  =  CCS(1,IRG)  »  J  *  COS(SVL(IRG)) 

Y  =  CCS(1,IRG)  »  SGRT(J  •  J  -  1.0)  •  SIN(SVL(IRG}) 
Z  =  CCS(2,IRG) 

CALL  PUT_POINT(IRG,  N,  X.  Y,  2) 

ENDPOINT  =  N 


C - - - 

C  Plot  the  outer  el  1  ipse 

C - 

I  =  GVL(IRG)  ♦  VSTEP 
DO  iCOUNT  =  1.  NUMVPOINTS 
X  =  CCS(1,IRG)  «  J  »  COS(I) 

Y  =  CCS(1,IRG)  .  SORT(J  .  J  -  10)  »  SIN(I) 
CALL  PUT_POINT(IRG,  N,  X.  Y,  Z) 
»RITE(CON_UNIT,»)  N  -  1.  N,  COLOR 
I  =  I  ♦  VSTEP 
END  DO 

J  =  J  ■.  USTEP 
END  DO 


Subdivide  the  surface  into  sections  to  show  solidity  and  exit. 


DO  ICOUNT  =  BEGINPOINT,  BEGINPOINT  ♦  NUMVPOINTS,  4 
WRITE(CON_UNIT,*)  ICOUNT,  ENDPOINT,  COLOR 
ENDPOINT  =  ENDPOINT  .  4 
END  DO 

RETURN 

END 
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CIRCULAR_CYLINDEfi()  -  written  by  D.  Mitchell  Bolling  Jr. 

TRI/TESSCO,  Inc, 

August  1989 

Purpose:  CIRCULAR  CYLINDER  is  called  by  PLOTCID  or  PLOKHIEF. 

It  generates  the  points  and  lines  required  for  drawing  a  cylinder 
or  section  of  a  cylinder. 

Subroutines  used:  PUT_P0INT() 

Modification  Log: 


SUBROUTINE  CIRCULAR_CYLINDER(IRG,  N,  COLOR,  C0N_UNIT,  SUBDIV) 
PARAMETER  PI  =  3.141S926536 

PARAMETER  MXSREG  =  500  !  Maxinuw  nuffiber  of  surface  regions 


List  of  subi„utine  paraneters: 


INTEGER  IRG  !  CHIEF  surface  region  id  nueber  (1  <  IRC  <  NSREG) 

INTEGER  N  !  coordinate  nueber  of  last  point  plotted 

INTEGER  COLOR  !  nueber  representing  color  to  draw  this  shape 
INTEGER  CONUNIT  !  lun  nueber  for  connections  file 
INTEGER  SUBDIV  !  if  set,  then  use  NSU  and  NSV  for  subdivisions 


List  of  global  variables: 


REAL 

SUL (MXSREG) 

!  lower  lie  its  of  U 

REAL 

SUU (MXSREG) 

!  upper  1 iaits  of  U 

REAL 

SVL (MXSREG) 

!  lower  1  in  its  of  V 

REAL 

SVU (MXSREG) 

!  upper  lie  its  of  V 

REAL 

CCS (10, MXSREG) 

!  constants  needed  for  CHIEF  equations 

COMMON/SVALS/NSREG, NSEQNS (MXSREG) , SUL (MXSREG) , SUU(MXSREG) , 

SVL(MXSREG) ,SVU(MXSREG) ,NSU(MXSREG) ,NSV(MXSREG) , 
CCS (10, MXSREG) , TRNSS (3, MXSREG) ,IZAX (MXSREG) , 
lORDU(MXSREG) , lORDV (MXSREG) .NCCEQS 


List  of  local  variables: 


INTEGER  ICOUNT 
INTEGER  JCOUNT 
INTEGER  BEGINPOINT 
INTEGER  ENDPOINT 
INTEGER  NUMVPOINTS 
INTEGER  NUilUPOINTS 


!  counter  used  to  draw  segeenting 
!  counter  used  to  draw  segeenting 
!  first  point  of  cylinder 
!  last  point  of  cylinder 

!  nueber  of  points  to  use  when  drawing  V  arcs 
!  nueber  of  points  to  use  when  drawing  U  lines 


REAL  I,  J 
REAL  X,  Y,  Z 
REAL  VSTEP 
REAL  USTEP 


!  scratch  variables 

!  current  calculated  x,y,z  values  (local) 
!  angle  between  each  V  arc  point  drawn 
!  distance  between  each  U  segeent  drawn 


C . - . - - - -  - 

C  Detereine  the  nueber  of  points  to  plot  for  each  arc  and  the  step  angle 
C  between  each  point.  For  a  full  2  ♦  PI  (380  degree)  arc,  48  points 
C  will  be  used.  The  einiiua  nueber  of  points  used  is  4. 

C . . . . . 

IF  (SUBDIV  Efl.  1)  THEN 
NUMVPOINTS  =  4  «  NSV (IRG) 

NUMUPOINTS  =  NSU(IRG) 

asE 

I  =  24.0  w  (SVU(IRG)  -  SVL(IRG))  /  PI 
NUMVPOINTS  =  4  ♦  INT(I  /  4.0  .  .99) 

NUMUPOINTS  =  1 
END  IF 

VSTEP  =  (SVU(IRG)  -  SVL(IRG))  /  NUMVPOINTS 
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USTEP  =  (SUU{IRG)  -  SUL{IRG))  /  NIMJPQINTS 

J  =  SUL(IRG) 

BEGINPOINT  =  N  ♦  1 
DO  JCOUNT  =  0,  NUMUPOINTS 

C - - - - - - - - - 

C  Calculate  the  first  point  of  the  circle  on  the  top  of  cylinder 
C - - - - - 

X  =  CCS(1,IRG)  ♦  COS(SVL(IRG)) 

Y  =  CCS(1,IRG)  «  SIN(SVL(IRG)) 

Z  =  J 

CALL  PUT_POINT(IPG,  N,  X,  Y,  Z) 

ENDPOINT  =  N 

C - - - - - 

C  Draw  the  rest  of  the  circle  on  the  top  of  the  cylinder 

C  -  - - - 

I  =  SVL(IRG)  ♦  VSTEP 
DO  ICOLNT  =  1,  NUMVPOINTS 
X  =  CCSU.IRG)  «  COS(I) 

Y  =  CCS(1,IRG)  ♦  SIN(I) 

CALL  PUT_PQINT(IRG,  N,  X,  Y,  Z) 

WRITE(CON_UNIT,«)  N  -  1,  N,  COLOR 
I  =  I  ♦  VSTEP 
END  DO 

J  =  J  •  USTEP 
END  DC 


C . - . 

C  Draw  vertical  lines  to  sieulate  a  3D  cylinder,  then  exit. 

C . - . . . - . 

DO  ICQUNT  =  BEGINPOINT,  BEGINPOINT  ♦  NUMVPOINTS,  4 
fRITE(CON  UNIT,*)  ICOUNT,  ENDPOINT,  COLOR 
ENDPOINT  =  ENDPOINT  ♦  4 
END  DO 

RETURN 

END 
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ELLIPTICAL_CYLINDER()  -  written  by  D.  Mitchell  Bolling  Jr. 

TRI/TESSCO,  Inc. 

August  1989 

Purpose:  ELLIPTICAL_CYLINDER  is  called  by  PLOTCID  or  PLOTCHIEF. 

It  generates  the  points  and  lines  required  fro  drawing  a  cylinder 
or  section  of  a  cylinder  where  circumference  is  elliptical  in  shape. 

Subroutines  used:  PUT_P0INT() 

Modification  Log: 


SUBROUTINE  ELLIPTICAL_CYLINOER(IRG.  N,  COLOR,  CON  UNIT,  SUBDIV) 
PARAMETER  PI  =  3.141S926S36 

PARAMETER  MXSREG  =  500  •  Maximum  number  of  surface  regions 


List  of  subroutine  parameters: 


INTEGER  IRG  !  CHIEF  surface  region  id  number  (1  (  IRG  <  NSREG) 

INTEGER  N  !  coordinate  number  of  last  point  plotted 

TNTPnFR  rni  QR  !  number  re,- resenting  color  to  draw  cylinder 
INTEGER  CONUNIT  I  lun  number  for  connections  file 
INTEGER  SUSOIV  !  if  set,  then  use  NSU  and  NSV  for  subdivisions 


List  of  global  variables: 


REAL 

SUL (MXSREG) 

!  lower  1 imits  of  U 

REAL 

SUU (MXSREG) 

!  upper  limits  of  U 

REAL 

SVL  (MXSREG) 

!  lower  limits  of  V 

REAL 

SVU (MXSREG) 

!  upper  1 imits  of  V 

REAL 

CCS (10, MXSREG) 

!  constants  needed  for  CHIEF  equations 

COMMON/SVALS/NSREG,NSEqNS (MXSREG), SUL (MXSREG), SUU(MXSREO), 

♦  SVL (MXSREG) , SVU (MXSREG) , NSU (MXSREG) , NSV (MXSREG) , 

♦  CCS(10, MXSREG) ,TRNSS(3, MXSREG) , IZAX(MXSREG) , 

.  lORDU(MXSREG) , lOROV(MXSREG) ,NCCEHS 

C . - . . . 

C  List  of  local  variables; 

C . - - - 

INTEGER  ICOUNT  I  counter  used  to  draw  segmenting 

INTEGER  JCOUNT  !  counter  used  to  draw  segmenting 

INTEGER  BEGINPOINT  !  first  point  of  cylinder 

INTEGER  ENDPOINT  !  last  point  of  cyl inder 

INTEGER  NUMVPOINTS  !  number  of  points  to  use  when  drawing  V  arc 

REAL  I,  J  !  scratch  variables 

REAL  X,  Y,  Z  !  current  calculated  x,  y,  z  values  (local) 

REAL  VSTEP  !  angle  between  each  V  arc  point  drawn 
REAL  USTEP  !  distance  between  each  U  segment  drawn 


C- . - . - - - - - - — 

C  Determine  the  number  of  points  to  plot  for  each  arc  and  the  step  angle 
C  between  each  point.  For  a  full  2  «  PI  (380  degree)  arc,  48  points 
C  will  be  used.  The  minimum  number  of  points  used  is  4. 

C- . . . . . 

IF  (SUBDIV  .Eq.  1)  THEN 
NUMVPOINTS  =  4  »  NSV (IRG) 

NUMUPOINTS  =  NSU(IRG) 

ELSE 

I  =  24.0  »  (SVU(IRG)  -  SVL(IRG))  /  PI 
NUMVPOINTS  =  4  »  INT(I  /  4  »  .99) 

NUMUPOINTS  =  1 
END  IF 

VSTEP  =  (SVU(IRG)  -  SVL(IRG))  /  NUJVPOINTS 
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USTEP  =  (SUU(IRG)  -  SUL(IRO))  /  NUUUPOINTS 

J  =  SUL(IRG) 

BEGINPOINT  =  N  ♦  1 
DO  JCOUNT  =  0,  NUMUPOINTS 

C  - . - - - 

C  Calculate  the  first  point  of  an  ellipse  of  cylinder. 

C . - . - . - - - 

X  =  CCS(5,IRG)  «  CDS(SVL(IRG)) 

Y  =  CCS(6,IRG)  »  SIN(SVL(IRG}) 

Z  =  J 

CALL  PUT_POINT(IRG,  N,  X,  Y,  Z) 

ENDPOINT  =  N 


C - - - - - 

C  Draw  the  rest  of  the  ellipse  of  cylinder. 

C - - - - - - - 

I  =  SVL(IRG)  ♦  VSTEP 
DO  ICCUNT  =  1,  NUMVPOINTS 
X  =  CCS(S,IRG)  *  COS(I) 

V  =  CCS (6, IRC)  •  SIN  (I) 

CALL  PUT_POINT(IRG,  N,  X,  Y,  Z) 
WRITE(CON_UNIT,«)  N  -  1,  N,  COLOR 
1=1*  VSTEP 
END  DO 

J  =  J  *  USTEP 
END  DO 


C . - . 

C  Draw  extra  vertical  lines  to  show  solidity,  then  exit. 

C . - . 

DO  ICOUNT  =  BEGINPOINT,  BEGINPOINT  .  NUMVPOINTS,  4 
»RITE(CON  UNIT,*)  I,  ENDPOINT,  COLOR 
ENDPOINT  i  ENDPOINT  .  4 
END  DO 

RETURN 

END 
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SPHEROIOO  -  written  by  0.  Mitchell  Bolling  Jr, 

TRI/TESSCO,  Inc. 

August  1989 

Purpose:  SPHEROID  is  called  by  PLDTCID  or  PLOTCHIEF  in  order  to 
generate  the  points  and  lines  for  drawing  a  SPHEROID  or  section 
of  a  SPHEROID. 

Subroutines  used:  PUT_P01NT() 

Modification  Log: 

June  1990  -  PLDICID  would  not  fully  close  circles  that  were 
slightly  less  than  PI  for  a  SPHEROID  exaaple.  Deterained  that 
it  was  due  to  floating  point  roundoff  error  within  the  increaental 
loops.  Changed  indexes  froa  REAL  to  INTEGER. 


SUBROUTINE  SPHEROID(IRG,  N,  COLOR,  C0N_UNIT,  SUBDIV) 


List  of  subroutine  paraaeters: 


INTEGER  IRG  I  CHIEF  surface  region  id  nuaber  (1  <  IRG  <  NSREG) 
INTEGER  N  !  coordinate  nuaber  of  next  point  to  be  drawn 

INTEGER  COLOR  !  nuaber  associated  with  color  to  draw  spheroid 

INTEGER  C0N_UNIT  I  lun  nuaber  for  connections  file 
IirEuEfi  ‘■iipniv  !  if  set,  then  use  NSU  and  NSV  for  subdivisions 

PARAMETER  PI  =  3.1415926536 

PARAMETER  MXSREC  =  500  !  Maxiaua  nuaber  of  surface  regions 


List  of  global  variables  used. 


REAL 

SUL(MXSREG) 

!  lower  liaits  of  U 

REAL 

SUU(MXSREG) 

!  upper  1  iaits  of  U 

REAL 

SVL(MXSREG) 

!  lower  1 iaits  of  V 

REAL 

SVU(MXSREG) 

!  upper  liaits  of  V 

REAL 

CCS(10,MXSREG3 

!  constants  needed  for  CHIEF  equations 

COMMON/SVALS/NSREG,NSEI}NS(MXSREG)  ,SUL(MXSREG)  ,SUU(MXSREG) , 

SVL(MXSREG) ,SVU(MXSREG) ,NSU(MXSREG) ,NSV(MXSREG) , 
CCS(10,MXSREG) ,TRNSS(3,MXSREG) ,IZAX(MXSREG) , 
lORDU  CMXSREG) ,  lOROV  (MXSREG) ,  NCCqS 


C . - - - - - . 

C  List  of  local  variables: 

C— . . . 

INTEGER  NUMUPOINTS  !  nuaber  of  points  used  when  drawing  U  arcs 

INTEGER  NUWPOINTS  !  nuaber  of  points  used  when  drawing  V  arcs 

INTEGER  UCOUNT  !  counter  for  U  angle 

INTEGER  VCOUNT  I  counter  for  V  angle 


REAL  USTEP 
REAL  VSTEP 
REAL  X,  Y,  Z 
REAL  U 
REAL  V 
REAL  T 


!  angle  between  each  point  around  arc  in  U  direction 
!  angle  between  each  point  around  arc  in  V  direction 
!  current  calculated  x,  y  and  z  coordinates 
!  current  angle  corresponding  to  current  UCOUNT 
!  current  angle  corresponding  to  current  VCOUNT 
I  teaporary  variable 


C . . - . .  . 

C  Deteraine  the  nuaber  of  points  to  plot  for  each  arc  and  the  step  angle 
C  between  each  point.  For  a  full  2  »  PI  arc,  48  points  will  bo  used, 

C  An  arc  of  less  than  PI  /  6  will  use  4  points. 

C . - .  - . . . 

IF  (SUBDIV  .EO.  1)  THEN 
NUMUPOINTS  =  4  *  NSU(IRG) 

NUMVPOINTS  =  4  »  NSV(IRG) 
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B.SE 

T  =  24.0  ♦  (SUU(IRG)  -  SUL(IRG))  /  PI 
NUMUPOINTS  =  4  »  INT{T  /  4.0  ♦  .99) 

T  =  24.0  »  (SVU(IRG)  -  SVL(IRG))  /  PI 
NUyVPOINTS  =  4  ♦  INT(T  /  4.0  «  .99) 

END  IF 

USTEP  =  (SUU(IRG)  -  SUL(IRG))  /  NUMUPOINTS 
VSTEP  =  (SVU(IRG)  -  SVL(IRG))  /  NUUVPOINTS 


C  - - - - - 

C  Draw  circles  in  U  direction  roughly  every  PI/8  radians  (30  degrees). 
C  If  SUBDIV  is  set,  then  draw  NSU(IRG)  circles  in  the  U  direction. 

C . . - - - - - - 

M  -  ri  B  rTpr\ 

DO  UCOUNT  NUMUPOINTS,  4 

C . . - - - 

C  Plot  the  f i rst  point. 

C— - - - 

X  =  CCF(1,IRG)  »  SIN(U)  ♦  COS(SVL(IRG)) 

Y  =  CfS(l,IRG)  «  SIN(U)  »  SIN(SVL(IRG)) 

Z  =  C.;S(1,IRG)  »  CQS(U) 

CALL  PUT_PQINT(IRG,  N,  X,  Y,  Z) 

V  =  SVL(IRG)  ♦  VSTEP 

DO  VCOUNT  =  1,  NUMVPQINTS 

C . - . — . . 

C  Connect  a  line  to  the  next  point. 

C . . . 

X  =  CCS(l.IRG)  *  SINfil)  ,  co^w 

Y  =  CCS(1,IRG)  *  SIN(U)  «  SIN(V) 

CALL  PUT_POINT(IRG,  N,  X,  Y,  Z) 

WRITE(CON_UNIT,»)  N  -  1,  N,  COLOR 

C . - . 

C  Draw  another  point  on  the  circle  in  the  U  direction. 

C . - . - . 

V  =  V  ♦  VSTEP 
END  DO 


C . - 

C  Move  down  and  start  drawing  the  next  circle. 

C— . 

U  =  U  •>  USTEP  »  4 
END  DO 


C . - . - . . . . . 

C  Draw  circles  in  V  direction  roughly  every  PI/6  radians  (30  degrees) 
C  If  SUBDIV  is  set,  then  draw  NSV(IRG)  circles  in  the  U  direction. 

C . . - - - - - 

V  =  SVL(IRG) 

DO  VCOUNT  =  0,  NUUVPOINTS,  4 


C - - - - - - 

C  Plot  first  point  of  circle  in  V  direction. 

C- .  - - - 

X  =  CCS(1,IRG)  »  SIN(SUL(IRG))  ♦  COS(V) 
Y  =  CCS(1,IRG)  »  SIN(SUL(IRG))  *  SIN(V) 
Z  =  CCS(1,IRG)  «  COS(SUL(IRG)) 

CALL  PUT_POINT(IRG,  N,  X,  Y,  Z) 

U  =  SUL(IRG)  ♦  USTEP 
DO  UCOUNT  =  1,  NUMUPOINTS 


C - - - - 

C  Connect  a  line  to  the  next  point. 

C- - - 

X  =  CCS(1,IRG)  *  SIN(U)  ♦  COS(V) 

Y  =  CCS(1,IRG)  «  SIN(U)  *  SIN(V) 

Z  =  CCS(I,IRG)  »  COS(U) 

CALL  PUT  POINT(IRG,  N,  X,  Y,  Z) 
»RITE(CQN  UNIT,*)  N  -  1,  N,  COLOR 


C - 

C  Move  down  the  circle  in  V  direction 
C - - - 
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U  =  U  -  USTEP 
END  DO 


C - - - 

C  Move  across  and  start  drawing  the  next  circle. 

C - - - 

V  =  V  .  VSTEP  *  4 
END  DO 

RETURN 

END 
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C - - - - - - 

c 

C  PnOLATE  OBLATE  SPHEROIL»(;  -  written  by  D.  Mitchell  Bolling  Jr, 

C  '  TRI/TESSCO,  Inc. 

C  August  1989 

C 

C  Purpose:  PROLATE_aBLATE_SPHEROID  is  called  by  PLOTCID  and  PLOTCHIEF  in 
C  order  to  generate  the  points  and  lines  for  drawing  a  prolate/oblate 
C  spheroid  or  section  of  a  prolate/oblate  spheroid. 

C 

C  SuOrouLint;  used:  PUT_POINT() 

C 

C  Modification  Log: 

C  June  1990  -  PLOTCID  would  not  fully  close  circles  that  we'-e 
C  slightly  less  than  PI.  Determined  that  It  was  due  to  floating 
C  point  roundoff  error  within  the  incremental  loops.  Changed  them 
C  from  REAL  to  INTEGER 
C 

c - - - - - 


SUBROUTINE  PROLATE_OBLATE_SPHEROID(IRG,  N,  COLOR, 

CON  UNIT,  SUBDIV) 


INTEGER  IRG  !  CHIEF  surface  region  id  number  (1  <  IRG  <  NSREG) 

INTEGER  N  !  coordinate  number  of  next  point  to  be  drawn 

INTEGER  COLOR  !  f  associated  with  color  to  draw  prolate/oblate  spheroid 

INTEGER  CON  UNIT  !  lun  number  for  connections  file 

INTEGER  SUBDIV  !  If  set,  then  use  NSU  and  NSV  for  subdivisions 


PARAMETER  PI  =  3.1415926536 

PARAMETER  MXSREG  =  500  !  Maximum  number  of  surface  regions 


of  global 

variables  used: 

REAL 

SUL (MXSREG) 

!  lower 

imits  of  U 

REAL 

SUU(MXSREG) 

•  upper 

imits  of  U 

REAL 

SVL (MXSREG) 

!  lower 

imits  of  V 

REAL 

SVU(MXSREG) 

!  upper 

is  its  of  V 

REAL 

CCS  (10, MXSREG) 

!  constants  needed  for  CHIEF  equations 

REAL 

NSU  (MXSREG) 

*  number  of  subdivisions  in  U  direction 

REAL 

NSV (MXSREG) 

1  number 

of  subdivisions  in  V  direction 

COMMaN/SVALS/NSREG,NSEQNS(MXSREG),SUL(MXSREG),SUU(MXSREG), 

SVL (MXSREG) ,SVU (MXSREG) , NSU (MXSREG) , NSV (MXSREG) , 
CCS ( 10 , MXSREG) , TRNSS (3 , MXSREG) , IZAX (MXSREG) , 

I OROU (MXSREG) , lORDV (MXSREG) ,NCCEhS 


C - - - - - 

C  List  of  local  variables: 

C - - - - 

INTEGER  NUMUPOINTS  !  number  of  points  used  when  drawing  U  arcs 

INTEGER  NUMVPOINTS  I  number  of  points  used  when  drawing  V  arcs 

INTEGER  UCOUNT  !  counter  for  U  angle 

INTEGER  VCOUNT  !  counter  for  V  angle 


REAL  USTEP 
REAL  VSTEP 
REAL  X,  Y,  Z 
REAL  U 
REAL  V 
REAL  T 


angle  between  each  point  around  arc  in  U  direction 
angle  between  each  point  around  arc  in  V  direction 
current  calculated  x,  y,  and  z  coordinates  (local) 
current  angle  corresponding  to  current  UCOUNT 
current  angle  corresponding  to  current  VCOUNT 
temporary  variable 


C - - - - - 

C  If  SUBDIV  is  set,  then  number  of  points  becomes  4  X  number  of 
C  subdivisions,  otherwise  use  the  default  subdivisions  and  determine 
C  the  number  of  points  to  plot  for  each  arc  and  the  step  angle 
C  between  each  point  For  a  full  2  •  PI  arc,  48  points  will  be  used 
C  An  arc  of  less  than  PI  /  6  will  use  4  points. 

C - - - 

IF  (SUBDIV  EQ  1)  THEN 
NUMUPOINTS  =  4  *  NSU(IRG) 

NUMVPOINTS  =  4  .  NSV(IRG) 
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ELSE 

T  =  24.0  *  (SUUCIRG)  -  SUL(IRG))  /  PI 
nUMUPOINTS  =  4  »  INT(T  /  4.0  *  .99) 

T  =  24  0  »  (SVU(IRG)  -  SVL(IRG))  /  PI 
NUMVPQINTS  =  4  »  INT(T  /  4.0  ♦  .99) 

END  IF 

USTEP  =  (S'JU(IRG)  -  SUL  (IRQ))  /  NUMUPOINTS 
VSTEP  =  (SVU(IRG)  -  SVL(IRG))  /  NUMVPQINTS 


C - - - 

C  Draw  circles  in  U  direction  roughly  every  PI  /  6  radians  (30  degrees). 
C  If  SUSDIV  is  set,  then  draw  NSU(IRG)  curcles  in  the  U  direction 
C - 

U  =  SUL(IRG) 

DO  UCOUNT  =  0,  NUMUPOINTS,  4 


C -  - 

C  Plot  the  f i rst  point. 

C - 

X  =  CCS(5,IRG)  .  SIN(U)  •  COS(SVL(IRG)) 

Y  =  CCS(5,IRG)  *  SIN(U)  v  SIN(SVL(IRGh 
Z  =  CCS(3,IRG)  «  COS(U) 

CALL  PUT_POINT(IRG,  N,  X,  Y,  Z) 

V  =  SVL(IRG)  •.  VSTEP 

DO  VCOUNT  =  1,  NUUVPOINTS 


C . — - 

C  Connect  a  line  to  the  next  point. 

C . - . - . 

X  =  CCS(5,IRG)  *  SIN(U)  «  COS(V) 

Y  =  CCS(S,IRG;  »  SIN(U)  «  SIN(V) 
CALL  PUT_POINT(IRG,  N,  X,  Y.  Z) 
WRITE(CGN  UNIT,»)  N  -  1.  N,  COLOR 


C . - . 

C  Draw  another  point  on  the  circle  in  the  U  direction 

C . 

V  =  V  ♦  VSTEP 
END  DO 

C . - . 

C  Move  down  and  start  drawing  the  next  circle. 

C . . 

U  =  U  .  USTEP  »  4 
END  DO 


C- . . .  - . . . . . . 

C  Draw  circles  in  V  direction  roughly  every  PI  /  6  radians  (30  degrees). 

C . - . - . — - - 

V  =  SVL(IRG) 

00  VCOUNT  =  0,  NUMVPQINTS,  4 

C . . . . - . . . . . 

C  Plot  first  point  of  circle  in  V  direction 

C- .  - - - - - 

X  =  CCS(5,IRG)  «  SIN(SUL(IRG))  «  COS(V) 

Y  =  CCS(5,IRG)  *  SIN(SUL(IRG))  *  SIN(V) 

Z  =  CCS(3,IRG)  «  C0S(SUL(IRG)) 

CALL  put_po:nt(irg,  N,  X,  Y,  Z) 

U  =  SX(IRG)  -  USTEP 
DO  UCOUNT  =  1,  NUMUPOINTS 


C - - - - - 

C  connect  a  line  to  the  next  point. 

C--- - - - - 

X  =  CCS(5,IRG)  «  SIN(U)  »  COS(V) 

Y  =  CCS(5,IRG)  »  SIN(U)  ♦  SIN(V) 

Z  =  CCS(3,IRG)  «  COS(U) 

CALL  PUT_PQINT(IRG,  N,  X,  Y,  Z) 
IRITE(C0N  UNIT,*)  N  -  1,  N,  COLOR 


C 
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C  Move  down  the  circle  in  V  direction. 

C - 

U  =  U  «  USTEP 
END  DO 

C - 

C  Move  across  and  start  drawing  the  next  circle. 

C - 

V  =  V  .  VSTEP  .  4 
END  DO 

RETURN 

END 
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TQRQID()  -  written  by  D.  Mitchell  Bolling  Jr. 

TRI/TESSCO,  Inc. 

August  1989 

Purpose:  TCROIO  is  called  by  PLOTCID  or  PLOTCHIEF  in  order  to 
generate  the  points  and  lines  for  drawing  a  toroid  or  section 
of  a  toroid 

Subroutines  used;  PUT  P0INT() 

Modification  Log: 

June  1990  -  PLOTCID  would  not  fully  close  circles  that  were 
slightly  less  than  PI  for  a  spherical  example.  Determined  that 
it  was  due  to  floating  point  roufidoff  error  within  the  incremental 
loops.  Changed  them  from  REAL  to  INTEGER 


SUBROUTINE  TQRaiD(IRG,  N,  COLOR,  C0M_JNIT,  SUBOIV) 


C . - .  -  . . . . . - . - 

C  List  of  subroutine  parameters: 

C . - . 

INTEGER  IRG  !  CHIEF  surface  region  id  number  (1  <  IRG  <  NSREG) 
INTEGER  N  •  coordinate  number  of  next  point  to  be  drawn 

INTEGER  COLOR  !  number  associated  with  color  to  draw  TOROID 

INTEGER  CONUNIT  !  lun  number  for  connections  file 
INTEGER  SUBOIV  I  if  set,  then  use  NSU  and  NSV  for  subdivisions 

PARAMETER  PI  =  3.1415926536 

PARAMETER  MXSREO  =  500  !  Maximum  number  of  surface  regions 

C . . - . 

C  List  of  global  variables: 


c 

REAL 

SUL (MXSREG) 

!  lower 

limits  of  U 

c 

REAL 

SUU (MXSREG) 

!  upper 

limits  of  U 

c 

REAL 

SVL (MXSREG) 

!  lower 

limits  of  V 

c 

REAL 

SVU (MXSREG) 

!  upper 

limits  of  V 

c 

REAL 

CCS (10, MXSREG) 

!  constants  needed  fo'  ChiiEF  equations 

COMMON/SVALS/NSREO, NSEqNS(MXSREG) , Sa (MXSREG) , SUU(MXSREG) , 

-  SVL (MXSREG) , SVU (MXSREG) , NSU (MXSREG) , NSV (MXSREG) , 

♦  CCS  00 ,  MXSREG) ,  TRNSS  (3 ,  MXSREG) ,  IZ  AX  (MXSREG) , 

«  lOROU(MXSREG) ,IORDV(MXSREG),NCCEqS 

C - - - 

C  List  of  local  variables: 

C - - - - - - 

INTEGER  NUMUPOINTS  !  number  of  points  used  when  drawing  U  arcs 

INTEGER  NUMVPOINTS  !  number  of  points  used  when  drawing  V  arcs 

INTEGER  UCOUNT  I  counter  for  U  angle 

INTEGER  VCOUNT  !  counter  for  V  angle 

REAL  USTEP  !  angle  between  each  point  around  arc  in  U  direction 

REAL  VSTEP  '  angle  between  each  point  around  arc  in  V  direction 

REAL  X,  Y,  Z  !  current  calculated  x,  y  and  z  coordinates  (local) 
REAL  U  I  current  angle  corresponding  to  current  UCOUNT 

REAL  V  f  current  angle  corresponding  to  current  VCOUNT 

REAL  T  I  temporai'y  variable 


C 

C 

C 

C 

C 

C 

r 


If  SUBOIV  is  set,  then  the  number  of  points  becomes  4  «  number  of 
subdivisions,  otherwise  use  the  default  subdivisions  and  determine 
the  number  of  points  to  plot  for  each  arc  and  the  step  angle 
between  each  point.  For  a  full  2  «  PI  arc,  48  points  will  be  used 
An  arc  of  IcjS  than  PI  /  6  will  use  4  points. 


IF  (SUBOIV  ER  1)  THEN 
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NUMUPOINTS  =  A  •  NSU(IRG) 

NLiMVPDINTS  =  4  »  NSV(IRG) 

ELSE 

r  =  24.0  »  (SUU(IRG)  -  SUL(IRG))  /  PI 
NUMUPOINTS  =  4  »  INT(T  /  4.0  .  .99) 

T  =  24.0  .  (SVU(IRG)  -  SVL(IRG))  /  PI 
NUWVPQINTS  =  4  .  INT(T  /  4.0  ♦  .99) 

END  IF 

USTEP  =  (SUU(IRG)  -  SUL(IRG))  /  NUMUPOINTS 

VSTEP  =  (SVU(IRG)  -  SVL(IRG))  /  NUMVPOINTS 


C - - - 

C  Draw  the  large  circles  in  U  direction  roughly  every  PI/8  radians 
C  (30  degree.s)  . 

C - 

U  =  SUL(IRG) 

DC  UCOUNT  =  0,  NUMUPOINTS,  4 


C - 

C  Plot  the  f  I  ."it  point. 

C - - 

X  =  '.CS(3,IRG)  V  COS(SVL(IRG))  / 

.  {CCS(1,IRG)  -  CQS(U)) 

Y  =  CCS(3,IRG)  «  SIN(SVL(IRG))  / 

♦  (CCS(l,IRa)  -  COS(U)) 

Z  =  CCS(3,IRG)  «  CCS(4,IRG)  •  SIN(U)  / 

♦  (CCS(1,IRG)  -  COS(U)) 

CALL  PUT_POINT(IRG,  N,  X,  Y,  Z) 

V  =  SVL(IRG)  ♦  VSTEP 

DD  VCOUNT  =  1,  NUMVPOINTS 


^ - 

C  Connect  a  line  to  the  next  point 


X  r  CCS(3,IRG)  »  COS(V)  / 

*  (CCS(1,IRG)  -  COS(U)) 

''  =  CCS(3,IPG1  •  SIN(V)  / 

♦  (CCS(1,IRG)  -  COS(U)) 

C.4LL  PUT  ?aiNT(IRG,  N,  X,  Y,  Z) 

*ri:e(con_unit,*)  n  -  i,  n,  color 

C . -  . . - . . - . 

C  Draw  arot''er  sent  on  the  circle  in  the  U  direction. 

C . - . . . . . - . 

V  =  V  •  VSTEP 
E.NO  CO 


C- - - - - 

C  Move  down  ard  start  drawing  the  next  circle. 
C - 

U  =  U  -  USTEP  .  4 
END  DO 


C - 

C  Draw  the  snail  circles  m  V  direction  roughly  every  PI  /  6  radians 
C  (38  degrees) 

C - 

V  SVL(!RG) 

DO  VCOUNT  =  0,  NUMVPOINTS,  4 


c - 

C  Plot  first  point  of  circle  in  V  direction. 

C - 

X  =  CCS(3,IRG)  .  cns(v)  / 

•  {CCS(1,IRC)  -  C0S(SUL(IRG))) 
Y  =  CCS(3,IRG)  .  SIN(V)  / 

(CCS(1,!RG)  -  COS(SUL(IRG))) 
r  -  CCS(3,IRG)  .  CCS(4,IRG)  »  S1N(5UL(IRG) )  / 

•  (CCS(’,IRG)  -  CCSfSUL(IRG))) 
CALL  PUT  PQINT(IRO,  N,  X,  Y,  Z) 

U  5UL(iRG)  •  i^SIEP 


f>5 


o  o  r>  o  <->  o 


Siders  and  Bolling 


DO  UCOUNT  =  1,  NUVUPOINTS 


Connect  a  line  to  the  next  point. 


X  =  CCS (3, IRQ)  «  COS(V)  / 

(CCS(1,IRG)  -  CQS(U)) 
Y  =  CCS(3,IRG)  «  SIN(V)  / 

(CCS(1,IRG)  -  COS(U)) 
Z  =  CCS(3,IRG)  »  CCS(4,IRG)  »  SIN{U)  / 
(CCS(1,IRG)  -  COS(U)) 
CALL  PUT  PDINT(IRG,  N,  X,  Y,  Z) 
WRITE(CGN_UNIT,0  N  -  1,  N,  COLOR 


Draw  another  point  on  the  circle  in  the  U  direction. 


U  =  U  .  USTEP 
END  DO 


C - 

C  Hove  across  and  start  drawing  the  next  circle. 

C - - 

V  =  V  ♦  VSTEP  »  A 
END  DO 

RETURN 

END 
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QUADRILATERAL ()  -  written  by  D.  Mitchell  Bolling  Jr. 

TRI/TESSCO,  Inc. 

August  1989 

Purpose:  QUADRILATERAL  is  called  by  PLQTCID  or  PLOTCHIEF.  It 
generates  the  points  and  lines  required  for  drawing  a  quadrilateral 

Subroutines  used;  PUT_PQINT{) 

Modification  Log: 


SUBROUTINE  QUADRILATERAL(IRG,  N,  CQLORN.  CON_UNIT.  SUBDIV) 

PARAMETER  MXSREG  =  506  I  Maxieue  number  of  surface  regions 

PARAMET'L-<  M^XCOR  =  1000  I  Max  mu*  number  of  finite  element  nodes 


C  i-ist  of  subroutine  paraieters: 

C - - 

INTEGER  IRQ  !  surface  region  id  number  (1  <  IRG  <  NSREG) 

INTEGER  N  !  current  coordinate  number 

INTEGER  COLQRN  !  number  representing  the  color  to  draw  plane 

INTEGER  CONUNIT  !  lun  number  for  connections  file 

INTEGER  SUBDIV  !  not  used  in  this  routine 

C . 

C  List  of  global  variables: 

C . - . 

C  REAL  CCS (10, MXSREG)  !  constants  needed  for  CHIEF  equations 

COMHON/SVALS/NSREG , NSEQN3 (MXSREG) , SUL (MXSREG) , SUU (MXSREG) , 

.  SVL (MXSREG) , SVU (MXSREG) ,NSU (MXSREG) ,NSV (MXSREG) , 

.  CCS (10 , MXSREG) , TRNSS (3 , MXSREG) , IZAX (MXSREG) , 

.  lORDU(MXSREG) ,IOROV(MXSREG) .NCCEQS 

COMMON/CORO/COORDS (MAXCOfi, 3) 

C . . . - . 

C  List  of  local  variables: 

C . - .  - - - - 

REAL  X,  Y,  Z  !  fork  variables 

INTEGER  COLOR  !  holds  color  and  line  type  (negative  for  dotted) 

C - - - . - . . . . . . 

C  if  IZAX  is  negative,  then  make  color  negative  to  show  inward  normal, 

C  The  graphics  driver  should  draw  a  dotted  line  when  colors  are  negative, 
C - 

COLOR  =  SIGN(COLORN,IZAX(IRG)) 

C - 

C  starting  at  node  1,  draw  lines  to  5,  2,  6,  3,  7,  4,  8  and  back  to  1. 

C 

C  1  < - 8  < - 4 

C  I 

C  V  I 

C  5  7 

C  I 

C  V  I 

c  2 - >  6 - >  3 

C 

C - - 

X  =  COORDS (CCS (1,  IRG),  1) 

Y  =  C00RDS(CCS(1,  IRG),  2) 

Z  =  COORDS  (CCS (1,  IRG) ,  3) 

CALL  PUT_P0INT(IRG,  N,  X,  Y,  Z) 

IF  (CCS(S,IRG)  NE  0)  THFN 
X  =  COORDS(CCS(5,  IRG),  1) 

Y  =  CQ0RUS(CCS(6,  IRG),  2) 


67 


Siders  and  Bolling 


Z  =  COOROS(CCS(5,  IRG).  3) 

cAa  PIT  n.  t.  ■’) 

»RITE(CON_UNIT,»)  H  -  1.  N,  COLOR 
END  IF 

X  =  COORDS (CCS (2,  IRG],  1) 

Y  =  COORDS (CCS (2,  IRG),  2) 

Z  =  COORDS (CCS (2,  IRG),  3) 

CALL  PUT  P0INT(IRG,  N,  X,  Y,  Z) 
fRITE(CON_UNIT,»)  N  -  I.  N,  COLOR 

IF  (CCS(5,IRG)  NE.  0)  THEN 
X  =  C00RDS(CCS(6,  IRG),  1) 

Y  =  COORDS (CCS (6,  IRG),  2) 

Z  =  C00RDS(CCS(6,  IRG),  3) 

CALL  PUT_POINT(IRG,  N,  X,  Y,  Z) 
WRITE(CON_UNIT,»)  N  -  1,  N,  COLOR 

END  IF 

X  =  C00RDS(CCS(3,  IRG),  1) 

Y  =  COORDS (CCS (3,  IRG),  2) 

Z  =  COORDS (CCS (3,  IRG),  3) 

CALL  PUT  POINT(IRG,  N,  X,  Y,  1) 
fRITE(CON_UNIT,*)  N  -  1,  N,  COLOR 

IF  (CCS(5,IRG)  .NE.  0)  THEN 
X  =  C00RDS(CCS(7,  IRG),  1) 

Y  =  COORDS (CCS (7,  IRG),  2) 

Z  =  COORDS (CCS (7,  IRG),  3) 

CALL  PUT_POINT(IRG,  N,  X,  Y,  Z) 
fRITE(CON_UNIT,»)  N  -  1,  N,  COLOR 

END  IF 

X  =  COORDS (CCS (4,  IRG),  1) 

Y  =  COORDS (CCS (4,  IRG),  2) 

Z  *  COORDS (CCS (4,  IRG),  3) 

CALL  PUT  P0INT(IRG,  N,  X,  Y,  Z) 
»RITE(CON_UNIT,«)  N  -  1,  N,  COLOR 

IF  (CCS(S,IRG)  NE.  0)  THEN 
X  =  COORDS (CCS (8,  IRG),  1) 

Y  =  COORDS (CCS (8,  IRG),  2) 

Z  =  COORDS (CCS (8,  IRG),  3) 

CALL  PUT  POINT(IRG,  N,  X,  Y,  Z) 
IRITE(CON_UNIT,«)  N  -  1,  N,  COLOR 

END  IF 

X  =  COORDS(CCS(l,  IRG),  1) 

Y  =  COORDS (CCS (1,  IRG),  2) 

Z  =  COORDS (CCS (1,  IRG),  3) 

CALL  PUT_POINT(IRG,  N,  X,  Y,  Z) 
WRITE(CON_UNIT,»)  N  -  1,  N,  COLOR 

RETURN 

END 
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AXISYUMETRICO  -  written  by  0.  Mitchell  Bolling  Jr 

TRI/TESSCO,  Inc. 

August  1989 

Purpose:  AXISYMMETRIC  is  ca I  led  by  PLOTCIO  or  PLOTCHIEF  It 
generates  the  points  and  lines  required  for  drawing  an  ax  i  syemetr  i  c  . 

Subroutines  used:  PtJT_PQINT() 

Modification  Log: 


SUBROUTINE  AXISYMMeTRIC(IRG,  N,  COLOR,  CONUNIT,  SUBOIV) 

PARAMETER  MXSREG  =  500  !  Maxiwuw  nueber  of  surface  regions 

FARAMETcd  MAXCOR  =  1000  !  Maxinuw  number  of  finite  element  nodes 

PARAMET^Jl  PI  =  3.141S926536 

C - 

C  List  of  subroutine  parameters: 

C- - - - - - - - 

INTEGER  IRG  !  surface  region  id  number  (1  <  IRG  <  NSREG) 

INTEGER  N  !  current  coordinate  number 

INTEGER  COLOR  I  number  representing  the  color  to  draw  plane 

INTEGER  CONUNIT  I  lun  number  for  connections  file 

INTEGER  SUBOIV  I  not  used  in  this  routine 

C .  - . 

C  List  of  global  variables: 

C . - . 

C  REAL  CCS (10, MXSREG)  !  constants  needed  for  CHIEF  equations 

C  REAL  COORDS (MAXCOR, 3)  !  holds  finite  element  coordinates 

COMMON/SVALS/NSREG.NSEQNS (MXSREG) , SUL (MXSREG) ,SUU (MXSREG) , 

♦  SVL (MXSREG) , SVU (MXSREG) ,NSU (MXSREG) ,NSV (MXSREG) , 

♦  CCS (10, MXSREG) ,TRNSS (3, MXSREG) ,IZAX (MXSREG) , 

♦  lORDU(MXSREG) , lORDV(MXSREG) ,NCCEQS 


COMMON/CORD/COOROS (MAXCOR, 3) 

C . - .  - 

C  List  of  local  variables: 

C . . . — . - . - . 

INTEGER  J  I  scratch  counter 

INTEGER  NIAIPOINTS  !  number  of  points  to  use  when  drawing  V  arc 

REAL  X,  Y,  Z,  R  !  work  variables 
REAL  I  I  angle  counter 

REAL  STEP  I  angle  between  each  V  arc  point  drawn 


C - - - - - 

C  Determine  the  number  of  points  to  plot  for  each  arc  and  the  steo  angle 

C  between  each  point.  For  a  full  2  ♦  PI  (360  degree)  arc,  48  points 

C  will  be  used.  The  minimum  number  of  points  used  is  4 

C - - - - - 

I  =  24.0  *  (SVU(IRG)  -  SVL(IRG))  /  PI 
NUMPOIVTS  =  4  »  INT(I  /  4.0  ♦  .99) 

STEP  =  (SVU(IRG)  -  SVL(IRG))  /  NUMPOINTS 

C . - - - - - 

C  Draw  first  points. 

C- . - . - - - 

R  =  C00R0S(CCS(1,  IRG),  1) 

Z  =  COORDS (CCS (1,  IRG),  2) 

X  =  COS(SVL(IRG))  *  R 
Y  -  SIN(SVL(IRG))  V  R 
CALL  PUT_POINT(IRG,  N,  X,  Y,  Z) 

R  =  C00RDS(CCS(2,  IRG),  1) 

Z  =  COOROS(CCS(2,  IRG),  2) 

X  =  COS(SVL(IRG))  »  R 


69 


Siders  and  Bolling 


Y  =  SIN(SVL(IRG))  »  R 

CAa  PUT_POINT(IRG,  N.  X.  Y,  Z) 


C  If  the  third  node  is  zero,  then  CCS(1,IRG)  and  CCS(2,IRG)  are  end 
C  nodes,  otherwise,  CCS(2,IRG}  is  a  midnode. 

C - - - 

IF  (CCS(3,IRG)  .NE.  0)  THEN 
R  =  COORDS (CCS (3,  IRG) ,  1) 

Z  =  COORDS (CCS (3,  IRG),  2) 

X  =  COS(SVL(IRG))  *  R 

Y  =  SIN(SVL(IRG))  »  R 

CALL  PUT_ POINT (IRG,  N,  X,  Y,  2) 

END  IF 


C - - - - - 

C  Draw  circles  perpendicular  to  the  Z  axis. 
C - 

I  =  SVL(IRG)  .  STEP 
DO  J  =  1,  NUMPOINTS 

R  =  COORDS (CCS (1,  IRG),  1) 

Z  =  COORDS (CCS (1,  IRG),  2) 

X  =  COS(I)  •  R 

Y  =  SIN(I)  .  R 

CALL  PUT_POINT(IRG,  N,  X,  Y,  Z) 

R  =  COOROS(CCS(2,  IRG),  1) 

Z  =  COORDS (CCS (2,  IRG),  2) 

X  =  COS(I)  X  R 

Y  =  SIN(I)  *  R 

CALL  PUT  P0INT(IRG,  N,  X,  Y,  Z) 


C . - . . 

C  If  the  third  node  is  zero,  then  CCS(1,IRG)  and  CCS(2,IRG)  are  end 
C  nodes,  otherwise,  CCS (2, IRG)  is  a  eidnode. 

C . - . 

IF  (CCS(3,IRO)  .NE.  0)  THEN 
R  =  C00RDS(CCS(3.  IRG),  1) 

Z  =  C00RDS(CCS(3,  IRG),  2) 

X  =  COS(I)  *  R 
Y  =  SIN(I)  «  R 

CALL  PUT_POINT(IRG,  N,  X,  Y,  Z) 

C . 

C  Connect  the  3  points  with  the  previous  3  points. 

C . - . . . 

WRITE(CON_UNIT,«)  N  -  5,  N  -  2,  COLOR 
IRITE(CON_UNIT,»)  N  -  ■«,  N  -  I,  COLOR 
WRITE(C0N  UNIT,*)  N  -  3,  N,  COLOR 
a.SE 

c . . . . .  . 

C  or  connect  the  2  points  with  the  previous  2  points. 

C . . - . - . - - - 

fRITE(CON_UNIT, »)  N  -  3,  N  -  1,  COLOR 
WRITE(CON_UNIT,«)  N  -  2,  N,  COLOR 
END  IF 
1=1*  STEP 
END  DO 


C - - - 

C  Draw  cross  sections  parallel  to  the  Z  axis  to  show  solidity 
C - 

I  =  SVL(IRG) 

DO  J  =  0,  NUMPOINTS,  A 
R  =  C00RDS(CCS(1,  IRG),  1) 

Z  =  COORDS (CCS (1,  IRG),  2) 

X  =  COS(I)  *  R 
Y  =  SIN(I)  *  R 

CALL  PUT_POINT(IRG,  N,  X,  Y,  Z) 

R  =  COORDS (CCS (2,  IRG),  1) 
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Z  =  C00RDS(CCS(2,  IRG) ,  2) 

X  =  COS(I)  ♦  R 
Y  =  SIN(I)  ♦  R 

CALL  PUT_POINT(IRG,  N,  X,  Y,  Z) 
»RITE(CON_UNIT,»)  N  -  1,  N,  COLOR 


C . - . - - - 

C  If  the  third  node  is  zero,  then  CCS(1,IRG)  and  CCS(2,IRG)  are  end 
C  nodes,  otherwise,  CCS(2,IRG}  is  a  lidnode. 

C - - - - 

IF  (CCS (3. IRG)  NE.  0)  THEN 
R  =  C00RDS(CCS(3,  IRG),  1) 

Z  =  COORDS (CCS (3,  IRG),  2) 

X  =  COS(I)  «  R 
Y  =  SIN(I)  *  R 

CALL  pi;t_point(irg,  N,  X,  Y,  Z) 

WRITE(CON_UNIT,»)  N  -  1,  N,  COLOR 
END  IF 

I  =  I  ♦  STEP  «  A 
END  DO 

RETURN 

END 
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TRIANGLEO  -  written  by  0.  Mitchell  Bolling  Jr. 

TRI/TESSCO,  Inc 
August  1989 


Purpose:  TRIANGLE  is  called  by  PLOTCID  or  PLOTCHIEF.  It 
generates  the  points  and  lines  required  for  drawing  a  triangle. 

Subroutines  used;  PUT_P0INT() 

Modification  Log; 


SUBROUTINE  TRIANGLE(IRG,  N,  COLOR,  C0N_UNIT,  SUBDIV) 

PARAMETER  MXSREG  =  580  *  Maxinuw  number  of  surface  regions 

PARAMETER  MAXCQR  =  1000  !  Maximum  number  of  finite  element  nodes 


List  of  subroutine  parameters: 

INTEGER  IRG  ! 

INTEGER  N  ! 

surface  region  id  number  (1  < 
current  coordinate  number 

IRG  <  NSREG) 

INTEGER  COLOR  ! 

number  representing  the  color 

to  draw  plane 

INTEGER  CON  UNIT  ! 
INTEGER  SUBDIV  ! 

lun  number  for  connections  fi 
not  used  in  this  subroutine 

le 

List  of  global  variables: 

C  REAL  CCS (10, MXSREG)  I  constants  needed  for  CHIEF  equations 

CQMMQN/SVALS/NSREG, NSEONS (MXSREG) , SUL (MXSREG) , SUU(MXSREG) , 
r  SVL (MXSREG), SVU (MXSREG), NSU (MXSREG), NSV (MXSREG) , 

♦  CCS (10 , MXSREG) , TRNSS (3, MXSREG) , IZAX (MXSREG) , 

♦  lORlXJ(MXSREG) ,  lOROV(MXSREG)  .NCCEQS 

C0IM0N/CaRD/C00RDS(MAXC0R,3) 


C . 

C  List  of  local  variables: 

C . 

REAL  X,  Y,  Z  !  Work  variables 


C  Starting  at  node  1,  draw  lines  to  4,  2,  5,  3,  6  and  back  to  1. 

C - - - - - 

X  =  C00RDS(CCS(1,  IRG),  1) 

Y  =  C00RDS(CCS(1,  IRG),  2) 

Z  =  COORDS (CCS (1,  IRG),  3) 

CALL  PUT_POINT(IRG,  N,  X,  Y,  Z) 

IF  (CCS(4,IRG)  ,NE.  0)  THEN 
X  =  COORDS (CCS (4,  IRG),  1) 

Y  =  COORDS (CCS (4,  IRG),  2) 

Z  =  COORDS (CCS (4,  IRG),  3) 

CAa  PUT_POINT(IRG,  N,  X,  Y,  Z) 

WRITE(CON_UNIT,»)  N  -  1,  N,  COLOR 

END  IF 

X  =  C00RDS(CCS(2,  IRG),  1) 

Y  =  C00RDS(CCS(2,  IRG),  2) 

Z  =  COORDS (CCS (2,  IRG),  3) 

CALL  PUT_POINT(IRO,  N,  X,  Y,  Z) 

WRITE(CON_UNIT,»)  N  -  1,  N,  COLOR 

IF  (CCS(4,IRG)  NE.  0)  THEN 
X  =  COOROS(CCS(5,  IRG),  1) 

Y  =  COORDS (CCS (5,  IRG),  2) 

Z  =  C00RDS(CCS(5,  IRG),  3) 

CALL  PUT„POINT(IRG,  N,  X,  Y,  Z) 

IRITE(CON_UNIT,»)  N  -  1,  N,  COLOR 


72 


NRL  Memorandum  Report  6709 


END  IF 

X  =  COORDS (CCS (3,  IRQ).  1) 

Y  =  COORDS (CCS (3,  IRO).  2) 

Z  =  COORDS (CCS (3,  IRQ),  3) 

CALL  PUT_POINT(IRC,  N,  X,  Y,  2) 
WRITE(CON_UNIT,«)  N  -  1,  N,  COLOR 

IF  (CCS (4, IRC)  .NE.  0)  THEN 
X  =  C00RDS(CCS(6,  IRQ),  1) 

Y  =  C00RDS(CCS(6,  IRQ),  2) 

Z  =  COORDS (CCS (6.  IRQ),  3) 

CAU  PUT_POINT(IRG,  N,  X.  Y,  Z) 
WRITE(CON_UNIT..)  N  -  1,  N,  COLOR 
END  IF 

X  =  COORDS (CCS (1,  IRQ),  1) 

Y  =  COORDS (CCS (1,  IRG),  2) 

Z  =  COORDS (CCS (1,  IRQ),  3) 

CALL  PUT_PaiNT(IRG,  N,  X,  Y,  Z) 

*rite(con_un:t,»)  n  -  i,  n,  color 

RETURN 

END 


73 


r>ooor»  ooo  oor><^r><^r>o  or»o  oooor>oor>ooooor>o 


Siders  and  Bolling 


C0NE()  -  written  by  0.  Mitchell  Bolting  Jr. 

TRI/TESSCD,  Inc. 

July  1990 

Purpose;  CONE  is  cslled  by  PLOTCID  or  PLOTCHIEF. 

It  generates  the  points  and  lines  required  for  drawing  a  cone 
or  section  of  a  cone. 

Subroutines  used:  PUT_PQINT() 

Modification  Log. 


SUBROUTINE  CONE(IRO,  N.  COLOR,  CONUNIT,  SUBDIV) 

PARAMETER  PI  =  3.1415926533 

PARAMETER  MXSREG  =  500  f  Maximus  number  of  surface  regions 


List  of  subroutine  paraseters; 


INTEGER  IRO  !  CHIff  surface  region  id  nuaber  (1  <  IRG  <  NSREG) 

INTEGER  N  I  coordinate  nuaber  of  last  point  plotted 

INTEGER  COLOR  I  nuaber  representing  color  to  draw  this  shape 
INTEGER  CONUNIT  !  lun  nuaber  for  connections  file 
INTEGER  SUBDIV  I  if  set,  then  use  NSU  and  NSV  for  subdivisions 


List  of  global  variables: 


REAL 

SUL (MXSREG) 

!  lower 

1  i a i ts  of  U 

REAL 

oUU (MXSREG) 

!  upper 

1  iaits  of  U 

REAL 

SVL (MXSREG) 

!  lower 

1 iaits  of  V 

REAL 

SVU (MXSREG) 

'  upper 

1 iaits  of  V 

REAL 

CCS (10, MXSREG) 

!  constants  needed  for  CHIEF  equations 

COMMON/SVALS/NSREG,NSECNS (MXSREG) , SUL (MXSREG) ,SUU (MXSREG) , 

SVL (MXSREG) , SVU (MXSREG) , NSU (MXSREG) , NSV (MXSREG) , 
CCS ( 10 , MXSREG) , TRNSS (3 , MXSREG) , IZAX (MXSREG) , 
lORDU (MXSREG) , lOROV (MXSREG) ,NCCEUS 


List  of  local  variables: 


INTEGER  VCOUNT 
INTEGER  UCOUNT 
INTEGER  BEGINPOINT 
INTEGER  ENDPOINT 
INTEGER  NUUVPOINTS 
INTEGER  NUMUPOINTS 


!  counter  used  to  draw  segmenting 
!  counter  used  to  draw  segmenting 
!  f i rst  point  of  ci rc le 
!  last  point  of  circle 

!  nuaber  of  points  to  use  when  drawing  V  lines 
•  nuaber  of  points  to  use  w.hen  drawing  U  arcs 


REAL  I,  U,  V 
REAL  X,  Y,  Z 
REAL  VSTEP 
REAL  USTEP 


scratch  variables 

current  calculated  x,y,i  values  (local) 
distance  between  each  V  segment  drawn 
angle  between  each  U  arc  point  drawn 


Deteraine  the  nuaber  of  points  to  plot  for  each  arc  and  the  step  angle 
between  each  point.  For  a  full  2  »  PI  (360  degree)  arc,  <8  points 
will  be  used.  The  ainiaua  nuaber  of  points  used  is  4. 


IF  (SUBDIV  ,Eq.  1)  THEN 
NUMUPOINTS  =  4  *  NSU(IRG) 

NUMVPOINTS  =  NSV (IRG) 
aSE 

I  =  24.0  *  (SUU(IRG)  -  SUL(IRG))  /  PI 
NUMUPOINTS  =  4  »  INT(I  /  40  .  .99) 
NUMVPOINTS  =  1 
END  IF 

VSTEP  =  (SVU(IRG)  -  SVL(IRG))  /  NUMVPOINTS 


74 


NRL  Memorandum  Report  6709 


USTEP  =  (‘^UUrifW)  -  SUL(IRG))  /  NUMUPQINTS 

V  =  SVL(IRG) 

BEGINPOINT  =  N  «  1 
DO  VCOUNT  =  t,  NUMVPOINTS 

C - 

C  Calculate  the  first  point  of  the  circle  on  the  top  of  cylinder 
C - - - 

X  =  CCS(1,IRG)  /  CCS(2,IRG)  ♦  (CCS(2,IRG)  -  V) 

♦  .  CQS(SUL(IRG)) 

Y  =  CCS(1,IRG)  /  CCS(2,IRG)  «  (CCS(2,IRG)  -  V) 

.  *  SIN(SUL(IRG)) 

2  =  V 

CALL  PUT  POINT(IRG,  N,  X,  Y,  2) 

ERDPOINT  =  N 

C - - 

C  Draw  the  rest  oT  the  circle  on  the  top  of  the  cylinder 
C - 

U  =  SUL(IRG)  .  'JSTEP 
DO  UCC'JNT  =  1,  NIMJHQINTS 

X  =  CCS(1,IRG)  /  CCS(2,IRG)  .  (CCS(2,IRG)  -  V) 

♦  ♦  CQS(U) 

Y  =  CCS(1,IRG)  /  CCS(2,IRG)  «  (CCS(2.IRG)  -  V) 

»  ?IN(U) 

CALL  PUT_POINT(IRG,  N,  X,  Y.  2) 

*RITE(CQN_UNIT,*)  N  -  1,  N,  COLOR 
U  =  U  ♦  USTEP 
END  DO 

V  =  V  .  VSTEP 
END  DO 

C . 

C  Draw  vertical  lines  to  siiulate  a  3D  cylinder,  then  exit. 

C . - . - . 

DO  UCOUNT  =  BEGINPOINT,  BEGINPOINT  .  NIMJPOINTS,  4 
IRITE(CON  UNIT,*)  UCOUNT,  ENDPOINT,  COLOR 
ENDPOINT  =  ENDPOINT  .  4 
END  DO 

RETURN 

END 
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c - - - 

c 

C  C0NNECT_ARR0WS{)  -  written  by  D.  Mitchell  Bolling  Jr. 

C  TRI/TESSCO,  Inc. 

C  P.O.  BOX  5684S8 

C  Orlando,  FL  32856 

C  January  1990 

C 

C  Purpose:  There  are  22  points  needed  to  draw  the  reference  axes, 

C  arrowheads,  and  the  axis  labels.  The  conne'-tivity  of  these  points 
C  are  added  to  the  connectivity  file  containing  the  CHIS-  figure 
r 

C  Subroutines  used:  NONE 
C 

C  Modification  Log: 

C 

C - 

SUBROUTINE  CONNECT_ARRO»S(CON_UNIT,  N) 

C - 

C  List  of  subroutine  paraeeters: 

C - -  - - 

INTEGER  CONUNIT  !  lun  nueber  for  connections  file 

INTEGER  N  !  keeps  track  of  the  next  point  to  be  created 

fRITE(CON_UNIT,»)  N  ♦  1,  N  •  2,  5  I  arrow  shafts  for  X,  Y,  I  axes 
•RITE(CON_UNIT,.)  N  *  1,  N  ♦  3,  5 

WRITE(CON_UNIT,»)  N  ♦  1,  N  .  4,  5 

fRITE(CON_UNIT,.)  N  .  5,  N  .  2,  5  '  X  arrowhead 

fRITE(CQN_UNIT,«)  N  ♦  2,  N  ♦  6,  5 

IRirE(CflN_UNIT,»)  N  ♦  6,  N  .  5,  5 

fRITE(C0N  UNIT,*)  N  •>  7,  N  .  3,  S  I  Y  arrowhead 

»RITE(C0N"UNIT,.)  N  .  3,  N  ♦  8,  5 

fRITE(CONlUNIT,«)  N  .  8,  N  ♦  7,  5 

WRITE(C0N  unit,*)  N  .  9,  N  .  4,  S  !  Z  arrowhead 

»RITE(CON_UNIT,*)  N  .  4,  N  ♦  10,  5 
»RITE(CQN_UNIT,*)  N  .  10,  N  ♦  9,  5 

IRITE(C0N  UNIT,*)  N  *  11,  N  .  12,  5  I  X  label 

»RITE(CON_UNIT,*)  N  ♦  13,  N  ♦  14,  5 

*RITE(CON_UNIT,*)  N  .  15,  N  -  16,  5  I  Y  label 

WRITE(CON_UNIT,*)  N  ♦  16,  N  .  17,  5 

WRITE(CON_UNIT,*)  N  *  18,  N  «  18,  S 

»RITE(CON_UNIT,*)  N  -  19,  N  .  20,  5  !  Z  label 

WRITE(CON_UNIT,*)  N  ■*  20,  N  .  21,  5 

WRITE(CON_UNIT,«)  N  ♦  21,  N  ♦  22,  5 

RETURN 

END 
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c — - - - - 

c 

C  R0T3D()  -  Iritten  by  Mitch  Bolling  of 
C  TRI/TESSCQ,  Inc. 

C  August  1989 

C 

C  R0T3D()  is  a  subroutine  that  inputs  3  dieensional  coordinates  and  3 
C  angles  of  rotation  and  returns  a  list  of  coordinates  that  have  been 
C  rotated  about  the  X,  Y,  and  Z  axes  by  the  aiount  specified  The 
C  rotations  are  done  by  3  separate  procedures.  The  first  to  rotate  arocd 
C  the  X  axis,  the  second  to  rotate  around  the  Y  axis  and  the  third  to 
C  rotate  around  the  Z  axis 

r 

c - 

SUeROUTINE  RaT3C(N,  AX,  AY,  AZ) 

PARAMETER  MAXCOORDS  =  10000 

C - 

C  List  of  parameters  used' 

C  - - 

INTEGER  N  I  number  of  coordinates  to  plot 

REAL  AX  !  angle  of  rotation  about  the  X  axis  (in  degrees) 

REAL  AY  !  angle  of  rotation  about  the  Y  axis  (in  degrees) 

REAL  AZ  !  angle  of  rotation  about  the  Z  axis  (in  degrees) 

REAL  A 
REAL  a 
REAL  C 
REAL  L 
REAL  V 

C . - . . 

C  List  of  global  variables  used: 

C . 

REAL  CDORD30  !  array  of  3d  points 

REAL  ROTATED  I  3d  points  ’fter  rotations 

COMMON  /  ROT_INFO  /  CaQRD30(3,  MAXCOORDS),  R0TAT£D(3, MAXCOORDS) 


C  List  of  local  variables  used: 

C 

C  Temporary  work  matrices  used  to  hold  rotated  homogeneous  coordinates. 

C . - .  - 

kEai.  u1(4,  MAXCOORDS) 

REAL  U2(4,  MAXCOORDS) 

REAL  TX(4,4)  I  describes  rotation  about  X  axis. 

REAL  TZ(4,4)  I  describes  rotation  about  Z  axis. 

REAL  RX(4,4)  !  describes  rotation  about  X  axis  until  axis  of 

I  rotation  is  m  the  XZ  plane. 

REAL  RXI(4,4)  '  describes  inverse  rotation  about  X  axis  until  axis 

I  of  rotation  is  back  in  original  position. 

REAL  RY(4,4)  i  describes  rotation  about  Y  axis  until  Z  axis 
I  corresponds  to  axis  of  rotation 

REAL  RYI(4,4)  !  describes  inverse  rotation  of  Y  axis 

REAL  TMP1(4,4)  !  temporary  work  matrices  used  to  ho'd  these  rotation 

REAL  TMP2(4,4)  !  matrices 

REAL  X  I  angle  of  rotation  about  the  X  axis  (in  radians) 

REAL  Y  !  angle  of  rotation  about  the  Y  axis  (in  radians) 

REAL  Z  !  angle  of  rotation  about  the  Z  axis  (in  radians) 

INTEGER  I,  J  I  scratch  counter  variables  (I  for  row,  J  for  column) 
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First  3  points  definad  in  U  will  represent  the  X,  Y,  end  I  sxes, 
respectively  snd  ire  used  to  keep  trick  of  the  orientition  of  these 
lies  IS  the  figure  is  rotited. 


Ul(l,  1)  =  1,« 
Ul(2,  1)  =  1.0 
Ul(3,  1)  =  0.0 
Ul(4,  1)  =  0.0 

Ul(l,  2)  =  0.0 
Ul(2,  2)  =  1.0 
Ul(3,  2)  >  0.0 
Ul(4,  2)  »  0.0 

Ul(l,  3)  =  0.0 
Ul(2.  3)  =  0.0 
Ul(3,  3)  =  1.0 
Ul(4,  3)  =  0.0 


Convert  the  coordinate  •>tri<!  into  a  hoiogeneous  latrix. 


DO  I  =  1,  N 
DO  J  =  1,  3 

U1(J,  I  *  3)  =  C00RD3D(J,  I) 
END  DO 

Ul(4,  I)  =  1.0 
END  DO 
N  =  N  '>  3 


Convert  rotation  ingles  free  degrees  into  radians 


X  =  -AX  /  57.20678 
Y  =  -AY  /  57.29578 
Z  =  -AZ  /  57.29678 


Define  the  4X4  iitrix  that  will  be  used  for  rotating  the  data 
about  the  X  axis 


TX(1,1)  =  1.0 
TX(2.1)  =  0.0 
TX(3,1)  =  0.0 
TX(4,1)  =  0.0 

TX(1,2)  •  0.0 
TX(2,2)  .  COS(X) 
TX(3,2)  »  -SIH(X) 
TX(4,2)  >  0.0 

TX(1,3)  =  0.0 
TX(2,3)  »  SIN(X) 
TX(3,3)  -  COS(X) 
TX(4,3)  =  0.0 

TX(1,4)  »  0.0 
TX(2,4)  =  0.0 
TX(3,4)  «  0.0 
TXC4,4)  w  1.0 


C . 

C  Wultiply  the  unrotited  coordinates  (Ul)  by  the  rotation  aatrix  (TX) 
C  to  produce  the  rotited  coordinates  (U2)  about  the  X  axis. 

C . 

DO  I  *  1,  N 
DO  J  =  1,  4 
U2(J.  I)  «  0.0 
DO  K  «  1,  4 

U2(J,  I)  =>  U2(J,  I)  *  U1(K,  1)  V  TX(K,  J) 
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END  00 
END  DO 
END  DO 

C - 

C  This  is  the  beginning  of  the  rotation  about  the  Y  axis  The 
C  coordinates,  after  a  rotation  about  the  X  axis,  are  stored  in  'J2 
C  The  variables  0,  C,  and  V  are  the  lengths  of  vectors  in  the  YZ  plane 
C  that  are  used  to  rotate  the  Y  axis  so  that  t  'S  a'igned  with  the 
C  Z  axis. 

C - 

8  =  U2(2,  2) 

C  =  U2(3,  2) 

V  =  Sq8T(B.8  *  C*C) 


C  Define  rotation  latrix  used  to  align  the  Y  axis  with  the  Z  axis. 


0X(1,1)  =  1  8 
RX(2,1)  =  0.0 
RX(3,i)  =  0.0 

RX(4,1;  =00 

RX(1,2)  =  0  0 
RX(2,2)  =  C  /  V 
RX(3,2)  =  -a  /  V 
RX(4,2)  =00 

RX(1,3)  =  0.0 
RX(2,3)  =  8  /  V 
RX(3,3)  =  C  /  V 
RX(4,3)  =  0.0 

RX(1,4)  =  0.0 
RX(2,4)  =  0.0 
RX(3,4)  =  00 
RX(4,4)  =  1  0 


f.  Define  the  4X4  aatrix  that  will  be  used  for  rotating  the  data 
C  about  the  Z  axis 

C . - .  - . - . 

:z(l,l)  =  CGS(Y) 

■z(2,ij  =  -s:n(y) 

'Z(3,l)  =  0  0 
■Z(4,l)  =00 

TZ(1,2)  =  SIN(Y) 

TZ(2,2}  =  CCS(Y) 

'Z(3,2)  =  00 
TZ(4,2)  =80 

TZ(l,3)  =  0  0 
^Z(2,3)  =  0.0 
TZ(3,3)  =  1.0 
TZ(4,3)  =  0.0 

TZfl,4)  =00 
TZ(2,4)  =00 
:Z(3,4)  =  0.0 
TZ(4,4)  =  1.0 


C  Th . s  4x4  aatrix  ;e(erses  the  rotation  about  the  X  ax'S 
C  (Inverse  of  RX)  This  p-ts  t^e  Y  axis  bach  to  its  previous  orientatio 


RXI(l,i;  =  I  0 
Rxirz.ii  =  0  0 
RX; (3, ;;  0  0 
Rxi(4,i;  =  0  0 

8X11,1,2)  -00 
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RXI(2,2)  =  C  /  V 
RXI(3,2)  =  8  /  V 
RXI(4,2)  =  8.8 

RXI(1,3)  =  8.8 
RXI(2,3)  =  -B  /  V 
RXI(3,3)  =  C  /  V 
RXI(4,3)  =  0  0 

RXI(1,4)  =  0.0 
RXI(2,4)  =  00 
RXI(3,4)  =  0.0 
RXI(4,4)  =  1.8 


C - 

C  To  siaplify  the  aatrix  iin i pu I  at  Ions,  first  Multiply  the  rotation 
C  Matrices  (RX,  T2  and  RXI)  together.  Then  Multiply  the  coordinates  (U2) 
C  with  the  coebined  product  (TUP2)  U1  =  U2  •  (RX  ♦  TZ  •  RXI) 

C 

C  Mu  1 1  i  p  ly  RX  by  TZ 

C - 

DO  I  =  1,  4 
DO  J  =  1,  4 

TMP1(I,J)  =  0.0 
DO  K  =  1,  4 

TMP1(I,J)  =  TMP1(I,J)  .  RX(I,K)  .  TZ(K,J) 

END  DO 
END  DO 
END  DO 


C . . . . 

C  Multiply  the  resulting  product  by  RXI  to  produce  TMP2 
C  (RX  .  TZ  *  RXI) 


DO  I  =  1,  4 
00  J  =  1,  4 

TMP2(I,J)  =  0.0 
DO  K  =  1,  4 

TMP2(I,J)  =  TMP2(I,J)  ♦  TUP1(I,K)  .  RXI(K,J) 
END  00 
END  DO 
END  DO 


C . - . . . 

C  Multiply  the  coordinates  (U2)  that  have  been  rotated  about  the  X  axis 
C  by  the  rotation  Matrix  (TMP2)  to  produce  coordinates  that  have  been 
C  rotated  about  the  Y  axis  (UI) . 

C - - - - - 

DO  I  =  1,  N 
DO  J  =  1,  4 

U1(J,  I)  =  0.0 
DO  K  =  1,  4 

U1(J,  I)  =  U1(J,  I)  •  U2(K,  I)  .  TUP2(K,  J) 

END  DO 
END  DO 
END  DO 


C - - - 

C  ''his  is  the  beginning  of  the  rotation  about  the  Z  axis  The 

C  coordinate  after  a  rotation  about  the  Y  axis  is  stored  in  Ul  The 

C  xariables  8,  C,  and  V  are  the  lengths  of  vectors  in  the  YZ  plane 
C  that  are  used  to  rotate  the  Y  axis  so  that  it  is  aligned  with  the  Z 
C  axis.  The  variables  V,  A,  and  L  are  lengths  of  vectors  m  the  XZ 

C  plane  used  to  rotate  the  Z  axis 

C - - - - - 

A  =  Ul(l,  3) 

B  =  Ul(2,  3) 

C  =  Ul(3,  3) 

¥  =  SQRTfB.9  •  C.C) 

L  =  S«RT(A.A  .  BtS  •  C.C) 
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C  Define  rotation  the  aatrices  used  to  align  the  Z  axis  with  the  global 
C  Z  axis  before  rotation. 


RX(1,1)  = 

10 

RX(2,1)  = 

0.0 

RX(3,1)  = 

0  0 

RX(4,1)  = 

0  0 

RX(1,2)  = 

0  0 

Ra(2,2)  = 

c  / 

V 

RX(3,2)  = 

-B  / 

V 

RX(4,2)  = 

0  0 

RX(1,3)  = 

0  0 

RX(2,3)  = 

8  / 

V 

RX(3,3)  = 

c  / 

V 

RX(4,3)  = 

0  0 

RX(1,4)  = 

0  0 

RX(2,4)  ^ 

0  0 

RX(3,4)  = 

0  0 

RX(4,4)  = 

1  0 

RY(1,1)  = 

V  / 

L 

RY(2,1)  = 

0  0 

RY(3,:)  = 

-A  / 

L 

RY(4,1)  = 

0  0 

RY(1,2)  = 

0  0 

RY(2,2)  = 

1.0 

RY(3,2)  = 

0  0 

RY(4,2}  = 

0  0 

RY{1,3)  = 

A  / 

L 

RY(2,3)  = 

0  0 

RY(3,3)  = 

V  / 

L 

RY(4,3)  = 

0  0 

RY(1,4)  = 

0  0 

RY(2,4)  = 

0  0 

RY(3,4)  = 

0  0 

RY(4,4]  = 

10 

[)cf  ( ne  the  4  X  - 

♦  aatrix  that  will  be  used  for  rotating  the  data 

about  the  Z  ax  is 

TZ(1,1)  =  COS(Z) 
TZ(2,1)  -  -SIN(Z) 
TZ(3,1)  =  0  0 
TZ(4,1)  =:  0  0 

TZ(1,2)  =  SIN(Z) 

rz(2,2)  =  cas(Z) 

TZ(3,2)  =  0  0 
TZ(4,2)  =  0  0 

TZ(1,3)  =  00 
TZ(2,3)  =  00 
TZ(3,3)  =  1  0 
TZ(4,3)  =  00 

TZ(1,4)  =  0  0 
TZ(2,4)  =  0.0 
TZ(3,4)  =  0  0 
TZ(4,4)  =  1  0 


C  - . - . - . . . . . . 

C  'base  2  inverse  rotation  aatr  cos  are  used  to  wove  me  Z  a«  s  Oacx 
C  to  its  original  position  (Inverse  o'  SX  and  RY) 
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RXI(1,1)  r  10 
RXI(2,1)  =  0.0 
RXI(3,1)  =  0.0 
RXI(4,1)  =  0.0 

RXI(1,2)  =  0.0 
RXI(2,2)  =  C  /  V 
RXI(3,2)  =  B  /  V 
RXI(4,2)  =  10 

RXI(1,3)  =  0.0 
RXI(2,3j  =  -B  /  V 
RXI(3,3)  =  C  /  V 
RXI(4,3)  =  0.0 

RXI(1,4)  =  00 
RXI(2,4)  =00 
RXI(3,4)  =  00 
RXI(4.4)  =10 


C - - - 

C  Inverse  of  RY. 

C-- - - - - 

RYI(1,1)  =  V  /  L 
RYI(2,1)  =  0.0 
RYI(3,1)  =  A  /  L 
RYI(4.1)  =  0.0 

RYI(1,2)  =  0.0 
RYI(2,2)  =  10 
RYI(3,2)  =  0.0 
RYI(4,2)  =  00 

RYI(1,3)  =  -A  /  L 
RYI(2,3)  =  0  0 
RYI(3,3)  =  V  /  L 
RYI(4,3)  =  0.0 

RYI(1,4)  =  0.0 
RYi(2,4)  =00 
RYI(3,4)  =  00 
RYI(4,4)  =  1.0 


C- . - . 

C  First  eultiply  the  rotation  eatrices  (RX,  RY,  TZ,  RYI,  and  RXI)  together. 
C  Then  eultlply  the  coordinates  (U2)  with  the  coebined  product  (TMP2) 

C  Multiply  RX  by  RY 

C - - - - - 

DO  I  =  1,  4 
DO  J  =  1,  4 

TMP1{I,J)  =  00 
DO  K  =  1,  4 

TMP1(I,J)  =  TMP1(I,J)  .  RX(I,K)  •  RY(K.J1 
END  DO 
END  DO 
END  DO 


- - 

C  Multiply  the  resulting  product  by  TZ  to  produce  TMP2  (RX  .  RY  •  TZ) 


DO  I  =  1,  4 
DO  J  =  1,  4 

TMP2(I,J)  =00 
DO  X  =  1,  4 

TMP2(I,J)  =  TMP2(I,J)  .  TMP1(I,K)  •  TZ(K,J) 
END  00 
END  00 
END  DO 
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C  Multiply  the  resulting  product  by  RYI  to  produce  TUPl 
C  (RX  ♦  RY  »  TZ  ♦  RYI) 

C - 

DO  I  =  1,  4 
DO  J  =  1,  4 

TMP1(I,J)  =  0* 

DO  K  =  1,  4 

TMP1(I,J)  =  TMP1(I,J)  •  TUP2(I,K)  *  RYI(K,J) 
END  DO 
END  DO 
END  DO 


C 

C 

C 

C 


Multiply  the  resulting  product  by  RXI  to  produce  TMP2 
(RX  *  RY  «  TZ  «  RYI  .  RXI) 


DO  I  =  1,  4 
DO  J  -  1,  4 

TMP2(I,J)  =00 
DO  K  =  1,  4 

TMP2(I,J)  =  rMP2(I,j)  .  TMP1(I,K)  .  RXI(K,J) 
END  DO 
END  DO 
END  DO 


C . - . . . - . 

C  Multiply  the  coordinates  (UI)  rotated  about  the  X  and  Y  aiies  by 
C  the  rotation  eatrix  (TMP2)  to  produce  coordinates  that  have  been 
C  rotated  about  the  Z  axis  (U2) . 

C . . . . 

DO  I  =  1,  N 
DO  J  =  1,  4 

U2(J,  I)  =  0.0 
DO  K  =  1,  4 

U2(J,  I)  =  U2(J,  I)  ‘  U1(K,  I)  .  TMP2(K,  J) 

END  DO 
END  DO 
END  DO 

C .  - . . . - . 

C  Output  the  latrix  of  rotated  X,  Y,  and  Z  coordinaves,  oeittlng  the 
C  hoiogeneous  coordinates. 

C .  . - . ---- . - . 

N  =  N  -  3 
DC  I  =  I,  N 
DO  J  =  1,  3 

ROTATED(J,  I)  =  U2(J,  I  -  3) 

END  DO 
END  DO 

RETURN 

END 
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c - 

c 

C  DRA»_ARROWS()  -  Written  by  Mitch  Bolling  of 
C  TRI/TESSCO,  Inc. 

C  January  1990 

C 

C  DRAW_ARROWS()  is  a  subroutine  that  generates  the  points  needed  for 
C  drawing  the  arrow  heads  and  letters  (X,  Y,  Z)  of  the  reference  axes. 
- - 


SUBROUTINE  DRAW^ARROWS (CRD  UNIT.  N.  AX,  AY,  AZ) 


PARAMETER  MAXCOORDS  =  10000 

- - 

C  List  of  paraeeters  used: 

- - 

INTEGER  CRDUNIT  I  lun  nueber  for  coordinate  file 

INTEGER  N  !  number  of  coordinates  currently  being  plotted 


REAL  AX 
REAL  AY 
REAL  AZ 


I  angle  of  rotation  about  the  X  axis  (in  degrees) 
!  angle  of  rotation  about  the  Y  axis  (in  degrees) 
!  angle  of  rotation  about  the  Z  axis  (in  degrees) 


C . - . - . - . 

C  List  of  global  variables  used: 

C . - . - . 

REAL  COOR03D  '  an  ay  of  unrotated  3d  points 

REAL  ROTATED  I  3c  points  after  rotations 

CONMON  /  R0T_INFQ  /  CaOR030(3,  MAXCOORDS),  ROTATED(3,MA:XOORDS) 


C  List  of  local  variables  used; 

C . . - . 

INTEGER  I,  J  I  scratch  counter  variables  (I  for  row,  J  for  column) 


REAL 

ARROWS (3,  4) 

!  4  points  that  define  the  reference  axes 

REAL 

MINX 

!  leftmost  point  of  rotated  object 

REAL 

MINY 

!  bottoe  point  of  rotated  object 

REAL 

MINZ 

1  deepest  point  of  rotated  object 

REAL 

MAXX 

I  rightmost  point  of  rotated  object 

REAL 

MAXY 

1  top  point  of  rotated  object 

REAL 

MAXZ 

!  highest  point  of  rotated  object 

REAL 

MAXDIFF 

1  length  of  side  of  box  drawn  around  CHIEF 

REAL 

MAXRFSIZE 

!  length  of  box  drawn  around  reference  axes 

REAL 

XCENTER 

!  global  x  position  for  tip  of  arrowhead 

REAL 

YCENTER 

!  global  y  position  for  tip  of  arrowhead 

REAL 

ZCENTER 

!  global  i  position  for  tip  of  arrowhead 

REAL 

HEADX 

1  X  offset  for  arrowhead  base  or  axis  label 

REAL 

HEADY 

1  y  offset  for  arrowhead  base  or  axis  label 

REAL 

ARROWLEN 

!  length  of  XY  vector  of  an  arrow 

DATA  ARROWS  /  0  0, 

0.0,  00,  !  center 

♦ 

0  8S, 

00,  00,  !  X  shaft 

* 

0  0, 

0  85,  00,  1  y  shaft 

♦ 

0  0, 

00,  0  85  /  '  I  shaft 

Sea  1 ing 

the  reference  axis  witn  respect  to  rotated  view  of  chief 

figure. 

MINX 

=  ROTATEDd,  1) 

MINY 

=  R0TATED(2,  1) 

MINZ 

=  R0TATED(3,  1) 

MAXX 

^  MINX 

MAXY 

=  MINY 

MAXZ 

=  MINZ 

DO  I 

=  2,  N 

IF  (ROTATEDd,  I) 

LT  MINX)  MINX  =  R0TATFD(1,  I) 
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IF  (R0TATED(2,  I)  LT. 
IF  (R0TATED(3,  I)  LT. 
IF  (R0TATED(1,  I)  .GT. 
IF  (R0TATED(2,  I)  GT. 
IF  (R0TATED(3,  I)  GT. 
END  OQ 


MINT)  MINT  =  R0TATED(2,  I) 

MINZ)  MINZ  =  RQTATED(3,  I) 

MAXX)  MAXX  =  ROTATED  0,  I) 

MAXY)  MAXY  =  ROTATED^Z.  I) 

MAXZ)  MAXZ  =  R0TATEC(3.  I) 


C - 

C  Find  the  largest  difference  of  rotated  coordinates 

C - 

MAXDIFF  =  MAXX  -  MINX 

IF  ((MAXY  -  MINY)  GT  MAXDIFF)  MAXDIFF  =  MAXY  -  MINY 
IF  ((MAXZ  -  MINZ)  GT  MAXDIFF)  MAXDIFF  =  MAXZ  -  MINZ 


C - - - 

C  sake  scale  of  axes  l/6th  of  size  of  CHIff  figure. 

C- - - 

MAXREFSIZE  =  MAXCIFF  /  6 

N  =  A 

DO  I  =  1,  N 

C0QRD3D(1,  I)  =  ARR0»S(1,  I)  .  MAXRFSIZE 

C00RD3D(2,  I)  =  ARR0fS(2,  I)  *  MAXREFSIZE 

COORD30(3,  I)  =  ARR0»S(3,  I)  ♦  MAXREFSIZE 

END  DO 


Rotate  reference  axes  in  the  sale  way  as  the  object  was 

rotated . 

CAU  R0T3D(N,  AX,  AY,  AZ) 

Now  write  rotated  reference  axes  in  bottoi  left  corner, 
points  will  becoie  the  tips  of  the  reference  arrows. 

These 

DC  I  =  I,  N 

fRITE(CRO  UNIT,  .)  ROTATED(I,  I)  ♦  MINX  -  MAXREFSIZE, 
.  R0TATED(2,  I)  ■>  MINY  -  MAXRFSIZE, 
.  R0TATED(3,  I)  .  MINZ  -  MAXRffSIZE 
END  DO 


C . . - . --- 

C  Create  points  used  to  draw  the  base  of  arrowhead  for  each  axis,  one 
C  point  on  each  side  of  the  shaft.  These  points  and  the  arrow  tip 
C  wi I  I  be  connected  in  a  triangle  to  fore  the  arrowhead. 

C . .  . - . - .  - 

DO  I  =  2,  4 

ARRO»LEN  =  SfiRT(RO'"ATED(l,I)**Z  .  R0TATE0(2,I)**2) 


C - - - 

C  First,  locate  tip  of  arrows  on  reference  axis  in  global  coordinates 
C - 

XCENfER  =  R0TATED(1,I)  ♦  MINX  -  MAXREFSIZE 
YCENTER  =  R0TATED(2,I)  -  MINY  -  MAXREFSIZE 
ZCENTER  =  R0TATED(3,T)  *  MINZ  -  M/XREFSIZE 

C--  - - - 

C  Decide  when  reference  axis  vector  is  scalier  than  arrowhead.  If  so 
C  then  write  the  two  arrowhead  base  points  on  top  of  the  arrow  tip 
C  (In  this  way,  the  arrowhead  will  not  appear  when  an  axis  is  seen 
C  a  I  lost  head  on.) 

C . - . . . - - - 

IF  (ARROfLEN  LT.  MAXRffSIZE  /  3)  THEN 

fRITE(CRD_UNIT,  .)  XCENTER,  YCENTER,  ZCENTER 
fRITE(CRD_UNIT,  «)  XCENTER,  YCENTER,  ZCENTER 
ELSE 

C - 

C  Dividing  by  ARROfLEN  lakes  it  a  unit  vector.  H-jIliply  by  MAXREFSIZE 
C  to  deteriine  how  large  arrowhead  will  be  froi  tip  to  base,  (this  way, 
C  the  arrowheads  for  all  three  axes  will  be  of  the  sale  length  ) 

C - - - - - 

HEADX  =  MAXREFSIZE/4  •  R0TATED(1,I)  /  ARROfLEN 
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HEADY  =  MAXREFSIZE/4  «  R0TATEIi(2, 1)  /  ARROWLEN 


C - - - 

C  Create  the  two  points  of  the  base  of  the  arrowhead. 

C - 

fRITE(CRD_UNIT,  »)  XCENTER  -  HEADX  ♦  HEADY/4, 

.  YCEHTEh  -  HE>DY  -  HEADX/4, 

.  ZCENTFR 

WRITE(CRD_UNIT.  ♦)  XCENTER  -  HEADX  -  HEADY/4, 

.  YCENTER  -  HEADY  .  HEAOX/4, 

♦  ZCENTER 

END  IF 
END  Ui, 


C - 

C  Create  points  used  to  draw  the  letter  ’X’. 

C - 

ARROWLEN  =  SqRT(RQTATED(l,2)M2  *  ROTATED (2. 2) »* 2) 
IF  (ARROWLEN  .£!(.  0.8)  ARROWLEN  =1.0 

XCENTER  =  R0TATED(1,2)  ♦  MINX  -  MAXREFSIZE 
YCENTER  =  ROTATED (2, 2)  •  MINY  -  MAXRffSlZE 
ZCENTER  =  R0TATED(3,2)  >  MINZ  -  MAXRffSIZE 


C . . . . . - . 

C  Move  the  saee  distance  away  frow  the  tip  of  arrow  as  the  arrowhead 
C  base  is,  but  love  in  the  opposite  direction. 

C . - . - . 

HEADX  =  MAXRffSIZE/4  .  RaTATED(l,2)  /  ARROWLEN 
HEADY  =  MAXRff SIZE/4  •  ROTATED (2, 2)  /  ARROWLEN 


C . 

C  Write  thj  points  describing  and  'X'  to  the  coordinate  file. 

Q  —  —  —  •  —  — 

WRITE(CRD  UNIT,  •)  XCENTER  ♦  HEADX  .  MAXREFSIZE/ J 2, 

*  '  YCENTER  ♦  HEADY  «  MAXREFSIZE/12, 

‘  ZCENTER 


WRITE(CRD  UNIT,  .)  XCEN.ER  ♦  HEADX  -  MAXREFSIZE/12, 
YCENTER  .  HEADY  -  MAXREFSIZE/ 12, 
ZCENTER 


WRITE(CRD  UNIT,  *)  XCENTER  ♦  HEADX  -  MAXREFSIZE/12, 
YCENTER  ♦  HEADY  .  UAXRB^SIZE/1?, 
ZCENTER 


WRITE(CRD_UNIT,  *)  XCENTER  •  HEADX  ♦  MAXREFSIZE/12, 
YCENTER  ♦  HEADY  -  MAXREFSIZE/12, 
ZCENTER 


C . - - - 

C  Create  points  used  to  draw  the  letter  ’Y’ . 

C-- . . . . . . 

ARROWLEN  =  SqRT(R0TATED(l,3)M2  .  R0TATED(2,3)»»2) 
IF  (ARROWLEN  Efl.  0.0)  ARROWLEN  =1.0 


C - - - - - 

C  Move  the  saee  distance  away  froe  the  tip  of  arrow  as  the  arrowhead 
C  base  is,  but  sove  in  the  opposite  direction. 

C - 

XCENTER  =  R0TATED(1,3)  ♦  MINX  -  MAXREFSIZE 
YCENTER  =  R0TATED(2,3)  •  MINY  -  MAXREFSIZE 
ZCENTER  =  ROTATED (3, 3)  .  MINZ  -  MAXRFSIZE 

HEADX  =  MAXRffSIZE/4  «  R0TATED(1,3)  /  ARROWLEN 
HEADY  =  MAXRff SIZE/4  »  ROTATED (2, 3)  /  ARROWLEN 


C -  - 

C  Write  the  poin.s  describing  and  'Y'  to  the  coordinate  file 

C - - -  - - - 

WRITE(CRD_UNIT,  *)  XCENTER  ♦  HEADX  -  MAXREFSIZE/12, 
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YCENTER  ♦  HEADY  «  «AXRffSIZE/12, 
ZCENTER 

WRITE(CRD_UNIT,  »)  XCENTER  »  HEAOX, 

YCENTER  -  HEADY, 

ZCENTER 

WRITE(CRD_UNIT,  «)  XCENTER  ♦  HEADX  -  UAXREFSIZE/12, 
YCENTER  .  HEADY  .  NAXREFSIZE/12. 
ZCENTER 

fRITE(CRD_UNIT,  »)  XCENTER  -  HEADX, 

►  YCENTER  .  HEADY  -  kAXRB^SIZE/12, 

.  ZCENTER 


Create  points  used  to  draw  the  letter  ’Z’. 


ARROWLEN  =  SURTCRC '\TED(1,4).«2  ♦  RQTATED(2,4)»*2) 

IF  (ARROWLEN  .ER.  B.8)  ARROWLEN  =1.8 

XCENTER  =  RCTATED(1,4)  ♦  MINX  -  UAXREFSIZE 
YCENTER  =  R0TATED(2,4)  .  WINY  -  MAXRB-SIZE 
ZCENTER  =  ROTATED (3, 4)  .  MINZ  -  MAXRffSIZE 

C . . . - . - . 

C  Move  the  saee  distance  away  froe  the  tip  of  arrow  as  the  arrowhead 
C  base  is,  but  eove  in  the  opposite  direction. 

C . . . 

HEADX  =  MAXRBFSIZE/4  *  R0TATED(1,4)  /  ARRUWLEN 
HEADY  =  UAXRff SIZE/4  «  ROTATED (2, 4)  /  ARROWLEN 

C . . 

C  Write  the  points  describing  and  'Z'  bo  the  coordinate  file. 

C . - . - . 

WRITE(CRO  UNIT,  »)  XCENTER  *  HEAOX  -  UAXRffSIZE/12, 

♦  ■  YCENTER  ♦  HEADY  ♦  MAXRffSIZE/12, 

ZCENTER 

WRITE(CRD  UNIT,  XCENTER  ♦  HEAOX  -  MAXRffSIZE/12, 

-  YCENTER  *  HEADY  •>  MAXRffSIZE/12, 

♦  ZCENTER 

WRITE(CRD  UNIT,  XCENTER  ♦  HEADX  -  MAXhJ^SIZE/lZ, 

.  YCENTER  .  HEADY  -  MAXREFSIZE/12, 

♦  ZCENTER 

WRITE(CRU_UNIT,  »)  XCENTER  ♦  HEADX  ♦  MAXRffSIZE/12, 

♦  YCENTER  .  HEADY  -  MAXRffSIZE/12, 

«  ZCENTER 

RETURN 

END 
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